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.his  report  describes  the  development  of  a finite-element  model 
for  analyzing  sheet-metal  forming  processes.  Materials  are  assumed 
to  be  rigid-plastis  with  the  view  that  the  usefulness  of  an  analysis 
method  depends  largely  upon  solution  accuracy  and  computation  effi- 
ciency. First,  the  variational  formulation  applicable  to  sheet-metal 
forming  is  described  by  considering  solution  uniqueness  and  the 
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effect  of  geometry  change  involved  in  the  forming  processes.  From 
this  variational  formulation,  a finite-element  process  model  based 
on  the  membrane  theory  is  developed.  Then,  three  basic  sheet-metal 
forming  processes,  namely,  the  bulging  of  a sheet  subject  to  hydro- 
static pressure,  the  stretching  of  a sheet  with  a hemispherical  head 
punch,  and  deep  drawing  of  a sheet  with  a hemispherical  head  punch, 
are  solved.  The  solutions  arrived  at  by  the  rigid-plastic,  finite- 
element  method  are  compared  with  existing  numerical  solutions  and 
the  experimental  data.  The  agreement  is  generally  excellent  and  it 
is  concluded  that  the  rigid-plastic,  finite-element  method  is  effi- 
cient for  analyzing  sheet-metal  forming  problems  with  reasonable 
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SECTION  I 


INTRODUCTION 

The  metal  forming  processes  basically  involve  large  amounts  of  elastic 
deformation,  and,  due  to  the  complexities  of  plasticity,  the  exact  analysis 
of  a process  is  infeasible  in  most  of  the  cases.  Thus,  a number  of  approxi- 
mate methods  have  been  suggested,  with  varying  degrees  of  approximation  and 
idealization.  Among  these,  techniques  using  the  finite-element  method  take 
precedence  because  of  their  flexibility,  ability  to  obtain  a detailed  solution, 
and  the  inherent  proximity  of  their  solutions  to  the  exact  one. 

A prime  objective  of  mathematical  analysis  of  metalworking  processes  is 
to  provide  necessary  information  for  proper  design  and  control  of  these 
processes.  Therefore,  the  method  of  analysis  must  be  capable  of  determining 
the  effects  of  various  parameters  on  metal  flow  characteristics.  Furthermore, 
the  computation  efficiency,  as  well  as  solution  accuracy,  is  an  important 
consideration  for  the  method  to  be  useful  in  analyzing  metalworking  problems. 

With  this  viewpoint  in  mind,  successful  efforts  have  been  carried  out 
in  analyzing  various  deformation  processes,  such  as  compression,  heading, 
piercing,  extrusion  and  drawing  by  the  rigid-plastic,  finite-element  method 
(matrix  method)  [1]  — [7] . 

The  formulation  of  the  matrix  method,  however,  cannot  be  extended  to 
the  sheet-metal  forming  analysis  due  to  the  following  reasons: 

(1)  The  classical  variational  formulation  which  is  the  basis  of  the 
matrix  method  does  not  necessarily  determine  a unique  deformation 
mode.  Physically,  there  is  no  inherent  indeterminacy  for  work- 
hardening solids,  but  this  indeterminacy  is  due  rather  to  the  fact 
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that  the  workhardening  rate  is  not  included  in  the  mathematical 
formulation  of  the  classical  variational  principle. 

(2)  The  kinematic  assumption  in  the  matrix  method  is  not  longer  valid 
for  the  sheet-metal  forming  process.  As  long  as  bulk  deformation 
or  in-plane  stretching  are  concerned,  this  kinematic  assumption 
that  the  magnitude  of  the  rate  of  rotation  is  negligible  compared 
to  the  strain  rate  does  not  deviate  much  from  the  real  situation 
and  yields  solutions  consistent  with  reality.  Geometric  nonlinear- 
ity in  sheet-metal  forming,  however,  invalidates  such  a simplification 
The  objective  of  the  present  investigation  is,  therefore,  to  develop  and 
establish  a finite-element  method  for  sheet-metal  forming  processes. 

In  Section  II  various  forms  of  variational  formulations  are  reviewed 
in  the  light  of  uniqueness  and  geometry  change  which  leads  to  a realization 
of  the  necessity  of  new  formulations.  In  Section  III  a new  formulation  is 
obtained  and  the  development  of  the  finite-element  model  from  it  is  described. 
With  the  particular  example  of  sheet-metal  forming  processes  in  mind,  the 
idealization  of  plane  stress  state  and  membrane  theory  is  implemented. 
Furthermore,  the  development  is  confined  to  the  case  of  axisymmetrical 
problems . 

To  establish  the  validity  of  the  proposed  method,  three  basic  sheet- 
metal  forming  processes  are  analyzed  and  the  solutions  are  compared  with 
other  available  experimental  data  and  numerical  solutions.  Hydrostatic 
bulging  is  treated  in  Section  IV.  Punch  stretching  with  a hemispherical 
punch  is  discussed  in  Section  V.  To  make  the  problem  tractable,  one  moving 
contact  boundary  is  considered  first  by  neglecting  die  profile;  then  the 
analysis  is  extended  to  include  two  moving  boundaries.  In  Section  VI  deep 
drawing  with  a hemispherical  punch  is  solved. 
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SECTION  II 


BACKGROUND 


1 . Uniqueness 

We  consider  the  quasistatic  deformation  of  a rigid-plastic  solid. 

On  a portion  Sy  of  the  surface  S of  this  body  are  prescribed  given  velocities, 
while  the  remainder  S^.  of  the  surface  S is  subjected  to  given  surface  trac- 
tions T^.  Assuming  that  these  surface  velocities  and  tractions  are  such 
that  the  entire  body  is  in  a state  of  plastic  flow,  we  want  to  determine 
the  stresses  CL  ^ and  strain  rates  throughout  the  body. 

The  conventional  formulation  of  variational  principle  for  this  problem 
is  that  among  all  kinematically  admissible  strain  rate  fields  e?L,  the  actual 
one  minimizes  the  expression  (Hill  [8]), 


itj  - f ae* 


dv  - 


T.v*  dS, 
1 l ’ 


CD 


where  a is  the  effective  stress,  e is  the  effective  strain  rate  defined  by 
5 = f \ /o!  .o’  . , 

2 ij  i]’ 

j .Ajtt:. 

3 ij  ij* 


respectively,  where  a!  ^ is  the  deviatoric  component  of  cl  ^ . Here  a strain 
rate  field  e?L , defined  throughout  the  body  under  consideration,  is  called 
kinematically  admissible  if  it  is  derivable  from  a velocity  field  vt  which 
satisfies  the  condition  of  incompressibility  v?  ^ = 0^  throughout  the  body 


The  comma  denotes  the  differentiation  with  respect  to  coordinates,  e.g., 
9yi 

Vi, i - 9x.  ' 
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and  the  boundary  conditions  on  Sy.  The  variational  principle  in  this  form 
has  been  successfully  ap*’ied  to  the  analysis  of  metal  forming  problems, 
such  as  extrusion  [6].  As  was  found  out  later,  and  we  will  discuss  this 
shortly,  the  success  is  related  to  the  type  of  boundary  conditions  prescribed 
on  the  surface  of  the  body  undergoing  deformation.  In  general,  with  the 
variational  formulation  of  tt^  in  Eq . (1),  there  is  a question  regarding 
uniqueness  of  deformation  mode  even  though  the  stress  field  is  uniquely 
determined  [8],  [9]. 

Consider  an  incipient  flow  in  a rigid-plastic  solid,  workhardening  or 
perfect  plastic,  governed  by  the  following  partial  differential  equations 
which  are,  of  course,  dual  to  the  variational  formulation  . With  respect 
to  Cartesian  reference  frame  the  following  equations  hold: 


Equilibrium  equations 

a. . . = 0 in  the  absence  of  body  force 
i],3  7 


(2a) 


Strain  rate-velocity  relationship 

e.  . = "x“(v . . + v.  .) 
ij  2"  i,3  3,i 


Constitutive  equation 


pa!j  = p being  an  arbitrary  constant 


(2b) 


(2c) 


Yield  criterion 

a = f\jo .'  . a . = H(e), 
2 xj  ij 


where  e is  the  effective  strain  defined 


by  £ = 


de  if  de 


4 /de. .de. . 
3 ij  ij 


(2d) 


Boundary  conditions 

n . a . . = T . on  ST, 

3 iJ  i T 

v.  = v.  on  S , 

1 l v 


(2e) 

where  n.  is  the  unit  normal  vector  to  the 
3 

surface  of  the  body;  and  v^  are  prescribed 
values 
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Suppose  that  is  the  solution  to  this  boundary  value  problem. 


Construct  a different  set  of  stress  fields  and  strain  rate  fields 

ij  iJ 


where  a C is  any  arbitrary  factor  and  may  vary 
ij  iJ  ij  iJ 


from  point  to  point  throughout  the  body.  Then,  it  is  easily  shown  that  this 


f 'y  # 

set  (cKj  ,eij  ^ satisfies  all  the  governing  equations  except  for  the  boundary 


:(2) 


conditions  on  Sy.  On  Sy  the  velocity  integrated  from  should  coincide 


with  the  prescribed  value  v^.  Since  strain  rate-velocity  relation  is  linear, 


integrating  would  yield  Cv^  if  is  integrated  to  give  \r,  and  there- 


fore C must  be  unity  on  S . With  tliis  and  the  compatibility  requirement  the 


deformation  mode  may  or  may  not  be  uniquely  determined.  One  example  of  a 
well-established  unique  kinematic  mode  is  in  the  plane-strain  problem.  In 
the  plane-strain  condition,  unless  one  family  of  the  characteristics  is 
straight,  the  governing  equation  of  the  velocity  field  becomes  the  telegraphy 
equation  which  is  hyperbolic  and,  therefore,  the  solution  is  uniquely  deter- 
mined if  the  boundary  curve  is  not  along  a characteristic. 


It  can  be  readily  shown  that  under  certain  boundary  conditions  the 


set  (o^,Ce?j^)  also  satisfies  the  boundary  conditions  on  Sy  and  therefore 


the  deformation  mode  is  clearly  not  unique.  The  following  is  a partial  list 
of  such  boundary  conditions. 


(1)  Sv  = 0,  i.e.,  all  the  boundaries  are  traction  boundaries; 


(2)  vi  = 0 on  Sy, 


(3)  On  S only  the  ratio  between  the  velocity  components  are  prescribed, 
v. 

e.g.,  ^ = a; 

j 


(4)  Mixed  boundary  condition;  e.g.,  a normal  component  of  v.  and  a tan- 

A 

gential  component  of  T^  are  prescribed  over  the  surface,  or  vice 
versa.  In  this  case,  the  additional  condition  of  whether  all  the 
characteristics  meet  on  a curve  in  the  region  should  be  checked  [10] 
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Concrete  examples  are  (1)  the  expansion  of  spherical  shells  [11]  or  cylindri- 
cal shells  [12]  under  internal  pressure,  and  (4)  the  indentation  of  a semi- 
infinite body  by  a flat  punch  under  the  plane-strain  condition  [13],  torsion 
of  a prismatic  bar  [10].  Among  sheet-metal  forming  processes,  hydrostatic 
bulging  belongs  to  case  (2)  and  punch  stretching  to  case  (4)  or  (3) . 

Note  that  the  physical  meaning  of  these  boundary  conditions  is  that 
the  plastic  flow  is  unconstrained  and  all  or  part  of  the  body  is  free  to 
deform.  Mathematically,  this  nonuniqueness  is  due  to  the  fact  that  the 
Levy-Mises  theory,  implied  in  the  variational  formulation  tt ^ and  also 
appearing  in  the  differential  equations  (2c),  does  not  include  the 
"viscosity  effect"  (in  Prager's  terminology  [9])  and,  therefore,  this 
indeterminacy  would  be  resolved  if  the  workhardening  effect  is  taken  into 
account.  In  fact,  for  the  workhardening  solid  there  is  no  inherent  indeter- 
minacy in  general;  the  apparent  nonuniqueness  is  due  simply  to  an  inadequate 
formulation  of  the  problem.  In  proper  formulation,  traction  rate  T\  must  be 
specified  on  S^.,  and  then  from  an  infinite  number  of  kinematically  possible 
modes  the  actual  mode  can  be  singled  out  by  the  additional  requirement  that 
there  must  exist  an  equilibrium  distribution  of  stress  rate  compatible  with 
the  implied  rate  of  hardening  everywhere  in  the  body  and  with  the  given 
traction  rate  on  S^,.  Besides,  the  workhardening  effect  is  explicitly 
brought  into  the  constitutive  equation  in  the  form  of 


a!  . 


he . . = — a (3) 

1J  5 

where  a is  the  time  rate  of  a,  h the  workhardening  effect  of  the  material 

being  equal  to  — . It  can  be  shown  that  the  constitute  equation  (3) 
de 

can  always  be  reduced  to  the  constitutive  equation  (2c),  but  not  necessari 
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vice  versa.  Therefore,  for  a perfectly  plastic  solid,  specifying  the  trac- 
tion rate  does  not  resolve  the  indeterminacy.  Hill,  then,  showed  that  among 
all  variational  modes  compatible  with  the  boundary  conditions  for  v..  on 
and  the  existing  stress  distribution  ^ , the  actual  mode  minimizes  the 
following  expression  when  geometry  changes  are  neglected  (Hill  [8]): 


1 • 2 

- 4 h(e*.)  dv  - 

2 iJ  Jc 


T.v. 

l l 


dS. 


(4) 


Note  that  the  virtual  mode  e* . in  tt0  should  be  normal  to  the  yield  surface 
at  the  existing  stress  point  in  the  stress  space  due  to  the  compatibility 
requirement  with  existing  stress  distributions.  For  statically  indeterminate 
problems,  however,  there  is  a coupling  between  stress  field  and  strain  rate 
field  and  we  have  to  solve  these  two  sets  of  variables  simultaneously. 


2.  Geometry  change 

When  the  effect  of  geometry  change  cannot  be  neglected  during  deforma- 
tion, it  is  necessary  to  reconsider  the  specification  of  the  loading  on  S^, 
and  the  stresses  since  the  changes  in  shape  and  area  of  surface  elements 
are  themselves  unknown. 

Let  be  the  position  vector  in  a Cartesian  reference  frame  at  time  t 
and  after  an  infinitesimal  time  fit,  x^  be  the  position.  Let  us  call  the 
configuration  at  time  t undeformed  configuration  and  the  one  at  time  t + fit 
deformed  configuration.  When  an  actual  force  dP^  acts  upon  the  area  element 
da  at  time  t + fit,  there  are  various  ways  of  reckoning  this  force. 

First,  the  actual  force  dP^  is  referred  to  the  deformed  configuration, 
or 

dP.  = n.o. .da,  (5a) 

x ) ij 
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where  is  the  unit  normal  vector  to  the  surface  element  of  area  da  in 
a deformed  configuration.  The  stress  tensor  defined  in  this  manner 
is  called  Cauchy  stress  tensor,  or  sometimes,  true  stress  tensor. 

Second,  the  actual  force  dP^  is  referred  to  the  undeformed  configura- 
tion, or 


dP.  = N.S. .dA, 
1 ) lj 


(5b) 


where  N.  is  the  unit  normal  vector  to  the  surface  element  of  area  dA  in 
J 

an  undeformed  configuration.  The  stress  tensor  defined  in  this  manner 
is  called  the  first  kind  of  Kirchhoff  stress  tensor,  or  sometimes,  nominal 
stress  tensor.  This  tensor  has  the  disadvantage  of  not  being  symmetric 
and  therefore  awkward  to  use  in  a constitutive  equation  with  a symmetric 
strain  tensor.  Nonetheless,  sometimes  this  stress  tensor  is  used  with 
nonsymmetric  velocity  gradients  [14] . 

Third,  to  obtain  a stress  tensor,  which  is  symmetric  and  referred  to 
the  undeformed  configuration,  we  proceed  as  follows.  Instead  of  the 
actual  force  dP.,  consider  a force  dP.  related  to  the  force  dP_  in  the 

l l 0 

same  way  that  a material  vector  dX^  is  related  by  the  deformation  to  the 


corresponding  vector  dx^.  That  is. 


3X. 

dP.  = 3-^  dP. . 

1 dX.  1 

J 


(5c) 


Refer  this  pseudo-force  dP^  to  the  undeformed  configuration  to  define  the 


second  kind  of  Kirchhoff  stress  tensor  t..: 

ij 


dP.  = N.T. .dA. 
i J iJ 


(5d) 


Using  the  expression  relating  the  area  change  of  the  same  element  during 
deformation  [15], 

8 


n.da  = — N.  -5—*-  dA,  (6) 

where  pQ  and  p are  densities  of  the  volume  element  before  and  after  the 
deformation,  the  relationship  between  different  stress  measures  is  obtained. 
From  Eqs . (6),  (5a),  and  (5b), 

P0  3xk 

dP.  = n.o.  .da  = o.  . — N,  3 — dA 
1 ] ij  ij  P k 3Xj 

(7a) 

= N.S. .dA 
J 13 


S . . " r.  O . « 

ij  p 3xr  lk 


and  from  Eqs.  (6),  (5a),  (5c),  and  (5d) , 


pQ  3X.  3)T 

Tij  p 3x.  °kr  3x 
k r 


All  these  different  stress  tensors  become  exactly  the  same  when  we 
bring  the  deformed  configurations  to  the  undeformed  configurations  and 
make  them  identical  in  the  limit.  Stress  rates,  however,  are  not  the  same. 
Let  6ui  be  the  increment  of  displacement  of  the  element;  then 


6S . . = 6a..  - a,  . 3—  6u.  to  the.  first  order, 
ij  tj  KJ  ox^  1 


neglecting  plastic  volume  change.  Or,  in  terms  of  rates. 


S. . = a. . - a..v.  . . 
ij  ij  kj  i,k 


Let  us  compare  the  magnitude  of  the  second  term  with  that  of  the  first  term 
in  the  right-hand  side  of  Eq.  (8).  Since 
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1 (3vi 

> 

H 

f3vi 

3vk] 

2 [axk 

axj  T I 

[axk 

' 9x.J 

i • 

i • 

2 eik  * 

2 “ik* 

where  is  the  rate  of  deformation  and  the  rate  of  rotation.  Then 


9v.  i . j 
°kj  alT  = 2 °kj£ik  * 2 °kjWik- 


Now,  if  workhardening  characteristics  are  given  by  the  relation 
a = H(e) , 

then 


H' 


da 

de  = FF 


For  an  order  of  magnitude  calculation  we  can  write  approximately 


Eik  = H' 


a,  .a., 
a e = JiLJJi 
°kjeik  H’ 


Then,  from  Eqs.  (9b)  and  (lOd),  Eq . (8)  becomes 

*ij  = ^ijf1  * 1v)  + I^ik- 


(10b) 
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From  liq . (11)  we  conclude  that  if  the  order  of  the  rate  of  rotation  w.  . 

ij 

is  the  same  as  or  less  than  the  order  of  strain  rate  e. and  if  the  work- 

ij 

hardening  rate  H'  is  greater  than  the  stress  level,  then  = cr  j . Other- 


wise, geometry  change  should  not  be  neglected. 

It  could  be  shown  [16]  that  when  geometry  change  is  taken  into  account, 
the  condition  for  continuing  equilibrium  requires  that 


3S.  . 

-s£-° 

1 


in  the  absence  of  body  force. 


Using  this  condition.  Hill  subsequently  derived  the  following  variational 
formulation  [17]: 


h(e?.)^  dv 
ij 


kj 


v*  , v*  .v  dv 

i,k  j,i 


T.vr  dS. 

1 1 


(12) 


Formulation  follows  essentially  the  same  line  of  formulation  except 
that  now  geometry  change  is  considered.  In  the  formulation  as  well  as 

in  7T-,,  virtual  mode  must  be  compatible  with  the  existing  stress  distribution 
and  the  boundary  condition  on  Sv.  As  has  been  discussed  earlier  for  statically 
indeterminate  problems  this  is  not  an  appropriate  formulation. 

Summarizing  the  development  so  far,  the  kinematic  mode  in  sheet-metal 
forming  of  a rigid-plastic  solid  is  not  uniquely  determined  by  considering 
the  first-order  expansion  of  the  potential  alone.  Consideration  up  to 
second-order  expansion  of  the  potential,  or  equivalent  consideration  of 
workhardening  rate  in  a physical  sense,  needs  stress  rate  terms  explicitly 
in  the  variational  formulation.  When  geometry  change  cannot  be  neglected, 
these  stress  rate-  are  related  to  stress  distribution,  which  is  not  known 
for  statically  indeterminate  problems.  The  approach  of  viewing  the  deformation 
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as  determining  the  incipient  flow  by  assuming  the  deformed  configuration 
coincident  with  the  undeformed  configuration  clearly  does  not  lead  to  a 
workable  variational  formulation  for  sheet-metal  forming  of  a rigid-plastic 
solid.  In  this  respect,  it  is  intended  to  develop  an  appropriate  variational 
formulation  in  the  next  section. 

t 


SECTION  III 


FINITE- ELEMENT  FORMULATION 


1 . Variational  formulation 

Let  be  the  position  vector  in  a Cartesian  frame  of  reference  at 

time  t,  the  moment  under  consideration.  Let  a. . be  the  true  stress  at 

ij 

time  t and  a. . + do. . the  true  stress  in  the  same  material  element  after 
1J  ij 

an  infinitesimal  time  dt,  both  tensors  being  associated  with  the  same 

Cartesian  axes.  Let  ds. . be  the  increment  in  nominal  stress  in  the  same 

1 3 

element  in  time  dt,  based  on  the  dimensions  at  time  t.  Let  du..  be  the 
increment  of  displacement  of  the  element,  then 


9(du.) 

ds . . = do . . - O . — » — — 
l]  13  j 8x. 


(13) 


Requiring  continuing  equilibrium  of  stresses,  the  virtual  work  principle 
gives 


3(du.) 


, 0 . . + do . . - 0.  . » — , 

l ij  iJ  kJ  9xk  J 


9(du.) 


i 9xi 


dV  = 


(E  ♦ dF.)6(du  ) ds, 
J q J J J 


(14) 


where  T.  = £.o.  . and  dT.  = £.ds..,  SL.  being  the  unit  normal  to  the  surface 
3 i i]  J i 13  i 

at  time  t.  The  variational  formulation  is  obtained  from  Eq.  (14)  as  follows: 

9(du. ) 3(du.) 


6<p  = 6 


o . . de . . dV  + 
ij  iJ 


1 2 

4 hdef . dV 
2 ij 


(J 


kj 


8xk  aXi 


(T.  + dT . ) du . ds  1 = 0, 
3 J J ' 
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where 


. f9(du  ) 9(du.) 

de  . . = -r-  — 5 + — -v-  — - 

il  2 I 3x.  9x. 

v J i 


and  h = j H' , with  H'  the  slope  of  the  stress  and  strain  curve.  The  first 
three  terms  of  the  functional  <J>  represent  the  energy  dissipated  during  the 
time  dt  up  to  the  second  order.  If  it  is  assumed  that  the  principal 
axes  of  true  strain-rate  keep  the  same  directions  in  the  element  and  the 
principal  components  of  strain-rate  maintain  the  constant  ratios  during 
the  time  dt,  the  dissipated  energy  can  be  expressed  directly  [18]  as 

T (a  dE  +4  hdE2) 

*•  p p 2 p 

per  unit  volume,  where  dEp  is  the  logarithmic  strain  components.  The  final 
form  of  the  functional  becomes 


_ - if  . 2 

<p  = odE  dV  4-  [ H'(dE)  dV  - (T.  + dT.)  du.  dS, 

. v 1 J c J J J 


where  dE  is  defined  by 


dE  = \ (dE  ) 2 


2.  Theory  of  the  finite-element  method 


An  important  step  in  finite-element  modeling  is  obtaining  approximate 
state  equations  in  a region.  The  weighted  residual  method  derives  the  state 
equations  directly  from  the  governing  differential  equations.  Let  us  write 
the  governing  differential  equation  as 


Lu  - f = 0, 


where  L is  the  differential  operator,  f is  the  known  function,  and  u is 
the  solution.  Withthe  trial  solution  u*,  Eq.  (17)  is  not  satisfied,  but 
there  remains  an  error  or  residual  R such  that 


R = Lu*  - f. 


(18) 


This  residual  is  multiplied  by  weight  function  w and  integrated  over  the 
domain  and  the  state  equations  are  derived  from  the  condition  that  this 
integral  vanishes  with  a given  choice  of  weight  function  w: 


wR  dv  = 0. 


(19) 


One  well-known  method  among  weighted  residual  methods  is  Galerkin's  approach. 

A more  frequently  used  approach  is  the  derivation  from  a variational 
principle  which  is  a dual  expression  of  the  governing  differential  equation. 
Assume  that  a functional  $,  which  is  equivalent  to  the  differential  equation, 
has  been  established.  Let  a continuum  be  divided  into  a finite  collection 
M of  subdomains  called  elements  interconnected  at  a finite  number  of  nodes 
N.  If  it  is  true  that  the  total  functional  is  equal  to  the  sum  of  the 
contributions  of  each  element  , then  we  may  write  as  follows: 


M 

* = l <f> 

m=l 


(m) 


(u). 


(20) 


In  each  element  let  us  approximate  the  solution  with  a linear  combination 
of  trial  functions  v^  such  that 


u 


* l 


a.  v. 
l l 


(21) 


holds,  where  ou  are  unknown  coefficients  to  be  determined  later.  By 
substituting  Eq.  (21)  into  Eq . (20),  we  have 
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(22) 


* = I <p(ra;)  (djV.) 
m=l 

M (m) 

= £ ^ J (aj)  since  v^'s  are  known 

m=l 

= $(a.)  • 

The  original  $ of  u is  now  discretized  with  a function  <p  of  parameters  a^, 
and  the  initial  variational  problem  reduces  to  determining  the  that 
minimizes  <p.  The  minimization  of  <f>  with  respect  to  or  may  be  written  as 

6<P  = f£"  6ai  = °»  (23) 

i 

where  6 denotes  the  first  variation.  Since  or  's  are  independent,  expression 
(23)  is  equivalent  to  a set  of  simultaneous  equations, 


(24) 


This  is,  in  fact,  the  classical  Ritz  technique.  It  is  the  choice  of  trial 
functions  that  makes  the  finite-element  method  different  from  the  Ritz 
method  and  renders  it  successful;  they  are  piecewise  polynomials.  Bsides, 
the  coefficients  or  , called  nodal  values  in  the  finite-element  literature, 
do  have  a definite  physical  meaning,  such  as  displacement  or  velocity. 

The  trial  function  v^  must  satisfy  certain  requirements  to  enable 
convergence  as  the  subdivision  into  ever  smaller  elements  is  attempted. 
First,  as  the  element  size  decreases,  the  functions  in  the  integral  must 
tend  to  be  single-valued  and  well  behaved  in  physical  problems.  This  is 
called  the  "completeness"  requirement  and  is  satisfied  if  the  trial  function 
is  of  class  CP  when  p is  the  highest  order  in  the  integrand  of  the 
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functional.  Second,  the  validity  of  the  summation  implied  in  Eq.  (20) 
must  be  preserved.  This  is  called  the  "compatibility"  requirement  and  is 
satisfied  if  v^  is  of  class  c^  * [19],  [20].  When  admissible  trial  functions 

are  used,  the  functional  converges  monotonically  with  an  increasing  number 

2 

of  elements  (or  decreasing  size)  at  a rate  proportional  to  h where  h is 
a characteristic  element  dimension. 


3.  Modeling  of  axisymmetric  problems 

The  general  outline  of  the  finite-element  modeling  stated  above  will 
be  expanded  in  detail  for  the  case  of  axisymmetric  thin  shells  subject  to 
axisymmetric  loading.  This  particular  problem  is  of  interest  since  some 
basic  sheet-metal  forming  processes  belong  to  this  category.  When  the  ratio 
of  thickness  to  the  radius  of  curvature  is  sufficiently  small,  bending 
moment  and  shearing  forces  may  be  neglected  without  serious  error  and  the 
membrane  theory  may  be  justified  [21].  Moreover,  the  state  of  stress  can 
be  treated  as  an  approximate  plane  so  long  as  — is  small  compared  with 
unity,  where  t is  the  local  thickness  and  s is  the  distance  in  any  direction 
parallel  to  the  surface.  We  now  may  rewrite  $ with  the  substitution  of 
t dA  = dv  to  Eq.  (16): 


$ = 5(dE)t  dA  + j H ' (dE) 2 (t  dA)  - 


(T  + dT)  dUj  dA 


for  the  unit  included  angle  of  the  element,  where  A is  the  area  of  the 
element  and  t is  the  sheet  thickness. 

From  the  symmetry  of  the  problem  it  is  easily  shown  that  the  circum- 
ferential direction  and  the  meridian  direction  are  the  principal  directions 
and  if  the  friction  between  the  shell  and  the  external  agent  is  negligible, 


the  thickness  direction  will  be  the  third  principal  direction.  Within  the 
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order  of  approximation  taken  in  the  formulation,  the  logarithmic  strain  incre- 
ment may  be  used  as  the  strain  increment  measure.  Then  the  definitions  of 
strain  increments  are 


dEl  dEr 


In  — 
s0 


In  — J 
r0 


if,  during  an  incremental  deformation,  an  element  of  undeformed  length  s0 
is  stretched  to  the  length  s and  the  point  currently  at  the  radial  distance 
rQ  moves  to  the  deformed  radial  location  r.  Subscripts  r,  0 refer  to  the 
meridian  and  the  circumferential  direction,  respectively. 

To  bring  the  model  closer  to  reality  in  the  present  investigation, 


normal  anisotropy  is  included  and  the  corresponding  stress-strain  increment 
relation  is  obtained,  using  Hill's  criterion  [13],  as 


(1  + R)a 


= dE 

(1  + R)a0  - R ar  (1  + R)5 


where  R is  the  planar  isotropy  parameter  which  is  the  ratio  of  width  strain 

to  the  thickness  strain  in  uniaxial  tension.  The  effective  stress  and  the 

f 

effective  strain  are  defined  as 


/2  2R  2 

a = /aQ  - o oa  + a , 

9 1 + R r 9 r 


dE  = -Ll±-  /dE*  + PVT  dE0dEr  + dE0' 
v 1 + 2R 


^Note  that  H'  = — must  be  consistent  with  these  definitions. 
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The  effective  strain,  dE,  may  be  written  in  matrix  form  as 


dE 


[dETD  dE]1/2, 


where 


D 


3(1  + R) 
= 2(1  + 2R) 


1 + R 
R 


R 

1 + R 


(29a) 


(29b) 


The  sheet  geometry  is  approximated  by  a series  of  conical  frustra,  as 
shown  in  Fig.  1.  Linear  trial  functions,  or  shape  functions,  as  they  are 
often  called  in  the  finite-element  literature,  are  enough  since  the  inte- 
grand in  the  functional  is  of  class  C* . The  unknown  coefficients,  or  nodal 
values,  are  taken  to  be  the  incremental  displacement  at  nodes.  Then  we  may 
write 


u(m)  = < dv^.dw^.dv^dw^) 


= ( du^,du2,du3>du^) 


(50) 


for  a representative  element  m,  where  dv^,  dw^  are  the  radial  and  the  axial 
components  of  incremental  displacement  of  the  i-th  node.  Then  the  incre- 
mental displacement  field  inside  the  element  may  be  written  as 


^dv 

1 + t' 

0 

1 - t' 

0 

1 

fdv] 

2 

2 

■ 

dw. 

1 

u = = 

[dwj 

0 

1 + t' 

0 

1 - t' 

dv2 

2 

2 

dw2- 

(31) 


= Nu 


(m) 


where  t'  is  the  local  coordinate  varying  from  the  value  of  -1  at  node  2 
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dW, 


Figure  1.  Approximation  of  the  Sheet  Geometry 
into  a Series  of  Conical  Frustra 

to  +1  at  node  1.  (See  Fig.  1.)  Due  to  this  incremental  displacement  field, 
an  element  of  length  sQ, 


so  ^(ro^i  “ 2)  * 

is  stretched  to  a new  length  s, 

s = >/^(r0)i  - (r0)2  + dvl  " dv2^2  + Uz'o^2  " + dw2  ‘ dwl^ 

= *^1  - r2)2  + {z2  - Zj)2, 


(32) 


where  (Tq)^  (Zq)^  are  radial  and  the  vertical  positions  of  the  i-th 
node  at  the  undeformed  configuration  and  (r)^,  (z)^  at  the  deformed 
configuration.  Since  the  element  is  straight,  any  point  of  t1  in  the  local 
coordinate  is  shown  to  have  a global  radial  position  rQ  determined  by 


r0  *•  2 ^r(Pl  + *•  2 ^r0^2' 


The  new  position  r of  the  same  particle  is  given  by 


(1  + t')  , (1  - t') 

r = r„  + -i — = dv,  + - — L dv 


(33) 


(34) 
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We  are  now  at  the  position  of  calculating  the  strain  increment  field. 
Recall  the  equation  (25)  and  substitute  Eqs.  (32),  (33),  and  (34)  into  it 
to  obtain 


dE  = 


1 , {(r0)l  - (r0}2  + dvl  ' dv2}  + {(Z0)1 

2 1" 2 


(ZQ)  2 + dw2  ' dwl^ 


In 


. (1  ♦ t')  J . (1  - t») 
ro  ♦ ~ 2 dvj  ♦ ' ' ' 2 dv2 


0 


(35) 


We  may  write  <J>^,  a contribution  from  the  m-th  element  to  the  total 
functional  $,  in  terms  of  nodal  values,  for  unit  angle  included: 

<f>W  =|  (odE  + i H ' (dE) 2 t dA  - | (T.  + dT.)v.  dA 

= | a(/j  t)  [dETD  dE]1/2  dA  + j j H'(|  t)  [dETD  dE]  dA  - j TTNu(m)  dA  , 


(36) 


where 


T = 


T1  + dTl 

T2  + dT2 

T3  + dT3 
t4  + dT4 


Minimization  gives  a set  of  simultaneous  equations: 


9(f) 


(m) 


9u 


(m) 


- T -1/2  9(d^T 
t)0[dElD  dE]  ' — Tjjjy  D dE  dA  ♦ 

9ul  ’ 


I <3 


t)H ' 


,(dE) 


D dE  dA 


9u 


/ \ U UL. 

(m)  - ~ 


- j NTT  dA  . 


(37) 


From  Eq . (35)  , 
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,(dE) 


9u 


(m) 


= a 


a^) 

3(dE2) 

ri  * r2 

1 * t' 

9uj 

9U1 

2 

s 

2r 

a(dEj) 

9(dE2) 

"(Z2  ' Zl} 

n 

9u2 

9u2 

2 

s 

u 

9(dEx) 

9(dE2) 

-(rl  ' r2} 

1 - f 

9u3 

9u3 

2 

s 

2r 

9(dE1) 

3(dE2) 

(Z2  ' ZP 

0 

9u4 

9u4 

2 

s 

(38) 


Therefore,  Eq.  (37)  becomes 

liii-  = f (|  t)0[|  dETD  dE]"1/2£D  dE  dA  + 
(mj  J •>  ~ 


(■|  t)H'(£D  dE  dA  - J NTT  dA  = 0. 


(39) 


These  equations,  being  valid  for  an  m-th  element,  are  now  to  be  combined 
undeT  the  condition  of  compatibility  that  the  first-order  derivative  of 
nodal  value  may  be  discontinuous  across  element  boundaries  but  the  nodal 
value  itself  must  be  continuous, 


3*  _ c 94> 

= 2 8u<»  = 

A/ 


(40) 


4.  Linearization 

Eqs.  (39)  and  (40)  are  nonlinear  equations  and  it  is  very  difficult 
to  solve  them  without  linearizing.  One  way  is  to  take  an  initial  guess  of 
the  solution  to  the  equation  as  u*  and  rewrite  Eq.  (39)  in  terms  of  the 
differences  between  this  initial  guess  and  the  correct  solution  Au,  where 
Ucorrect  = u*  + Au,  ar*d  expand  it.  Where  the  initial  guess  is  sufficiently 
close  to  the  correct  solution,  we  may  neglect  higher-order  terms  of  Au  and 
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thereby  linearize  successfully.  This  can  be  done  mathematically  in  a 
systematic  way  and  is  called  the  Newton-Raphson  method  [22].  Say  we  have 
a nonlinear  equation  tKu)  = 0,  then  we  may  expand  into  a series  with  respect 
to  the  correct  solution  u^  such  that 


<p(u) 


If  u and  Ug  are  sufficiently  close,  we  may  neglect  the  higher-order  terms 
and  write 


= ip*  + 


du 


Au  = 0. 


In  our  formulations  the  equations  to  be  minimized  are 


8u 


(m) 


0 , and , 


therefore,  the  expressions  corresponding  to  Eq.  (41)  are 


9u(m)8uWl 
i J 


(Au)  = 


9<pv 


9u. 


(41) 


(42) 


It  may  be  shown  that 


92C^m) 

3uW9uW 
i J 


= P 


(m) 


_1_ 

dE 


rv  2 bbT]  2 H'bbT\.  ,. 
O + H'  dE)  (K  - -=• + •=■ -K  dA 

^ dF  dE 


(43a) 


where 


b = QD  dE 
K = QDQT 
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and  that 


39(m)  = H(m)  _ p(m) 


1 


(43bj 


where 


,(m) 


= 2 f J_ 

3 J dE 


(a  + H'dE)bt  dA, 


pW  = nTx  dA. 


By  assembling  the  equations  obtained  for  an  element,  we  finally  have 


P*Au  = F - H* 


(44) 


We  evaluate  the  integrals  with  the  Gaussian  quadrature  formulation. 

We  have  yet  to  introduce  the  boundary  conditions  for  solving  a physical 
problem.  For  an  incremental  displacement  prescribed  boundary,  the  corres- 
ponding perturbations  should  vanish  and,  for  a traction  prescribed  boundary, 
the  prescribed  traction  value  will  enter  into  the  £ vector.  The  solution 
procedure  is  as  follows: 

(1)  Assume  an  initial  guess  u^,  and  compute  P,  H,  F corresponding 
to  this  guess. 

(2)  Solve  Eq.  (3.31)  and  obtain  Au. 

(3)  Obtain  a new  initial  guess  u2  = Uj  + Au. 

Repeat  this  process  until  convergence  is  achieved.  Convergence  is  checked 
by  the  fractional  norm.  A norm  is  defined  by  a square  root  value,  i.e., 

„ „ /2  2 

Hull  = ✓u1  + u2  + 

and 
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this  value  reaches  the  magnitude  smaller  than  a predetermined  value,  say, 
lO-6,  the  iteration  stops  and  the  solution  is  thus  obtained. 


SECTION  IV 


HYDROSTATIC  BULGING 

1.  Introduction 

The  ductility  of  sheet  metal  under  biaxial  stress  is  often  examined 
by  means  of  the  so-called  bulge  test.  A uniform  plane  sheet  is  placed  over 
a die  with  an  aperture  and  is  firmly  clamped  around  the  perimeter.  An  in- 
creasing hydrostatic  pressure  is  applied  to  one  side  of  the  sheet,  causing 
it  to  bulge  through  the  aperture.  From  the  measured  profile  and  thickness 
of  the  plastically  deformed  sheet  near  the  pole,  it  is  possible  to  calculate 
the  local  state  of  stress  in  terms  of  the  applied  pressure.  If,  in  addition, 
the  state  of  strain  is  measured  by  means  of  a grid,  the  stress-strain  char- 
acteristics of  the  metal  under  biaxial  tension  are  obtained.  The  advantage 
of  this  test  over  any  other  simple  one  is  that  a greater  range  of  pre- 
instability strain  can  be  obtained. 

Hydrostatic  bulging  is  not  only  important  as  a material  property  test, 
but  also  as  a forming  operation.  Thus,  a number  of  theoretical  investiga- 
tions, dealing  with  axisymmetric  hydrostatic  bulging  (Fig.  2)  has  appeared 
in  the  literature. 

The  classical  analysis  of  bulging  is  the  one  by  Hill  [23].  His  solu- 
tions are,  however,  special  ones.  Instead  of  analyzing  deformation  with 
a given  stress-strain  characteristic,  Hill  first  adopted  special  kinematic 
assumptions  and  from  them  deducted  the  necessary  stress-strain  character- 
istics which  satisfy  all  the  governing  equations  under  the  prescribed 
kinematic  mode.  The  kinematic  assumptions  are  first,  that  any  material 
element  describes  a circular  path  which  is,  moreover,  orthogonal  to  the 
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Figure  2 Schematic  view  of  hydrostatic  bulging. 

momentary  profile,  and  second,  that  circumferential  strain  is  numerically 
equal  to  the  tangential  strain.  The  required  stress-strain  characteristic 
is  found  to  be  an  exponential  type.  Hill's  other  solution  on  a linear 
workhardening  solid  uses  the  method  of  successive  approximation  by  adopting 
a yield  criterion  which  is  neither  von  Mises  nor  Tresca,  for  the  purpose 
of  mathematical  simplicity. 

Analyses  of  work  by  Woo  [24],  Yamada  [25],  and  Wang  [26]  are  based 
upon  the  realistic  choices  of  stress-strain  characteristics  and  the  yield 
criterion.  In  applying  the  deformation  theory  of  rigid  plasticity,  Wang 
experiences  a mathematical  difficulty  and  attributes  this  to  the  fact  that 
the  differential  equations  associated  with  the  deformation  theory  possesses 
a singularity  which  has  the  effect  of  restricting  the  range  of  calculation 
within  a certain  value  of  the  polar  strain.  Besides,  the  agreement  of 
deformation  theory  predictions  with  the  experiment  is  rather  poorer  than 
the  incremental  theory  prediction  [27] . 

In  applying  the  incremental  theory  of  rigid  plasticity,  researchers 
experience  a difficulty  in  satisfying  the  boundary  condition  at  the  fixed 
edge,  i.e..  Eg  = 0.  To  avoid  this  difficulty.  Woo  uses  the  deformation 
theory,  while  Yamada  reasons  that  introducing  an  elastic  strain  component 
into  the  formulation  will  resolve  this  "mathematical  difficulty"  (in 
Yamada' s terms)  and  turns  to  the  elasto-plastic  constitutive  law.  Another 
theoretical  work  of  interest  comes  from  Wang,  using  the  parametric  repre- 
sentation of  the  stresses. 
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The  only  published  solution  on  hydrostatic  bulging  using  the  finite- 
element  method  is  the  one  by  Iseki  et  al . [28],  with  the  incremental  theory 
of  elasto-plasticity . 

2.  Computational  procedures 

In  adopting  the  finite-element  model  to  hydrostatic  bulging,  it  is 
necessary  to  reconsider  the  external  work  increment  term,  since  the  pressure 
is  uniform  over  the  entire  surface  of  a closed  shell.  In  this  case  the 
increment  of  external  work  may  be  written  as  [29],  [30], 

Aw  = pVV,  (.45) 

where  VV  is  the  increase  of  the  volume  enclosed  by  the  deformed  sheet  and  p 

is  the  pressure  acting  on  the  deformed  configuration. 

As  an  initial  condition.  Hill's  special  solution  is  utilized.  In 

other  words,  the  initial  profile  of  the  bulge  is  assumed  to  be  a part  of 

1 a^ 

a sphere  whose  radius  is  given  by  r = ^ + h) , where  a is  the  radius 

of  the  original  blank  and  h is  the  polar  height  at  the  moment.  With  this 
geometry,  a pressure  p is  prescribed.  This  pressure  should  be  greater, 
at  least,  than  the  pressure  which  makes  the  sheet  having  initial  geometry 
everywhere  plastic.  The  initial  guess  on  the  incremental  displacement 
is  also  obtained  from  Hill's  special  solution  by  assuming  normal  trajectory 
of  the  element  particle  to  the  bulge  profile.  The  program  for  computing 
the  initial  guess  is  given  in  Appendix  A. 

When  a converged  solution  is  obtained  for  the  given  pressure,  a new 
bulge  profile  is  determined  from  the  initial  bulge  profile  and  incremental 
displacement  grid.  Then  the  pressure  is  assigned  a higher  value  and  the 
converged  solution  for  the  previous  step  is  used  as  the  initial  guess  for 
the  incremental  displacement  field  and  the  computation  continues  in  this 
way.  The  program  for  the  analysis  of  hydraulic  bulge  is  given  in  Appendix  B. 
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3.  Results  and  discussion 


To  examine  the  validity  of  the  present  FEM  for  hydrostatic  bulging, 
the  solution  is  compared  with  those  achieved  by  the  elasto-plastic  FEM 
and  the  experiment. 

The  following  conditions  were  employed  for  the  comparison  with  the 
elasto-plastic  FEM: 

- 0 2 9 2 

Workhardening  characteristics:  a = 105 (.0019  + e)  x 10  kg/m 

= 1 .036 ( .0019  + e)°’2  x 109N/m2 
-4 

Thickness:  3.0  x 10  m (=  0.3  mm) 

_2 

Radius  of  the  sheet:  2.4  x io  m (=  24  mm) 

Anisotropy  parameter:  1.0 

An  identical  problem  was  also  solved  by  Yamada  [25],  using  the  finite- 
difference  method  with  the  elastic-plastic  theory.  Fig.  3 shows  the  rela- 
tionship between  hydrostatic  pressure  and  the  polar  thickness  strain.  The 
solid  line  represents  the  elasto-plastic  FEM  (and  also  the  finite-difference 
method)  and  the  points  indicate  the  solution  given  by  the  rigid-plastic  FEM. 
The  deviation  of  the  first  point  by  the  rigid-plastic  FEM  is  thought  to 
reflect  the  approximation  involved  in  the  initial  condition  that  the  sheet 
is  everywhere  plastic  and  that  the  initial  geometry  is  a part  of  a sphere. 
The  solution  can  be  improved  numerically  by  taking  a smaller  value  of  h 
in  generating  the  initial  condition.  Nevertheless,  the  solutions  after 
this  first  step  are  in  extremely  good  agreement  with  the  elasto-plastic  FEM 
and  any  disturbance  in  the  initial  conditions  does  not  matter  after  an 
initial  deformation  of  a small  magnitude.  The  pressure  increment  is  raised 
by  twice  after  some  deformation  and  it  is  to  be  noted  that  the  solutions 
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with  the  larger  pressure  increment  size  are  still  accurate.  This  means 
that  the  method  is  computationally  economical  with  a reasonable  accuracy. 

After  the  last  point  in  the  diagram  the  solution  diverges  and  it  is  thought 
that  the  pressure  maximum  has  been  reached.  The  convergence  is  excellent; 
in  every  step,  five  to  seven  iterations  seem  to  be  sufficient.  Fig.  4(a), 

(b)  show  the  comparisons  of  strain  distributions.  The  circumferential  strain 
distributions  are  in  good  agreement.  The  tangential  strain  distribution  by 
the  rigid-plastic  FEM  deviates  somewhat  at  the  edge  from  that  by  the  elasto- 
plastic  FEM.  The  tangential  strain  is  more  sensitive  to  the  method  employed 
than  the  circumferential  strain,  but  this  deviation  of  tangential  strain  is 
not  serious  because  the  solution  closely  follows  that  by  the  finite- 
difference  method  and  we  may  conclude  that  the  strain  distribution  is 
accurately  predicted.  Fig.  5 shows  the  distributions  of  stresses  when  the 
polar  thickness  strain  is  (-0.4).  Fig.  6 shows  the  bulge  profile  at  some 
stages  of  deformation.  A number  of  material  elements  are  traced  during 
derormation  and  are  shown  on  each  bulge  profile. 

Next,  the  solution  is  compared  with  Mellor's  [31]  experiment  on  half- 
hard  aluminum. 

Workhardening  characteristics:  a = 15,460(1  + 0.76e)  psi 

= 1.066(1  + 0.76 e)  x 108N/m2 

Radius  of  the  sheet:  5.0  inches  = 1.27  m 

-4 

Thickness:  .035  inch  = 8.89  x 10  m ' 

Anisotropy  parameter:  1.0 

One  thing  to  be  mentioned  is  that  in  the  actual  experiment,  the  die  has 

3 

a round  profile  of  radius  g-  in.,  but  in  the  analysis  this  profile  has  been 

31 
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neglected.  Fig.  7 shows  that  the  agreement  of  the  relation  between  pressure 
and  polar  height  is  good.  The  agreement  in  the  bulge  profile  is  also  excel- 
lent, as  in  Fig.  8.  As  shown  in  Fig.  9(a),  (b) , the  theoretical  circum- 
ferential strain  still  closely  predicts  the  experimental  one,  but  there 
is  some  discrepancy  in  thickness  strain  distribution.  As  has  been  mentioned, 
the  actual  die  has  a round  profile  which  has  been  neglected  in  the  analysis, 
and  it  is  thought  that  the  thickness  strain  is  more  sensitive  to  the  profile 
than  is  the  circumferential  strain.  Initially,  there  is  virtually  no  dis- 
crepancy, but  increases  at  later  stages.  This  may  be  explained  by  the  fact 
that  initially  the  sheet  is  not  in  contact  with  the  profile,  but  as  deforma- 
tion continues,  more  of  the  sheet  is  biought  into  contact  with  the  profile 
and  makes  the  actual  situation  different  from  the  one  used  in  the  analysis. 

In  general,  the  theoretical  prediction  by  the  rigid-plastic  FEM  is 
in  good  agreement  with  both  the  experiments  and  the  analyses  by  the  elasto- 
plastic  FEM  and  the  finite-difference  method. 
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Figure  9.  (a)  Circumferential  Strain  Distribution;  (b)  Thickness  Strain  Distribution 


SECTION  V 


STRETCHING  OF  A SHEET  WITH  HEMISPHERICAL  PUNCH 

1 . Introduction 

Punch  stretching  is  commonly  used  to  assess  the  pressing  quality  of 
sheet  metals.  A circular  sheet  is  clamped  firmly  along  the  periphery  and 
is  stretched  by  a rigid  punch  of  hemispherical  shape.  The  depth  of  the 
deformed  sheet  when  it  fractures  is  usually  taken  as  a measure  of  ductility. 
See  the  schematic  diagram  in  Fig.  10. 

An  experimental  investigation  of  punch  stretching  as  a forming  problem 
dates  back  to  Loxely  and  Freeman  [32],  who  demonstrated  that  the  interfacial 
friction  between  the  punch  and  the  sheet  has  a significant  effect  on  the 
strain  distribution  in  the  sheet  and,  consequently,  on  the  location  of 
fracture  and  dome  height  when  the  sheet  fractures.  Keeler  and  Backofen  [33], 
in  characterizing  the  limit  stretching,  followed  the  strain  history  of  each 
element  with  the  progress  of  deformation  and  observed  the  occurrence  of 
discontinuity  in  tangential  strain  at  a certain  element,  which  was  subse- 
quently interpreted  as  the  onset  of  diffuse  necking  [34].  Based  upon  Hill's 
analysis  [35],  they  believed  that  localized  necking  is  not  possible  in  punch 
stretching,  but  that  only  diffuse  necking  takes  place,  increasing  the  over- 
all nonuniformity  of  straining. 

The  observation  of  localized  necking  in  situations  where  Hill's 
analysis  denies  one  has  been  well  established  in  the  case  of  in-plane 
stretching  and  has  prompted  the  development  of  Marciniak  and  Kuczynski's 
theory  [36],  [37],  In  punch  stretching.  Gosh  and  Hecker  [38]  observed 
localized  necking  and  reported  that  local  necking  sets  in  even  though 


38 


I 


the  plane-strain  condition,  which  is  thought  to  be  responsible  for  local 
necking  in  in-plane  stretching,  is  not  achieved.  This  is  attributed  to  the 
fact  that  in  punch  stretching  an  increment  in  tangential  strain  is  geometri- 
cally tied  to  an  increment  in  circumferential  strain  and,  therefore,  the 
approach  to  plane-strain  condition  becomes  slower.  Another  experimental 
investigation  of  punch  stretching  is  the  one  by  Alexander  and  Kaftanoglu 
[39].  They  observed  that  the  deformation  is  limited  by  the  "strain 
propagation  instability"  or,  local  necking  in  common  terminology,  and  not 
by  "maximum  load  instability"  or,  diffuse  necking. 

From  the  viewpoint  of  the  deformation  analysis,  punch  stretching  is  a 
complicated  problem  because  a moving  boundary  separates  the  region  in  contact 
with  the  punch  head  from  the  unsupported  one.  The  friction  over  the  punch 
head  gives  rise  to  additional  complications.  One  special  solution  is  by 
Chakrabarty  [40].  Following  the  line  of  Hill's  special  solution  on  hydro- 
static bulging  he  obtained  an  analytical  solution  for  a special  material 
having  exponential  type  stress-strain  characteristics.  For  more  general 
materials  the  only  solutions  available  are  the  numerical  ones.  Numerical 
solutions  of  importance  are  those  by  Woo  [41]  and  by  Wang  [42],  [43]. 

Woo's  and  Wang's  solutions  were  obtained  by  the  finite-difference 
method.  The  only  solution  by  the  finite-element  method  on  punch  stretching 
is  one  by  Wifi  [44].  His  elasto-plastic,  finite-element  model  does  not 
neglect  the  bending  moment  nor  the  effect  of  shear  stress  and  uses  two- 
dimensional  triangular  elements  to  take  the  thickness  variation  into 
account.  Friction,  which  is  of  primary  significance  compared  with  the 
secondary  effect  of  bending  and  thickness,  is  assumed  to  be  perfect,  meaning 
that  once  the  element  touches  the  punch  head,  it  does  not  slide  over  the 
punch  but  sticks  to  it. 
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2.  Computational  procedures 

In  applying  the  finite-element  method  to  punch  stretching,  a thought 


should  be  given  to  the  implementation  of  boundary  conditions.  The  boundary 
conditions  in  punch  stretching  are  stated  not  only  by  prescribing  tractions 
and  incremental  displacements  but  sometimes  by  their  ratios.  In  this 
report,  the  problem  is  similar  to  the  ball  indentation  problem  (Lee  et  al . 
[45]). 

The  radial  and  vertical  positions  of  the  material  elements  in  the 
contact  region  are  not  independent  but  they  are  related  to  each  other 
through  a mathematical  expression  for  the  geometrical  requirement  that  they 
must  be  actually  on  the  surface  of  the  punch  head.  The  expression  is 

(rQ  + v)2  + (c  + zQ  + w)2  = r2,  (46) 

where  r^,  z^  are  radial  and  vertical  positions  of  the  element  at  the  present 
undeformed  configuration;  v,  w are  the  increments  of  horizontal  and  vertical 
displacements,  and  c is  a parameter  related  to  the  punch  height  h by  the 
expression 

c = r - h. 

P 

See  Fig.  11.  Recall  that  the  finite-element  formulation  in  Chapter  IV  has 
already  been  linearized  and  what  it  really  solves  for  are  the  perturbation 
terms.  Therefore,  we  also  linearize  the  boundary  condition  (46)  to  obtain 

2(rQ  + v*)Av  + 2(c  + z^  + w*)Aw  = r2  - (rQ  + v*)2  - (c  + z^  + w*)2,  (47) 

where  starred  (*)  quantities  are  initial  guesses,  and  Av,  Aw  are  perturba- 
tions. By  rearranging  (47),  we  have 
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Figure  11.  Geometrical  Requirement  for  the 
Node  on  the  Contact  Region 
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(48a) 


Av  = — Aw  + 6, 
a ’ 


where 


a = 


(r0  + v*>  _ 1 
(c  + Zy  + w*)  tan  0 


and 


2 2 2 
rp  - (c  + zQ  + w*)  - (rQ  + v*) 

6 = 2(r0  + v*) 


(4  8b) 


(4  Sc) 


When  the  finite-element  model  is  implemented,  all  the  tractions  are  trans- 
formed into  generalized  nodal  forces.  Therefore,  it  is  convenient  to  write 
the  boundary  condition  in  terms  of  the  generalized  nodal  forces  tt^  and  tt^  , 
the  horizontal  and  vertical  components,  respectively.  See  Fig.  11. 

Now 

ir.  . = N cos  0 - S sin  0 

(r)  (49) 

it,  . = N sin  0 + S cos  0 
(z) 

where  N and  S are  generalized  forces  normal  and  tangential  to  the  punch 
head.  We  eliminate  N through  the  relation 

2 2 
cos  0 + sin  0=1 

and  obtain 

tt  cos  0 = tt  sin  0 + k,  (50) 

z r 

where  k is  the  frictional  force  at  nodes.  However,  from  geometry  we  know 
that  the  following  holds: 
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cos  0 


r0  + V 


z„  + w"  + c 


sin  0 = 


So,  (50)  may  be  written  as 


tt  ( , kr 

* ♦ -ill  = E 

(z)  a (rQ  + v*) 


If  the  die  has  a round  profile  of  the  radius  rD>  then  the  requirement 
for  a material  element  to  lie  geometrically  on  the  profile  is  similar  to 
the  requirement  to  be  satisfied  on  the  punch  head.  Therefore,  we  have 
(similar  to  Eq.  (46)), 


U - r0  - v)2  ♦ (rD  - zQ  - w)2  = r2. 


where  a is  the  radius  of  the  sheet.  Linearization  of  Eq.  (S3)  gives 
2(a  - rQ  - v*)Av  + 2(rQ  - zQ  - w*)Aw  = -r2  + (a  - rQ  - v*)2  + (rD  - - w*) 


or,  rewriting. 


a Aw  „ 
Av  = — + ft, 
Y 


where 


(a  - rQ  - v*) 
(rD  - z0  ' W*} 


(a  - r0  - v*)2  * (rp  - zQ  - w*)2  - r2 
2(a  - rQ  - v*) 


The  tractional  boundary  condition  over  the  die  profile  can  be  written 
similarly  as 


For  the  portion  of  the  sheet  which  is  not  in  contact  with  the  punch 
head  nor  with  the  die  profile,  the  displacement  increment  in  the  radial 
direction  and  the  displacement  increment  in  the  axial  direction  are  not 
bound  to  each  other,  as  is  the  case  for  the  contact  region,  but  remains 
as  independent  variables.  Tractions  are,  however,  given  the  value  of  zero. 

With  the  advancement  of  the  punch  head,  the  portion  of  the  sheet  in 
contact  with  the  punch  or  die  profile  increases  and,  consequently,  the 
boundary  separating  this  "contact  region"  from  the  "unsupported  region" 
changes.  The  presence  of  this  moving  boundary  is  always  a source  of  compli- 
cations in  the  numerical  analysis  of  punch  stretching  because  it  requires  a 
basically  trial-and-error  approach.  The  treatment  of  the  moving  boundary 
used  in  the  present  analysis  for  punch  stretching  without  round  die  corners 
is  explained  in  detail  as  follows: 

First,  assume  the  position  of  the  boundary  in  future  configurations. 

In  the  FEM  this  means  assuming  which  nodes  will  be  in  contact  with 
the  punch  head  in  the  future  configuration.  Then,  obtain  a converged 
solution  based  upon  this  assumption  and  check  to  see  if  it  is  true. 
Since  the  position  of  the  boundary  is  already  known  in  the  current 
configuration,  in  practice  we  assume  and  check  how  much  this  boundary 
advances . 

(1)  Check  whether  the  boundary  is  assumed  to  advance  too  fast. 

Compute  the  normal  component  of  the  generalized  nodal  force 
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to  the  punch  head  for  the  nodes  in  contact  with  the  punch  head. 

If  every  generalized  normal  force  is  directed  outward  from  the 
punch  head,  then  all  the  nodes  which  are  assumed  to  be  in  contact 
with  the  punch  head  actually  do  so.  On  the  other  hand,  the  gener- 
alized normal  force  in  the  direction  toward  the  punch  head  for 
any  node  means  that  external  force  other  than  the  one  exerted 
by  the  punch  is  necessary  for  this  particular  node  to  conform 
with  the  punch  geometry.  Since  there  is  physically  no  source 
of  applied  force  other  than  the  punch,  the  assumption  that  this 
particular  node  is  in  contact  with  the  punch  head  is  not  correct 
and  the  position  of  the  boundary  should  be  re-assumed  to  exclude 
this  node  from  the  contact  region. 

(2)  Check  whether  the  boundary  is  assumed  to  advance  too  slowly. 

Compute  the  distance  between  the  nodes  in  the  unsupported  region 
and  the  center  of  the  hemisphere  of  the  punch  head.  If  this 
distance  is  shorter  than  the  radius  of  the  punch  head  for  any 
node,  it  means  that  this  particular  node  is  inside  the  punch 
head.  Since  this  is  physically  impossible,  the  assumption  that 
this  particular  node  is  not  in  contact  with  the  punch  head  is 
not  correct  and  the  position  of  the  boundary  should  be  re- 
assumed to  include  this  particular  node. 

Although  this  basically  trial-and-error  approach  seems  to  be  very 
time  consuming,  in  actual  computation  we  can  predict  the  movement  of  the 
boundary  fairly  accurately  based  upon  the  distance  between  the  free  nodes 
and  the  punch  surface.  Furthermore,  since  we  already  know  the  position  of 
the  boundary  in  the  current  configuration,  it  is  enough  to  check  the 
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boundary  assumption  for  only  a few  nodes  neighboring  the  previous  position 
of  the  boundary,  not  for  whole  nodes. 

The  procedure  described  above  for  the  contact  region  on  a punch  head 
is  also  applicable  for  the  contact  region  at  the  die  corner.  Handling  two 
moving  boundaries  simultaneously  really  does  not  invoke  any  new  theoretical 
difficulties  but  only  takes  more  computation  time  and  may  be  impractical 
for  inefficient  numerical  methods. 

In  order  to  implement  Coulomb  friction  between  sheet  and  punch  or  die, 

we  first  prescribe  a tangential  friction  force  S and  obtain  a converged 

solution  and  then  compute  generalized  nodal  forces.  From  Eqs . (49)  we  then 

are  able  to  compute  the  normal  component  N and  the  friction  coefficient 
S 

p = ^ corresponding  to  the  initially  prescribed  value  of  S.  If  the  com- 
puted friction  coefficient  is  not  what  is  intended,  then  we  modify  the 
S value  and  repeat  the  process.  It  should  be  noted  here  that  the  correction 
of  frictional  force  S needs  the  necessary  modification  only  in  the  F matrix 
(Eq.  (44)),  while  the  stiffness  matrix  P,  which  is  the  most  time-consuming 
part,  remains  the  same. 

The  deformation  step  is  controlled  by  the  punch  head  increment,  which 
is  designed  in  the  present  codes  to  yield  the  maximum  increment  of  effec- 
tive strain  roughly  equal  to  a preset  value.  In  the  present  work  the 
optimum  size  is  shown  to  be  a 0.04  increment  of  effective  strain.  The 
solution  generally  converged  after  10  15  iterations  for  a single  step 

within  the  fractional  norm  of  10  ^ . The  actual  program  is  shown  in  Appendix  C. 

3.  Results  and  discussion 

The  present  rigid-plastic  FEM  is  compared  with  the  finite-difference 
methods  by  Wang  [43]  and  Woo  [41],  and  also  with  the  experiment  by 
Kaftanoglu  and  Alexander  [39]. 
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(1)  Comparison  with  the  finite-difference  solution  by  Wang 
The  parameters  used  in  Wang's  example  are  as  follows: 

Material : copper 

0 326  2 

Stress-strain  characteristics:  a = 30.5eu'  ton/in. 

. 0.326  . _8  . 2 

= 4.6361e  x io  N/m 

Anisotropy:  R = 1.0 
Friction:  y = 0.04 

Thickness:  t = 0.035  in.  = 8.89  x 10~2  m 

_2 

Punch  radius:  r = 1.0  in.  = 2.54  x 10  m 
P 

Radius  of  sheet:  a = 1.15  in.  = 2.921  x 10  2 m 
Initial  radius  is  sometimes  denoted  by  r^. 

The  two  methods  are  in  excellent  agreement  in  predicting  the  punch 
head  for  a given  punch  travel.  See  Fig.  12;  the  solid  line  represents 
Wang's  solution  and  the  points  represent  the  rigid-plastic  FEM.  Fig.  13 
shows  the  thickness  strain  distribution.  Again,  a good  agreement  between 
the  two  solutions  is  apparent. 

The  second  example  has  the  following  parameters: 

0 2 

Stress-strain  characteristics:  a = ke 

Anisotropy:  R = 1.0 

Friction:  y = 0.2 

Punch  radius:  r =1.0 
P 

Radius  of  sheet:  r^  = 1.0 

In  Wang's  work  all  the  results  are  reported  in  the  dimensionless  number. 

Figs.  14  and  15  show  the  circumferential  strain  distribution  and  thickness 
strain  distribution,  respectively.  The  solid  line  represents  Wang  and  points 
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Punch  load 


Punch  travel 

Figure  12.  Punch  Head  vs.  Punch  Travel 
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Thickness  strain 


represent  the  rigid-plastic  FEM.  Excellent  agreement  of  the  two  solutions 
is  demonstrated. 

The  step  size  has  an  important  effect  upon  the  accuracy  and  efficiency 
of  the  solution.  The  smaller  the  step  size,  the  better  the  accuracy,  al- 
though more  computation  time  is  required.  Fig.  16  demonstrates  that  there 
is  a limit  to  increasing  efficiency  while  maintaining  accuracy.  For  example, 
solutions  with  a step  size  of  0.08  in  the  effective  strain  increment  deviates 
considerably  from  the  solutions  obtained  with  step  sizes  of  .02  or  .04.  In 
the  remainder  of  the  work  the  step  size  of  .04  is  most  often  used. 

Compared  with  this  significant  effect  of  step  size,  the  mesh  size  does 
not  exert  a great  influence  upon  the  solution,  as  is  demonstrated  in  Fig.  17. 
The  solution  with  a coarse  mesh  (10  elements)  is  essentially  the  same  as 
the  one  with  a finer  mesh  (40  elements),  even  though  the  latter  will  be 
helpful  in  pinpointing  the  exact  location  of  peak  strain. 

In  the  examples  above,  there  is  only  one  moving  boundary,  that  between 
punch  and  sheet,  since  the  presence  of  the  round  die  profile  is  neglected. 

In  practice,  the  die  always  has  a round  profile  and  as  the  radius  of  the 
profile  gets  larger,  it  becomes  necessary  to  include  the  die  profile  in 
the  analysis.  In  this  case  there  are  two  moving  boundaries,  the  second 
being  the  one  between  sheet  and  die.  The  only  work  reported  which  includes 
the  die  profile  into  the  analysis  is  the  one  by  Woo. 

(2)  Comparison  with  the  finite-difference  solution  by  Woo 

The  parameters  in  Woo's  example  are: 

0 504  2 

Stress-strain  characteristics:  a = 5.4  + 27. 8e  ‘ ton/in. 

0 504  9 2 

for  e < 0.36:  = (0.08208  + 0.422569e  ) x 10  N/m 

_ . ...  0.375  _ 2 

= 5.4  + 24. 4e  ton/in. 

0 375  9 2 

for  e > 0.36:  = (0.08208  + 0.37089  ) x 10  N/m 
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Material : copper 

_2 

Punch  radius:  1 in.  = 2.54  x 10  m 

_ X 

Die  profile  radius:  0.3  in.  = 7.62  x 10  m 
Radius  of  sheet:  1.3  in.  = 3.302  x 10  2 m 
Coefficient  of  friction:  0.04 
Thickness  of  sheet:  0.035  in.  = 8.89  x lo”4  m 

Figs.  18  and  19  are  the  thickness  strain  distribution  and  the  circumferential 
strain  distribution.  Solutions  by  Woo  are  represented  by  solid  lines  and 
the  solutions  by  the  rigid-plastic  FEM  are  represented  by  points.  Agree- 
ment between  the  two  solutions  is  excellent  for  most  of  the  deformation. 
However,  at  later  stages  of  deformation,  a discrepancy  is  observed  around 
the  edges.  Re-examining  Woo's  computational  procedure  reveals  that  in 
order  to  avoid  the  difficulty  of  satisfying  boundary  conditions  exactly 
along  the  fixed  edge  =0),  he  allowed  a small  increment  of  circumfer- 
ential strain  along  the  edge  at  each  stage.  In  the  present  rigid-plastic 
FEM  such  difficulty  does  not  exist,  so  there  is  no  need  to  relax  the 
boundary  condition.  The  discrepancy  observed  at  later  stages  of  deformation 
may  be  attributed  to  this  difference. 

With  regard  to  the  instability.  Woo  stated  that  it  occurs  when  the 
resultant  tangential  stress  determined  from  the  strain  hardening  character- 
istics cannot  obtain  the  value  required  for  the  equilibrium  and  at  that 
instant  he  stopped  the  computation.  In  the  present  rigid-plastic  analysis 
such  an  instability  is  not  observed  at  the  point  reported  by  Woo,  and  the 
computation  continues. 


The  parameters  of  Kaftanoglu  and  Alexander's  experiment  on  soft  copper 

are: 
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0 3789 

Stress-strain  characteristics:  a = 68,394(0.0122  + e)  ' psi 

a 7iri  rn  run  -.0.3789  , 8 . 2 

= 4.7156  (0.0122  + e)  x 10  N/m 

Thickness:  0.048  in.  = 1.219  x 10-3  m 

Friction  condition:  PTFE  film  lubricant 

Radius  of  the  sheet:  0.717  in.  = 1.821  x 10~3  m 

_3 

Punch  radius:  0.65  in.  = 1.651  x 10  m 

Kaftanoglu  reports  that  the  friction  condition  changes  with  deformation 
and  measures  three  different  friction  coefficients:  p = 0.2  at  stage  1, 
p = 0.135  at  stage  2,  and  p = 0.07  at  stage  3.  To  include  the  changing 
friction  coefficient  into  the  analysis,  we  need  more  information  on  the 
friction  history,  which  is  difficult  to  obtain  experimentally.  Therefore, 
as  a representative  value,  we  use  the  mean  of  three  values  of  the  friction 
coefficient,  p = 0.135,  for  our  computation.  Figs.  20  and  21  show  the 
distribution  of  the  circumferential  strain  and  the  thickness  strain.  The 
agreement  between  the  experimental  data  and  the  numerical  solution  is  a 
reasonable  one  considering  the  fact  that  the  exact  friction  condition  is 
not  known. 

(4)  Influence  of  formulation  of  constitutive  relation 

Various  formulations  have  been  given  for  plastic  stress-strain  rela- 
tionships of  workhardening  materials.  Among  them,  the  parabolic  hardening 
law  has  been  used  extensively  for  sheet  metals  because  of  the  ease  with 
which  it  characterizes  workhardening  properties  of  materials.  However,  it 
was  suggested  recently  [46]  that  the  Voce  equation  [47]  is  a better  repre- 
sentation of  materials  behavior  when  solving  plasticity  problems  involving 
workhardening  rate.  The  forming  limit  curves  were  compared  using  the 
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parabolic  hardening  law  and  the  Voce  equation  [48] , and  the  result  indicates 
an  importance  of  the  choice  of  workhardening  representation  of  materials. 

Because  the  term  containing  the  rate  of  workhardening  appears  in  the 
finite-element  formulation  of  sheet-metal  forming,  it  is  of  importance  to 
examine  the  influence  of  workhardening  representation  on  the  mechanics 
computed  by  the  finite-element  method.  The  material  is  aluminum  alloy 
2036-T4.  The  parameters  are  as  follows: 

0 222 

Stress-strain  characteristics:  a = 86, 000(e)  ' psi  for  the  parabolic 

hardening  law 

a = 65,000(1  - (1  - 0.508) exp(-8.51e)}  psi 

for  the  Voce  equation 

Fig.  22  shows  the  two  stress-strain  curves  together  with  tension  test  data 
from  the  specimens  cut  in  the  three  directions  (0°,  45°,  90°).^ 


(1)  The  stress  and  strain  values  in  tension  tests  in  the  three  directions 
were  converted  to  values  of  the  effective  stress  and  effective  strain 
according  to 

(a)  Tension  in  the  0°  (rolling)  direction: 


5 ■ /?/ . TV.  •..XSL 


r0  + r90  + r0r90 


V 


- = /2  /r0  * r90  * V90 
3 r90  + r0r90 


C0* 


(b)  Tension  in  the  45°  direction: 


- = /I  AO  * r90)(1  * r45^~ 

2 2^r0  + r90  + r0r90}  a45’ 


/2  /2(-r0  + r90  + r0r90^ 
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(c)  Tension  in  the  90°  direction: 
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(Footnote  continued  on  next  page) 
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Figure  22.  Stress-strain  Curve  for  A1 . 2036-T4 
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r-value:  r„.  = 0.66.  r._  = 0.69,  r..  = 0.70,  and  r = 0.685 
0 45  90  a 

Radius  of  die  opening:  (0.80  in.) 

Blank  thickness: 

Radius  of  punch  head:  0.75  in.  and  0.45  in. 

Coefficient  of  friction:  0 and  0.2 

Punch  stretching  was  performed  on  a horizontal  hydraulic  press.  Tests  were 
interrupted  for  strain  measurements  (thickness  and  circumferential  strains) 
from  the  grids  photoprinted  on  the  specimen.  Load-displacement  relationships 
were  also  recorded.  First,  the  experimental  strain  distributions  were  com- 
pared with  computed  results,  using  the  parabolic  hardening  law  in  Fig.  23. 

In  the  experiment  Johnson's  wax  was  used  as  the  lubricant  and  was  applied 
at  each  stage.  In  comparison,  two  discrepancies  are  apparent:  (i)  the  co- 
efficient of  friction  does  not  stay  constant;  particularly,  at  the  last  stage, 
the  experimental  strain  distributions  indicate  that  the  coefficient  of  fric- 
tion is  less  than  0.2,  which,  however,  gives  good  agreement  for  other  stages, 
and  (ii)  the  measured  thickness  and  circumferential  strains  for  a given 
punch  depth  do  not  follow  the  corresponding  theoretical  curves.  This  is 
attributed  to  the  fact  that  the  accurate  strain  measurements  is  extremely 
difficult  for  critical  comparison  between  theory  and  experiment.  The  load 
values  summarized  in  Table  1 show  an  excellent  agreement  between  the  two. 


where  r„,  r,r,  r„„  are  the  r-values  obtained  from  the  tension  of  speci- 
0 45  90  r 

mens  cut  in  the  0°,  45°,  and  90°  directions,  respectively. 


(2)  The  effective  stress  and  effective  strain  defined  in  the  formulation  of 
this  report  differ  from  the  definition  above  by  a factor  such  as 
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a = V -zr  -Z , 
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where  r is  the  average  r-value  defined  by  rg 
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Table  1 

PUNCH  LOAD  AND  DISPLACEMENT  RELATIONS 


Punch  head  radius  = 19.05  mm  (0.75  in.) 


Theoretical 

Displacement  Punch 

mm  (in.)  N 

load 

(lb) 

Punch 

(lb) 

Experimental 

load  Displacement 

N 

4, 

.06 

(0. 

.160) 

6,330 

(1,423) 

(970) 

4,315 

2. 

79 

(0. 

.110) 

6, 

.10 

(0. 

.240) 

10,889 

(2,448) 

(1,730) 

7,695 

4 . 

.83 

(0. 

.190) 

7. 

.54 

(0. 

.297) 

14,483 

(3,256) 

(2,990) 

13,300 

7. 

11 

(0. 

.280) 

9, 

.80 

(0. 

.386) 

22,059 

(4,959) 

(4,940) 

21,974 

10. 

08 

(0. 

.397) 

12. 

.45 

(0. 

.490) 

30,301 

(6,812) 

(6,580) 

29,269 

12. 

45 

(0. 

,490) 

Punch  head  radius 

= 11.43  mm  (0.45 

in.) 

4 

.06 

(0. 

,160) 

5,124 

(1,152) 

(920) 

4,092 

2. 

72 

(0. 

,107) 

6. 

.10 

(0. 

,240) 

8,131 

(1,828) 

(2,000) 

8,896 

6. 

30 

(0. 

,24  8) 

8. 

.53 

(0. 

,336) 

12,237 

(2,751) 

(2,770) 

12,322 

8. 

18 

(0. 

322) 

9. 

.58 

(0. 

,377) 

13,963 

(3,139) 

(3,130) 

13,923 

9. 

68 

(0. 

381) 

For  a smaller  punch  size,  the  strain  distributions  are  compared  in  Fig.  24. 
The  same  observations  as  those  in  Fig.  23  apply.  Again,  the  punch  load  is 
in  good  agreement . 

The  influence  of  workhardening  representations  on  the  detailed  mechanics 
is  examined  in  Figs.  25,  26,  27,  and  28.  Referring  to  Fig.  25,  the  general 
trend  of  strain  distributions  is  not  altered  by  the  workhardening  representa- 
tion. However,  the  magnitude  of  strains,  particularly,  peak  strains,  differ. 
With  the  Voce  equation,  the  peak  strains  are  larger  than  those  computed  by 
the  parabolic  workhardening  law.  This  difference  becomes  larger  as  the 
punch  penetrates. 

It  appears  that  the  difference  of  the  two  is  more  significant  for 
higher  friction  in  the  larger  punch  size.  However,  in  the  smaller  punch 
size,  the  difference  of  the  two  strain  distributions  is  about  the  same  for 
the  two  coefficients  of  friction,  0 and  0.2,  as  shown  in  Fig.  26. 

It  is  rather  surprising  to  find  in  Figs.  27  and  28  that  the  punch  load 
for  the  same  punch  displacement  is  higher  with  the  parabolic  workhardening 
law  than  that  with  the  Voce  equation.  The  difference  becomes  significant 
for  large  punch  penetration.  From  these  results,  it  is  concluded  that  the 
representation  of  the  workhardening  characteristics  of  the  material  does 
have  an  influence  on  the  computed  strain  distributions  and  load-displacement 
relationships.  The  difference  becomes  critical  for  large  punch  displacement 
in  predicting  both  peak  strains  and  the  punch  load.  In  order  to  determine 
which  representation  is  preferable,  however,  more  experiments  with  improved 
accuracy  and  control  are  needed. 
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Figure  27.  Comparison  of  Theoretical  Load  Displacement 
Curves  Using  (1)  the  Parabolic  Workhardening 
Law  and  (2)  the  Voce  Equation  for  p = 0 


SECTION  VI 


DEEP  DRAWING  OF  A SHEET  WITH  HEMISPHERICAL  PUNCH 

1.  Introduction 

In  a deep  drawing  test  a circular  sheet  of  metal  is  placed  between 
the  blank  holder  and  the  die  and  then  fully  drawn  into  the  shape  of  a cup. 

The  formability  is  then  measured  by  the  maximum  size  of  the  blank  which  can 
be  drawn  without  a failure,  or,  more  often,  by  its  ratio  to  the  punch  diameter. 
This  ratio  is  called  the  limiting  drawing  ratio  and  this  particular  kind  of 
test  is  called  the  Swift  test. 

Deep  drawing  is  not  only  a useful  method  of  material  testing,  but  also 
one  of  the  basic  operations  in  sheet-metal  stamping.  In  practice,  various 
shapes  are  possible  for  the  bottom  of  the  punch;  however,  most  past  investi- 
gations are  on  deep  drawing  with  a flat-bottomed  punch  [49] -[56]. 

Among  the  earlier  works  on  deep  drawing  are  those  by  Hill  [13]  and  by 
Chung  and  Swift  [52]  using  the  incremental  theory  of  plasticity.  More 
refined  analyses  are  the  finite-difference  solutions  by  Chiang  and  Kobayashi 
[57],  b-  Wang  and  Budiansky  [51],  and  by  Chakrabarty  and  Mellor  [49].  Even 
though  such  a refinement  improves  the  understanding  of  the  deep  drawing 
process,  their  works  are  not  complete  because  they  treat  the  deep  drawing 
problem  as  an  in-plane  pure  radial  drawing  and  are  concerned  mostly  with 
the  deformation  mechanics  on  the  flange.  However,  it  has  been  observed 
experimentally  (Chung  and  Swift  [52])  that  the  die  profile  and  the  punch 
profile  significantly  affect  the  punch  load  and  the  strain  distributions 
and  therefore  a further  refinement  is  necessary  by  considering  these 
parameters  in  the  analysis.  Woo  [53]  performs  such  an  analysis  and  then 
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is  able  to  show  that  the  solution  obtained  by  extrapolating  the  strain 
distribution  over  the  flange  to  the  die  throat  predicts  more  straining 
than  the  one  obtained  by  taking  the  profiles  into  consideration. 
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Contrary  to  these  numerous  investigations  on  deep  drawing  with  a flat- 
bottomed  punch,  very  few  works  are  reported  on  the  deep  drawing  of  a sheet 
with  a hemispherical  head  punch  (Fig.  29).  Woo  [58]  analyzes  this  problem 
by  breaking  down  the  deep  drawing  process  into  two  component  processes  of 
the  pure  radial  drawing  over  the  flange  and  the  punch  stretching  over  the 
hemispherical  punch  head.  He  first  obtains  solutions  for  pure  radial  drawing 
in  the  flange  and  then  uses  this  solution  at  a point  initially  situated  near 
the  die  lip  as  the  boundary  condition  for  the  stretching  problem,  and  there- 
by essentially  matched  the  punch  stretching  component  with  the  pure  radial 
drawing  component  at  a particular  point  in  the  die  profile  region. 

Instead  of  this  tedious  process  of  boundary  matching,  it  is  desirable 
to  have  a numerically  efficient  and  reliable  method  which  can  treat  the 
problem  in  a unified  manner.  The  FEM  is  such  an  alternative.  The  finite- 
element  model  developed  for  the  deep  drawing  problem  is  the  one  by  Wifi 
[44]  with  a limited  treatment  of  friction.  Also,  Levy  et  al . [59]  developed 
the  elasto-plastic  finite-element  program  for  cupdrawing  based  on  the  defor- 
mation theory  of  plasticity. 

2.  Computational  procedure 

The  entire  sheet  undergoing  the  deep  drawing  process  can  be  divided 
into  four  regions:  the  contact  region  with  the  punch  head,  the  unsupported 
region,  the  contact  region  with  the  die  profile,  and  the  flange  over  the 
die.  Different  kinds  of  boundary  restrictions  are  imposed  depending  upon 
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the  regions.  For  example,  the  flange  is  constrained  to  move  only  horizon- 
tally along  the  die  face,  while  the  contact  region  with  the  die  profile  or 
punch  head  should  satisfy  the  kind  of  boundary  conditions  discussed  in 
Section  V. 

The  only  difference  in  deep  drawing  with  a hemispherical  head  punch 
from  the  punch  stretching  with  a round  die  corner  is  the  presence  of  the 
flange  which  is  free  to  slide  over  the  die.  The  addition  of  this  moving 
flange  is,  in  effect,  equivalent  to  the  addition  of  the  third  moving  boundary, 
because,  even  though  the  boundary  separating  the  flange  from  the  die  profile 
remains  stationary  in  the  space,  it  continues  to  move  from  the  viewpoint  of 
the  deforming  sheet.  To  treat  this  we  make  an  assumption  on  this  third 
moving  boundary  and  see  if  it  is  true  by  checking  the  radial  positions  of 
the  nodes.  If  the  new  radial  position  of  any  node  which  is  assumed  to  lie 
on  the  flange  or  the  die  profile  does  not  fall  on  the  expected  region  after 
converged  solution  is  obtained,  then  the  boundary  assumption  is  modified. 

Another  point  to  be  mentioned  is  the  blank  holding  condition  of  which 
there  are  two  types:  clearance  holding  and  force  holding.  The  idealization 
of  the  deformation  state  corresponding  to  the  force  blank  holding  is  the 
plane  stress  state  and  the  one  corresponding  to  the  clearance  holding  in 
the  plane-strain  state.  The  present  rigid-plastic  FEM  is  built  to  handle 
the  plane  stress  state  deformation  and  therefore  a modification  is  necessary 
to  handle  the  clearance  blank  holding.  No  reported  work  on  deep  drawing 
with  a hemispherical  head  punch  under  clearance  blank  holding  is  available 
and  therefore  in  the  present  work  only  the  deep  drawing  with  the  force 
holding  is  analyzed.  The  blank  holding  force  is  implemented  in  the  formula- 
tion as  a tangential  friction  force  acting  on  the  last  node  located  at  the 
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rim  of  the  sheet.  The  distribution  of  the  blank  holding  force  over  a finite 
area  near  the  rim  can  be  handled  without  difficulty  in  the  present  FEM,  but 
this  distributional  effect  turns  out  to  be  insignificant  [53] . Therefore, 
tangential  frictional  force  is  confined  to  the  last  node  at  the  rim  of  the 
sheet.  The  increment  of  deformation  is  controlled  by  the  punch  head  move- 
ment. The  program  is  in  Appendix  D. 

3.  Results  and  discussion 

The  only  available  work  on  the  complete  analysis  of  deep  drawing  with 
the  hemispherical  head  punch  is  one  by  Woo  [58] . Along  with  the  numerical 
solution  by  the  finite-difference  method,  he  also  conducted  an  experiment. 
The  parameters  are: 


Material:  soft  copper 

Stress-strain  characteristics:  a = 5.4  + 27. 8e0,504  ton/in." 

0 504  9 2 

for  e < 0.36:  = (0.08208  + 0.422569e  ) * 10  N/m 

c a -..0.375  , . 2 

= 5.4  + 24.4s  ton/in. 

for  £ > 0.36:  = (0.08208  + 0 . 37089£° ' -;>75)  * 109N/m2 

Blank  radius:  2.2  in.  = 5.588  x 10  “ m 

-2 

Radius  of  the  die  throat:  2.123  in.  = 5.392  x 10  m 

-2 

Radius  of  die  profile:  0.5  in.  = 1.27  x 10  m 

-2 

Radius  of  punch  head:  1 in.  = 2.54  x 10  m 
Blank  holding  force:  0.5  ton  = 500  kg 
The  solution  by  the  rigid-plastic  FEM  is  in  excellent  agreement  with 
the  experiment  for  the  flange  part;  however,  over  the  punch  head  it  predicts 
more  straining  than  the  experiment  when  the  friction  coefficient  of  0.04  is 
assigned  for  the  contact  region  over  the  punch  head  and  over  the  die  in  the 
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numerical  analysis.  When  the  friction  coefficient  is  increased  to  a value 
of  0.1  over  the  punch  head,  while  the  same  friction  coefficient  of  0.04  is 
used  for  the  flange,  the  analysis  predicts  less  straining  over  the  punch 
head  than  the  experiment.  See  Figs.  30,  31,  32,  and  33.  The  deviation  of 
the  numerical  solution  from  the  experimental  data  gets  larger  as  deformation 
progresses,  which  is  reflected  in  the  punch  load  vs.  punch  depth  relation- 
ship in  Fig.  34. 

The  lubricant  used  in  the  experiment  is  graphite  in  tallow  and  Woo  sug- 
gested the  friction  coefficient  to  be  0.04.  In  the  analysis  the  practical 
difficulty  always  lies  in  the  assignment  of  a reasonable  value  of  friction 
coefficient  because  friction  coefficient  under  a real  sheet-metal  forming 
condition  is  hard  to  measure  and  it  may  even  change  during  deformation. 

Comparison  of  Woo's  numerical  solution  with  the  experimental  data  does 
not  yield  any  better  agreement  than  the  present  rigid-plastic  FEM.  In  com- 
paring his  numerical  solution  with  the  experiment  Woo  made  the  correction 
on  the  circumferential  strain  based  upon  the  argument  that  the  strain  value 
obtained  from  the  analysis  is  the  value  at  the  neutral  surface  of  the  sheet, 
while  experimental  data  are  obtained  from  the  outside  surface  and  therefore 
a compensation  for  the  thickness  difference  is  necessary.  There  could  be 
a question  about  Woo's  correction  because  the  ratio  of  the  punch  radius 
or  die  profile  radius  to  the  sheet  thickness  is  sufficiently  large  in  his 
experiment  that  the  membrane  theory  is  justifiable.  Besides,  it  seems  a 
more  consistent  way  to  consider  the  problem  in  the  three-dimensional  stress 
state  instead  of  the  plane  stress  condition,  which  is  the  case  used  in  Woo's 
analysis,  if  the  variation  of  the  strain  across  the  thickness  is  to  be  taken 
into  account . 
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figure  31.  Distribution  of  Circumferential  Strain  for  y = 0 


■imt’erentiRl  .strain 


Figure  33.  Distribution  of  Circumferential  Strain  for  u = 0.1, 
Ud  = 0.04  P 
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Figure  34 . Punch  Load  vs . Punch  Depth 


It  is  necessary  to  have  more  numerical  solutions  and  experimental  data 
with  a known  friction  state  to  assess  the  validity  of  the  present  rigid- 
plastic  FEM  for  deep  drawing  problems.  However,  the  present  rigid-plastic 
FEM  had  dealt  with  other  sheet-metal  forming  problems  in  a unified  and 
consistent  manner  and  therefore  it  seems  reasonable  to  expect  its  validity 
for  deep  drawing  problems  when  it  is  established  for  other  problems. 
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SECTION  VII 


SUMMARY  AND  DISCUSSION 

It  has  been  made  clear  that  classical  variational  formulations  for  the 
rigid-plastic  solid  are  not  appropriate  for  solving  the  sheet-metal  forming 
problems.  This  is  due  to  the  nonuniqueness  of  the  deformation  mode  under 
certain  boundary  conditions.  This  nonuniqueness,  however,  can  be  resolved 
by  taking  the  workhardening  rate  into  consideration.  Such  an  introduction 
of  the  workhardening  rate  into  the  formulation,  on  the  other  hand,  necessi- 
tates the  consideration  on  the  geometry  change.  The  available  classical 
formulation  in  which  these  two  aspects  are  considered  is  not,  however, 
applicable  to  the  statically  indeterminate  problems,  sheet-metal  forming 
being  one,  because  it  is  formulated  in  such  a way  that  knowledge  of  stress 
distribution  is  necessary. 

Within  the  framework  of  Eulerian  descriptions  and  the  hypothetical 
identity  of  the  deformed  configuration  with  the  undeformed  configuration, 
further  improvement  in  the  applicability  of  the  variational  formulations 
to  the  statically  indeterminate  problems  is  not  possible.  Therefore,  an 
incremental  deformation  at  a generic  stage  is  considered  by  separating  the 
deformed  configuration  from  the  undeformed  configuration.  The  relevant 
equations  are  expressed  with  the  undeformed  configuration  at  each  step  as 
the  reference  frame  and  the  variational  formulation  is  established. 

From  this  variational  formulation  a finite-element  model  is  developed 
for  the  sheet-metal  forming  prlblems.  In  many  sheet-metal  forming  processes 
the  membrane  theory  is  justifiable  and  therefore  this  idealization  is 
introduced  in  building  the  model. 
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Three  basic  sheet-metal  forming  processes,  i.e.,  the  bulging  of  a 
sheet  subject  to  the  hydrostatic  pressure,  the  stretching  of  a sheet  with 
a hemispherical  head  punch,  and  deep  drawing  of  a sheet  with  a hemispherical 
head  punch  are  solved  by  the  proposed  method  and  its  solutions  are  compared 
with  the  existing  numerical  solutions  and  the  experimental  data.  The  agree- 
ment is  generally  excellent  and  therefore  the  prime  objective  of  the  present 
investigation  has  been  achieved. 

In  hydrostatic  bulging  the  strain  distributions  and  the  pressure  vs. 
polar  height  relationship  predicted  by  the  present  rigid-plastic  FEM  are 
in  excellent  agreement  with  the  available  numerical  solution  by  the  elasto- 
plastic  FEM  and  experimental  data.  The  difficulty  of  satisfying  the 
boundary  condition  along  the  fixed  periphery  experienced  in  the  finite- 
difference  method  does  not  appear  in  the  present  rigid-plastic  FEM. 

In  punch  stretching,  to  make  the  problem  more  tractable,  the  presence 
of  the  die  profile  is  neglected  first  so  that  there  is  only  one  moving 
boundary.  This  problem  is  successfully  solved.  Taking  the  die  profile 
into  consideration  is  equivalent  to  introducing  another  moving  boundary, 
and  while  handling  two  moving  boundaries  simultaneously  could  be  time 
consuming,  the  present  rigid-plastic  FEM  again  proves  to  be  efficient  and 
reliable.  The  strain  distributions  and  the  punch  load  vs.  punch  depth 


relationship  predicted  by  the  present  rigid-plastic  FEM  are  in  excellent 


agreement  with  the  numerical  solutions  by  the  finite-difference  method  and 


the  experimental  data. 

We  then  investigate  the  influence  of  workhardening  representation  by 
comparing  solutions,  computed  by  both  the  parabolic  workhardening  law  and 
the  Voce  equation  methods.  The  two  workhardening  representations  result 


81 


in  the  difference  of  peak  strains  and  load-displacement  relationships,  and 
the  difference  becomes  increasingly  significant  as  punch  displacement  in- 
creases. It  is  concluded,  however,  that  the  selection  of  a proper  work- 
hardening representation  requires  more  experiments  with  improved  accuracy 
and  control. 

The  present  method  is  further  extended  to  the  deep  drawing  problem. 

The  strain  distribution  predicted  by  the  present  rigid-plastic  FEM  is  in 
excellent  agreement  with  the  experimental  data  over  the  flange  of  the  sheet; 
however,  over  the  punch  head,  agreement  is  not  as  good.  By  assigning  two 
different  values  of  the  friction  coefficient  over  the  punch  head,  two 
strain  distributions  are  obtained;  one  predicts  more  straining  than  the 
experimental  data,  and  vice  versa.  Therefore,  an  improvement  in  the  predic- 
tion seems  possible  by  giving  the  friction  coefficient  a proper  value  which 
is  between  these  two  bounds;  however,  the  validity  of  the  present  rigid- 
plastic  FEM  for  deep  drawing  analysis  remains  inconclusive  at  this  stage 
mostly  because  of  the  lack  of  comparable  numerical  solutions  and  experi- 
mental data.  This  is  apparently  due  to  the  increased  sophistication  and 
accompanying  computation  time  when  three  moving  boundaries  are  treated 
simultaneously  and  to  the  practical  difficulty  of  determining  proper  fric- 
tion coefficients. 

It  is  concluded  that  the  present  rigid-plastic  FEM  can  treat  the  sheet- 
metal  forming  problems  with  efficiency  and  reasonable  accuracy. 
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APPENDIX  A 

PROGRAM  FOR  THE  INITIAL  GUESS  FOR  HYDROSTATIC  BULGING  ANALYSIS 

This  program  is  to  provide  the  initial  guess  and  initial  geometry 
for  Appendix  B.  It  is  based  upon  the  analysis  by  Hill  [23]. 

(I)  Data  preparation 

1.  Read  NUMNP  (IS) 

NUMNP:  Total  number  of  nodal  points  to  be  generated 

2.  Read  RADIUS,  DIS1,  DIS2  (3F  10.0) 

RADIUS:  Radius  of  the  sheet  to  be  bulged 

DIS1:  Polar  height  of  the  bulge  in  the  initial  geometry 
DIS2:  Polar  height  of  the  bulge  in  the  new  configuration 


\ 


( 


GPIO 
GO  IO 
G»T  C 
GOTO 
GO  IO 
GOTO 
G°  T O 
GO  I 0 
cep 

GOTO 
GOTO 
GO|0 
GO  I C 
GOTO 
GOTO 
GOI  0 
GOTO 
GOTO 
GO  TO 
GO  I 0 
GOTO 
GOTO 
GOTO 
GOI  0 
GRID 
GOTO 
GO  10 
GOTO 
GOTO 
GRID 
GO  I 0 
coin 
GRID 
GOTO 
GOTO 
GOTO 


1 PROGRAM  GO IlM  INPUT .OUTPUT ,T  AOS5* INPUT ,TAOc6*0UTOUT .PUNCH » 

2 C 

3 C****» A************************** **•***•*•***********»*•****-***"** 

A C THIS  PROGRAM  IS  tq  GENERATE  TFF  INITIAL  0E0MET  PV  ANC  VELOCITY 

5 C FIELD  TOR  HVOOPST/TIC  8LLGR  PROBLEM,  FOLLOWING  HILL 

6 C*********** *************  A A A* A ******* **** •♦»»»****»** A*** A* ******** 

7 C 

A COMMON  *(2000) 

9 C 

1 0 C 

1 1 <?*•****»-*****•** ******* a* *******•••*< ******** »*«**.* *»*»»*** *>*»*«• 

12  C NUHNOsNUMBEP  OF  NCOAL  POINTS  TC  of  CENEOATFC 

I?  f 
1 4 C 
15 

1 6 C 
IT  r 
1 8 
1 O 
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24 

25 

26 
PT 
?« 

29  C 

30  C 

?1  CALL  GtJEOSfAfNl  > . A ( N2  1 . A ( NT  ) , A C N4  > . A < US  ) * A C N6 1 , A ( N7  > , A C N « ) . A ( NC  J , 

22  1 A(  MO!  .NUMNO) 

33  C 

T4  1001  F0QMATIT5) 

35  STO° 

36  ENC 


GO  ID 

38 

SUROOUT INE  GUESS! OR .77 • CODE . SLOP , P . Z . UP , U 2 .UUP ,UU7. NU*F c \ 

GOT  0 

30 
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GOTO 

44 
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GOTO 

45 

K=  0 

GOTO 

46 

CI?sR  IS  1 

GOTO 

47 

OOaOADI US/FLOAT ( NUMNC-l ) 

GO  ID 

48 
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40 

c 

GO  10 

50 

50 

KsK+1 

GOf  0 

51 

IFCK  .FO.  21015*0152 

GOI  O 

5? 
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GOTO 

53 

DC  100  IM.NUWNF 

GOTO 

54 

T F ( IT  .FC.  1 »c  ( I )*PO 

GOTO 

55 
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GOTO 

56 
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57 
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62 
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APPENDIX  B 


PROGRAM  FOR  THE  ANALYSIS  OF  HYDROSTATIC  BULGING 


This  program  is  for  the  analysis  of  hydrostatic  bulging. 


(I)  Data  preparation 

1.  Read  HED  (A  12) 

Output  title 

2.  Read  RVALUE,  T,  ACOEF  (5F  10.0) 

RVALUE:  Normal  anisotropy  parameter 

Set  1.0  for  isotropic  material 

T:  Initial  thickness  of  blank 

ACOEF:  Accelerating  coefficient 
To  start  with,  set  1.0 

3.  Read  ITER,  NREAD,  ITCONT,  NFORM,  NPUNCH,  NPRINT , FLIMIT  (615,  F 10.0) 
The  program  control  card 

ITER:  Number  of  iterations  to  be  executed 

NREAD:  1,  if  new  data  are  to  be  supplied; 

0,  otherwise 

ITCONT:  0,  if  computation  starts  at  the  very  beginning  and  first/ 
second  steps  are  included  in  the  steps  to  be  computed; 

1,  otherwise 

NFORM:  Number  of  steps  to  be  computed 

NPUNCH:  1,  if  solution  is  to  be  punched  at  the  end  of  each  step; 

0,  otherwise 

FLIMIT:  Value  of  (error  norm) / (solution  norm)  required  for 
convergence.  To  start  with,  set  this  .000001 

4.  Read  NUMNP  (6  I 5) 

NUMNP:  Number  of  nodal  points 

5.  Read  YVALUE,  PRESTN,  EXPNT,  PRESTS  (4F  10.0) 

Material  characteristics  are  specified. 

Stress  = YVALUE*  (Strain  + PRESTN) **EXPNT  + PRESTS 
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6.  Read  PRES,  DPRES  (4F  10.0) 

PRES:  Current  pressure  value 
DPRES:  Increment  of  pressure 

7.  Read  N,  CODE(N),  R (N) , Z(N),  UR(N),  UZ(N),  SLOP(N),  (15,  F5.0,  5F  10.0) 
Nodal  information 

N:  Node  number.  Node  number  1 is  at  the  rim  of  the  blank  and  the 
last  node  is  at  the  pole 

R(N) : Radial  position  of  the  node 

Z(N):  Axial  position  of  the  node 

UR(N) : Increment  of  displacement  in  radial  direction 

UZ(N):  Increment  of  displacement  in  axial  direction 

SLOP(N):  Slope  of  the  element 
Set  this  0.0 

C0DE(N):  Type  of  boundary  conditions: 

1.0,  if  magnitude  of  UR(N)  is  fixed; 

2.0,  if  magnitude  of  UZ(N)  is  fixed; 

3.0,  if  magnitudes  of  UR(N)  and  UZ(N)  are  fixed; 

0.0,  if  neither  the  magnitude  of  UR(N)  nor  UZ(N)  are  fixed 

In  subroutine  PRELIM  the  interpolation  of  data  is  built  in. 

8.  If  NREAD  = 1,  the  input  data  is  to  be  placed  behind  nodal  information 
cards 
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"ULGF 

9?6 

400 

BIN l*B(N)-A(M)4P( | ) 

PULGE 

927 

GO  TP  300 

9ULGE 

920 

c 

PULGE 

9 29 

ENG 

BULGE 
PULGF 
BULGE 
eUL  Cf 
BULGF 
BULGE 
PUL  GF 
BULGE 
0V.1L  GF 
PULGF 
BULGE 
BULGE 


| 

944  SUBROUTINE  HAPD?(EDS«V| 

945  C 

946  C44**4*4**444*4*4*******#*******4***********4*i 

947  C COMPUTE  WORK  HARDENING  BATE 

940  C4*****5********4*»4444A»**4445**4*#*******44*i 

949  C 

950  COMMON/M AT EPL/VV ALU E, PO* STN, F XFNTt OPE STS 

951  C 

05?  Y*E*OMTAVVALUE4(ODE«TN4EPSI44(F  >pmt-1, 1 

953  C 

954  PET  L<>N 

955  ENC 


BULGE 
PULGF 
BULGE 
BUL  GF 
PULGE 
BULGE 
PULGE 
BULGE 
PULGE 
BULGE 
BULGE 
BULGE 


931  SUBROUTINE  HARC(EFS,V| 

93?  C 

933  C4»**4**4*44*4 5*4 4444 444***44**44 4*44 44*444444: 

934  C WORKNAPOFNTNG  CHAPACTEPI STIC  CLP VE 

935  C4444444************** 4 4*4 444*44* 4*4* 4* ******* 

936  C 

937  COBMHN/MATFPL /V  VALUE .PPESTN, E XFNT . OPE  STS 
93B  C 

939  7 *V  VALUE  * ( PR5£TK*EdE145EXPN'T*P9FSTS 

940  C 

941  RETURN 

94?  ENC 
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APPENDIX  C 


PROGRAM  FOR  THE  ANALYSIS  OF  PUNCH  STRETCHING 


This  program  is  for  the  analysis  of  the  stretching  of  a sheet  with 
hemispherical  punch,  where  the  die  profile  is  neglected. 


(I)  Data  card  preparation 

1.  Read  HED  (A  12) 

2.  Read  RVALUE,  T,  ACOEF  (5F  10.0) 

3.  Read  ITER,  NREAD , ITCONT,  NFORM,  NPUNCH,  NPRINT,  FLIMIT  (615,  F10.0) 

4.  Read  NUMNP  (615) 

5.  Read  PNRAD,  RADIUS,  FRITN  (4F  10.0) 

PNRAD:  Radius  of  the  hemispherical  punch 
RADIUS:  Radius  of  the  blank 

FRITN:  Friction  coefficient  between  the  punch  head  and  the  blank 

6.  Read  YVALUE,  PRESTN,  EXPNT,  PRESTS  (4F  10.0) 

7.  Read  ECONST,  TDIST  (4F  10.0) 

ECONST:  Step  size  in  terms  of  the  maximum  magnitude  of  the  effective 
strain  increment.  To  start  with,  set  this  0.04 

TDIST:  Criterion  distance  of  the  contact  of  the  sheet  with  the 
punch  head.  To  start  with,  set  this  0.008 

8.  Read  N,  CODE(N) , R(N) , Z(N),  UR(N) , UZ(N),  SL0P(N),  (15,  F5.0,  5F  10.0) 
Code  (N)  =4.0  for  the  contact  zone  of  the  sheet  with  the  punch  head 

9.  If  NREAD  = 1,  the  new  input  data  is  to  be  placed  behind  the  nodal 
information  card 
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STOCK 

1 

STB  CM 

2 

STOCK 

3 

STPCH 

4 

STOOh 

s 

STPCH 

ft 

STPCH 

7 

STPCH 

A 

STPCH 

0 

STRCw 

10 

STB  CM 

1 1 

STPCH 

1 2 

STPCH 

1 3 

STRCW 

14 

STPCH 

IS 

STRCw 

1ft 

STPCH 

1 7 

STOCK 

1 8 

STPCH 

I o 

STPCH 

20 

stoch 

21 

stock 

2? 

STPCH 

23 

STPCH 

24 

STo^h 

2* 

STBCH 

?ft 

STOCh 

?7 

STPCH 

2° 

STOCH 

2C 

STOCh 

30 

STRCw 

31 

STCCw 

32 

STRCw 

3 3 

STOCH 

34 

STPCH 

3* 

STPCH 

36 

STPCH 

37 

STPCH 

38 

STRCH 

3P 

STPCH 

40 

STOCH 

4 1 

STB  CM 

4? 

STRCw 

43 

STOCH 

44 

STOCH 

45 

STOCH 

4* 

STPCH 

4? 

STOCH 

4P 

STOCH 

49 

STRCH 

so 

STPCH 

51 

STOCH 

52 

STOCH 

53 

STPCH 

SA 

STPCH 

55 

STOCH 

56 

STOCH 

57 

STPCH 

58 

STPCH 

50 

STOCH 

ftO 

STPCH 

61 

ST®CH 

6? 

STOCH 

6.1 

STPCH 

64 

STPCH 

65 

STRCH 

6* 

STPCH 

67 

STOCH 

66 

STOCH 

60 

STOCH 

70 

STPCH 

71 

STOCH 

7? 

STOCH 

73 

stoch 

74. 

STOCH 

75 

STPCH 

76 

STPCH 

77 

c 

C! 


c 

e> 

c 

c 

r 

r 

C> 

c 


c 
r i 
C 
C 

o 

c 

c 


COCCt»«  STtCH  INPUT  • "*UTCUT  ,TAPf  «r  lNruT,T*or  *i:t'TOUT 


Tw*I  c 

HPR  F 


PROGRAM  JC  Cf-3  T^c 

T Mc  RAOJUS  PF  THt 


ANALVi!«  C*P  ’ HF  0|»NCW  ST3cTf^!W(.,  PYJ.KJM 
OIc.  CCfFTLc  ! S NFrt_crT5C. 


CP«MPN/rcNriN  /NLMMP  ,Nl)KEL  ,h?o(  1 ?)  .DLL  .NEC  ,Vc°BW  • v I Flr  , Tc  SI  t ITER  . 
1 NOE  AD.  NPUNC-f  KTP  TNT  ,C  VALUE  ,T  , a»RAND.  PNP  AC,  OAD  »US  , I TN  , 

PFC3NST ,PNME D . TP1 ST 

Civ  M"’  K/MATE  PL  / vViLUc*P®ScTM,e*PNTtD~[  CTS 
C^MMOK/J SO/PRVAl 


OPTS3AM  IS  rpc  CCNTT  DLL  I Nr'  THf  OlMFNSTCN  rp  TK  CD^OLcTE 

DOPfQAM,  ITS  PlIODre?  IS  TO  P5T  VENT  ASS’GNJNG  A C ABGE  0 THAN 
NECESSARY  DIMENSION  POP  ANY  ACCAY  THCOUC^  r »r  USE  HP  -rw  = 
p OLID*  I NG  STATFMEn-t 


COMMON  a < 5C00  ) 


NPIELT  IS  THE  ClMrN'STCN  QP  AOCAY  A.  ITS  VALUr  CAN  PE  DPT£SMlN!.n 
®P  :C  I SCL  Y RV  PUNNING  THF  PP^GCAM  PnCT  . 


NF  I EL DrEOOC 


«r  tr  • r * 


c 


e 

c 

c 

c 

c 

c 

c 

r 

c 

c 

r 

C 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 


TEST* l . 
PPADfG.l 000) 
PE  AD( 5 « I 004 ) 
PE  AC(S,  1003  > 
PE  ADf  s, 1 0 C 3 ) 
PEACf  E , 1 004  ) 
BE  ACCS  ,1 004  > 
READf  S , 1 0 04  > 


PEC 

rvalue  ,“r • acoep 

TTER.NBEAC,  ITCONT  , NP-TCA  ,N=UN  CH  or  I NT,  pl  I v;  T 
NUMnP 

PNBAC,CADILS,RRI'rN 
W ALU E . PRESTN. £x CNT , ORESTS 

cccnst . tc i st 


HE0*PUTPUT  TTTLE 

RVALUE = VALUE  CE  THF  AKISrTPCPV  CtttwrTSC 

ACTPF  = ACCELET  AT  INC  OB  C c C r.L  E = A T I N C !C  I ENT  cp  C CK'VL  OGENCS 

NCFAOsO,  IP  TC  PYCACC  thf  CSATING  STATPMPNT  IN  SU°  = CUTI'ic  PL  APT 
ITCONTrC.  TP  CDMOLTATICN  S7ACTS  at  THE  VF  = V »£GG1NM6G  ANC  PI^ST/ 
SECCNC  STFCS  APE  tNCLUOEn  IN  Thp  3 T F o 5 Tr  CCMOUTEC 

* 1 • 0 THE  P * I SF 

THIS  INCH*  I s PELATE0  Tr  t*-f  rFT EC m I NA T I ON  DP  STE°  SI  ?E 
NPDCV*NU«*0cR  CP  STEP®  ASSIGNEr  pun 

NPUNC  H~  I , IP  C A 7 A ARP  TC  5E  PINCHED 
*0,  CTFFP»»ISF 

PL  IMII*  VALUE  CP  (ERRCF  NCPM1 / f SCLUT ION  NC*om » RFCUICFP 
POP  CCNVEPGENCE 

NPPfNTsl  , IP  N 00  A U PCINT  OA-A  /CP  TC  PE  PRICED 

*C,  otnepmsf 

NUMNCsNUMPER  CP  NODAL  POINTS 
ONBADaBAOIUS  CP  H? M ? epnro i c AL  PUNCH  h£  a n 
OAC tUS*P AD IUS  OP  THE  RLANK 

rptTN  *PBICTTCF  CCERPICIEN"^  “E’wEFN  PLAN*  AND  PUNCH 
cCONST*  STEP  SI?E  IN  MOIMUM  PPPECTIVE  STCAIN  INCREMENT 


YVALUE,  PRESTN,  ryPNT,  ORESTS  ApF  TO  EX°°eSS  THE  vpC KMA POEN I NG 
CH  A C A C Tc  PI S T I C S CP  THE  PLAN* 

STRESS*YVALUc*(OPFSTN»STRA I N 1 a *c YPNT4 GREETS 

NEOeNUHPEC  OP  FCUATICNS  TC  «F  SOLVEC 
NUMEL*NUKPSR  CP  ELEMENTS 
*#BAN0*PAN0  • I C T h 
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STPCH 

T*  C 

STBCm 

BPV  4L  «PV4LUE 

STPCH 

1 

SO 

KU*SL  »KU*NP-1 

STOfH 

1 

»1 

MMANOa* 

STPCH 

( 

1? 

K50«NU*4KP*3 

STPCH 

( 

> 3 

NOiNEQ 

STPCH 

( 

14 

N5L *NUMgL 

STPCH 

1 

If 

C 

STPCH 

1 

S6 

c 

STPCH 

t 

IT 

N 1 * 1 

STPCH 

( 

se 

N2«N  1 ♦NU*4NP 

STOCH 

( 

IO 

N3*N2  4NUMNP 

STPCH 

>n 

N4bK3*KUMNP 

STPCH 

« 

N**N44-NU*4NP 

STPCH 

>? 

K6  *E54NUMNP 

STPCH 

< 

>3 

K7«K6*NtMNF 

STPCH 

«; 

>4 

N%«M7«NU*4EL 

STPCH 

\ 

J5 

*4P*KS»KVJPEL 

STPCH 

>0 

N 1 0«N<J*NLMEL 

STPCH 

< 

>7 

M 1 -MO  4NUMEL 

STPCH 

< 

JS 

NI2«M 14KUPEL 

STPCH 

»o 

N1 3*N 124NUMEL*? 

STPCH 

i< 

)0 

M4«N1  34KU»SL*4 

STPCH 

i ( 

)1 

N 1 3*N 144NUMEL44 

STPCH 

1 c 

>2 

M6*N154NEC 

STPCH 

l < 

53 

M7«N154KEQ4MF4KC 

STPCH 

i< 

>4 

N19*N 1 74NUMEL 

STPCH 

1C 

>5 

Nl  <3»N1  94NUMNP 

STPCH 

i< 

>6 

N?0>N 1 5*NIMNP 

STPCH 

it 

>7 

N21 *N20 ♦NUPNP 

STPCH 

i ( 

58 

K22  *N21  4NU»*KP 

STPCH 

it 

>o 

N2?«N2?4NUHNO 

STPCH 

1 1 

! 0 

K24«N?34NcC 

STPCH 

1 1 

1 1 

N 25 >N?44NUMNP 

STPCH 

1 1 

2 

N26*N254NUMNP 

STPCH 

1 1 

13 

*27  >N26+ NUMKP 

STPCH 

1 1 

14 

c 

STPCH 

i i 

5 

c 

STRCH 

1 1 

16 

c 

STPCH 

1 1 

7 

IPCN27  ,L  E • NFIEL01G0  TO  100 

STPCH 

1 1 

i e 

C 

STPCH 

1 1 

ip 

WP I TF t 6 , 1 00  1 ) N?7 

STPCH 

1 c 

•0 

ST  C P 

STPCH 

i ; 

?i 

C 

STPCH 

l < 

*? 

100  CONTINUE 

STPCH 

i ; 

*3 

■PITS (6 .J 002 » K?7 

STPCH 

u 

>4 

C 

STPCH 

i < 

►5 

r 

STPCH 

i ; 

>6 

CALL  CDEL  I ••  ( A IN  1 1 * A f N2  1 • A ^ N3  1 « A ( NA  1 « 4 t NS  ) * A 1 Nf  ) ) 

STPCH 

l ; 

>7 

C 

STOCH 

1 1 

» 8 

CALL  PLASH  AtM  I*  Af  NT)  , 4|  N3>,  / CN4  »,  ACN5  1 , 4CN*»  , A(N7  1 %ACN9  1 ,A  CNP1  • 

STPCH 

i ; 

>0 

lA(N10),AtNll),A(M2».AtM3),A(M4»,A(Nl5>.A(N16),A<M7).A<Nl«). 

STPCH 

i : 

JO 

2A(N1Q>#  A(N20),A(N?1  ) , A < N2?  ) . A t N 23  ) , A t N24 ) , A ( N2 * ) , A ( N2 6 ) 

STPCH 

i : 

n 

3 v NC.N5L*cLIMIT#ITCCNT.ACCEF) 

STPCH 

i : 

«2 

c 

STPCH 

l : 

S3 

1000  FCBHATC12A6) 

STPCH 

i: 

)4 

1001  C1°M*T(  ✓ //*  THE  DIMENSION  no  THE  APCAY  (4)  IS  TCC  ?4*Al.L*/ 

CTPCM 

i : 

15 

1*  TEE  SIZE  OP  TEF  APP4Y  (A)  ML  ST  PE  *.  17) 

STPCH 

i : 

J6 

1002  ^OPMATt//*  THE  NECCESSABV  SIZE  C'  TE-f  AOray  (A)  IS*.  17) 

STPCH 

12 

J7 

1003  FOBMATtfcl5,FlC.C) 

STPCH 

i : 

J 9 

1004  FCCE4T(4P10.0) 

STOCH 

1 2 

•9 

1005  FrSMAT(4I5.FlC.C) 

STPCH 

14 

lO 

C 

STPCH 

14 

u 

$TQP 

STDCH 

14 

12 

ENC 

STPCH 

1 

14 

SUP POUT  I me  P PEL I*f c .7. LB • U 7 . CCOE . SLOP ) 

STPCH 

1 

15 

C 

STPCH 

1 

16 

STOCH 

1 

!** 

C PE  40  ANO  POINT  0C  CONTROL  IKFrrUATlCN  AND  M A Tc  B I A L PPOBEPTIES  MATN0013 

STPCH 

1 

19 

STPCH 

1 

IP 

C 

STPCH 

50 

C 

STOCH 

51 

COMhoE/GENCOK/NUHNP.nueEL. H£C  tl 2 ) . DLL • NFO , N*C PM , v i EL  0 • TE ST , I TEC . 

STOCH 

52 

1 NPcAn.MOUNCH,  NPP  INT.C  VALtJE  . T,MPANO  , BNP AO,  PAOlUS.PPITN, 

STPCH 

53 

2ECCNST, CMMEC.TC  1ST 

STPCH 

54 

C 

STPCH 

55 

DIMENSION  P ( 1 ) . ?( 1 ) .COCEt 1 ) . UP t 1 ) .U7C 1 ) , SLOP! 1 ) 

STPCH 

54 

C 

SToch 

| 

57 

c 

STOCH 

59 

50  CONTINUE 

K 


p 


i 


I 


STPTH 

1*3 

W=>I  tp  (6.2000  HFr  tKuvKC,K.jv«-L 

STBCH 

1*0 

CALL  HAPCf C .. V IFir  1 

S^SCH 

1 6 1 

WP I TE ( 6 ,201 O ) vlEir 

stpcw 

1 *2 

I 3c  ( f>  , i 000  ) ITFP 

ST®  CH 

163 

r 

«;rorH 

1*4 

r***v 

**  * 

4 *M  A N 0 0 3 0 

STOCH 

16* 

C 

R ■ A ” AN  r PRINT  0*  4*00  A L °PINT  D A ▼ A 

MAI  NO  031 

STPCH 

1 66 

f#t*t 

*** 

* »ma  JN0C2? 

STDCH 

1*7 

r 

CTOCM 

IfP 

L * 0 

MA 1N0C  34 

S’PCH 

l*c 

IP ( KPCI NT. “C. C 1 GC  V“  60 

STSC* 

170 

WRIT*  (6.11141 

STPCH 

171 

»oi Tc  (6  ,2004  ) 

ma  11,0033 

STPCH 

17? 

60 

READ  ( 5.10021  N ,CrPr  ( N ) , c ( N)  , 3 ( N » 

,UF(N) ,U?(N> ,SLC*(N> 

STPCH 

1 73 

C 

STPCh 

174 

NL  = L ♦ 1 

ms  tnoo?6 

STPCH 

17* 

ZX*N-L 

MAT  NO  0 37 

STPCH 

1 76 

T*(L  .FC.  0)  CC  y~  70 

STP  Ch 

1 73 

no: ( p (N  l-P (LI  1 /? x 

m a ! NO 0 39 

STPCH 

1 7fi 

07*  ( 7 < N )-Z(L  ) >/  7X 

MA ! NOC  3C 

STPCH 

1 7 R 

n u°  * ( Uc  ( N )-UF (L  11/7* 

STPCH 

1 *0 

fVJZ=(U7(*4  1-UZ(L  ll/ZX 

STPCH 

1 *1 

OS*  (SLCF(M-SLr.P(L  1 >/7x 

STPCH 

16? 

c 

STPCH 

1 9? 

30 

L*L  ♦ 1 

M A ! NO  04 1 

STP  CH 

1 64 

IF(N-L)  130,  PC, 9 0 

ma  TN004  2 

STRCH 

1 65 

C 

STPCH 

1 Sfc 

90 

CPOF (L ) *0.0 

MA I NO  04  3 

STPCH 

1 S7 

3 ( L 1 *P ( L — 1 )*OR 

M A I NO  0 4 4 

STPCH 

1 09 

SLOP ( L )=SLCP( L-  1 14PS 

STPCH 

1 *<3 

7(L 1 -71 L- 1 14P7 

STPCH 

1 00 

U9(  L ) *'JC  (L-  1 1 *0  0° 

STPCH 

l <51 

U7( L 1 *U 7 (L-l 1 4ru7 

STPCh 

1 0? 

GO  70  3 c 

M A IN004C 

STPCH 

1 03 

c 

STPCH 

1 04 

oo 

JF(NLVNP-N)  100,110.60 

STPCH 

1 os 

c 

STPCH 

1 06 

1 00 

WPITF  (6.2000)  N 

mat N0062 

STOCH 

1ST 

CALL  FXIT 

M a IN  00  5 3 

STPCH 

l 09 

c 

STPCH 

1 SO 

1 1 0 

CONTI NUF 

MA 1 N 0 0 54 

STPCH 

POO 

I*( nop INT.FO. Cl  GO  TO  120 

STPCH 

poi 

w°ITE  (6,2003)  (K , CO^F (K  ) , O ( K ) , 7( k 

),UP ( K ) ,U7( K l.SLHPfK  ) ,*  * 1 .NUMN3  » 

STOCH 

202 

1 20 

cont: nu* 

STPCH 

203 

c 

C TROW 

?0A 

N^C  *3*NL'MNP 

STPCH 

20* 

W 3 I T* (6,1122)  N£0,voand 

STOCH 

2 06 

c 

STPCH 

207 

1002 

PQPVAT  ( 15,55.0,5*1 0. Cl 

MA  IN  0 1 22 

STPCH 

209 

1 0 03 

FQCMiT ( 1615 ) 

STPCH 

200 

1 004 

*CPVAT( 19,211 1,2*10.9) 

STPCH 

210 

1005 

rODVAT(2I*,4*  10.0 

ST»CH 

21  l 

1 00* 

Fp C«AT(//  * TFF  NrpAL  PCINTS  AT  WHICH 

*nsre  CALCUL4'r!r7'c 

AF  F n*  s 

13 

STPCH 

21  2 

1 EO*  //  2 r I 5 ) 

STPCH 

213 

1 003 

FOFMAT  ( 1 HI  , 1 5X  , 3RH  LINEARLY  r*  I S3C  l PUT5C  p^UNOAOV 

STRESS* 

s/ 

STPCH 

>14 

1 / 60H  NCDE  I NODE  J ®F*SSUF* 

1 PRESSURE 

J 

SHE  A* 

T 

STOCH 

?!  S 

2 14H  SHF  A O J) 

STPCH 

216 

1 OOP 

FCFMAT  (210, 4E  1*.*  ) 

STPCH 

217 

1 noc 

c-pcwAT(///*  MAX  I ML  V NUWP*P  0*  JT-CATir  NS  ALLOWED 

*r*  EACH 

INCPEM 

rN 

STPCH 

219 

IT  = * , I?) 

STPCH 

210 

1114 

* OP  V A T ( J H 1 , 3 OH  NOPAL  F C I NT  INFCFAATICN  « * F-O  p £ SCALING//) 

STOCH 

220 

1 1 22 

Ffl PMAT(  / //  * NUV0Ec  PF  5 QUA  7 1 rs  e 

SB  , 

14/ 

STOCH 

221 

l * FANCwIOtf 

X * 

. 14/ 

STOCH 

22? 

2 * CT AGONAL  ELEMENT* 

X* 

• 14  ) 

STRCH 

223 

200C 

*OPMAT  (1H  1 2 A C / 

STPCH 

2 24 

1 30H0  NUMBER  C*  NCOAL  FCIN’S 

I 7 

/ 

MA  INO  1 27 

STOCH 

225 

2 3 OHO  NUMBER  C*  ELF  WENT  S — 

13 

/) 

STOCH 

225 

2002 

FpoMAT  ( I 12.* 12.2. 2*1 2 .3,  3= 24.7 ) 

STPCH 

227 

2003 

FOCMAT  (1113,416,1112) 

MA IN  0 1 27 

STPCH 

229 

2004 

FpOMAT  (/  * NODAL  POINT  TYPE 

P-PPOINATE 

Z-OPCTNATF  c 

L C 

STOCH 

229 

1 Ap  fo  oisdlacfmfnt  7 lCAD  CP  displacement 

B * 7 6 

-SLOPE 

* ) 

STPCH 

230 

200* 

cnBMA*(  //,»*PRCES  ?o*CI*I*r'  AT  NOPAL 

POINT* ,// . 

STOCH 

231 

| * NCOAL  PT • FLFMENTl  ELEMENT? 

PR*  SSU*e 

CMCAP* , /) 

STP  CH 

232 

2000 

* P°  MAT  (26H0NPPAL  PCINT  CARO  FP*rF 

Nx 

I*  ) 

M A INC  1A* 

STRCH 

233 

20  10 

FORMAT  ( / / * INITIAL  Y I*LD  STRESS  * 

* , 

* 15.7//) 

STPCH 

234 

RFTLCN 

STPCH 

23* 
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IL  «N 

STOCH 

1 267 

NsN-l 

STOCH 

1 26e 

IF(N*E0.01  RETURN 

STOCH 

1 260 

IH  = VI  NO  ( NN , N + VMM  1 

STOCH 

1270 

VsN 

QT-DCM 

1 27  1 

C 

cto-CH 

1 27? 

CO  400  I*lL  » 1^ 

«TOCM 

127? 

MsM ANN 

STOCH 

1274 

400 

9(N)=“(M-A(V1*F(  I> 

STOCH 

1 27? 

r 

STOCH 

1 ?7* 

GO  TO  30  C 

STOCH 

1 277 

r 

STOCH 

1 276 

ENO 

STOCH 

1 260 

SU9  POUT  IKE  PCt*IX(A,°,NEC,MP  4K'Cf  N,  CONST  A , CONST®.  CON  ST  F,  COK  STG  > 

C T P C H 

1 ?si 

C 

STP  OH 

1 ?6? 

r< 

**  « 

STOCH 

1 293 

r 

THIS  SUBROUTINE  HANDLES  TW'  vjxPn  5CUNDAPV  CONDITION 

STRCH 

l ?«4 

C***tty»t4*»»M**iMtHtt»»t****t**t»*>'A«**(*’>*,'***>»***t>t,,tt*lt>*:*>* 

r « a 

STOCH 

1265 

C 

STOCH 

l ?P6 

DIMENSION  A ( NFO  * U *P(  11 

STOCW 

1 2 97 

IR*l*N-2 

STOCH 

1 26® 

17=10* 1 

STP  CH 

1 280 

A ( 17,1  1 = 4(17*1  »♦?  **4f  IP,2l/rnN?TA*4(rc,  1)  /CHNC TA  /CONSTA 

STOCH 

1 200 

5(171  =5 (I  7>*6<IP)/C  CNST  4— CONS  T C 4 ( 4 ( JC*  1 >/C0NS-A*A(  TP , 2 ) > 

STOCH 

1 20  1 

0 ( I 7 1 * 5 ( I7>*CCNETE 

S TP  CH 

1 ?P? 

A < I F,  1 ) r l .0 

STOCH 

1 203 

A(  T S, 21  =0. 0 

STOCH 

1 204 

P ( TP  1*0  .0 

STOCH 

I 295 

c 

STOCH 

1 206 

DO  100  m*3,v«an0 

STPCH 

12  97 

Vs  I 7-VA  J 

STP  CH 

1 209 

I®(K  .LT.  1 ) GC  TO  100 

STPCH 

1 299 

UWtU. 1 

STOCH 

1300 

A(K,M>=A(K,V)*A  ( K « V V 1 /C  CNST  A 

STPCH 

1 301 

5(K)*«(Kl-A(Kfvvl*CCNSTG 

STPCH 

1 302 

1 00 

AIK  ,MM  > sO  ,0 

STOCH 

1 303 

c 

STOCH 

1 304 

MP  1 xM  5 ANIV  1 

CTOCH 

1 -*05 

c 

STOCH 

1 306 

DO  200  Vsg.MPl 

STPCH 

1 3 07 

V Vs  V ♦ 1 

STPCH 

1 ?0P 

A(I7,*0«A(IZ*V)*A(1  P. VV 1 /CONST  A 

STOCH 

1 309 

I 7 7 * I P *M 

STPCH 

1 310 

P(I?71*S(I7Z)-A(IC.vm ) 4 CONST  G 

STPCH 

131  1 

200 

4 ( 10.MM1 *0*C 

STPCH 

1 312 

c 

STPCH 

1313 

oettON 

STOCH 

1314 

END 

STOCW  1316 


SU°  POUT  I NE  HACDIEPf  .VI 
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c 

C' 

c 


«ORKHAOnfN!Nr<  CHAPA  rTE  s I 5T’ C CUt-Vf 


STPC* 

1 TIT 

STPCH 

1 31  P 

STD  CM 

1 MR 

STFCM 

1 320 

STPCH 

1 ■*  2 l 

STPCH 

1 ??2 

STPCH 

1 ??3 

STRCH 

1 7 24 

STPCH 

1 325 

STPCH 

13  2C 

STC  CM 

1 3?  T 

<-"MMnN/M*TEKL/VVALl''  iPPESTN#f  XPKT  • °FE  £T  S 

V*  V VA  L UF  * ( PR*?  «TNa€PSI**?KFA'T*PPE  STS 

PF  TURN 
ENC 


STPCH 

1 329 

SUBROUTINE  HAPP2(EPe.VI 

CTCCH 

1 330 

C 

STPCH 

1 331 

r *** 

STPCH 

1 3*2 

C 

COMPUTE  WORK  HAPPENING  PAT- 

STPCH 

1 33? 

c**< 

STPCH 

1 3 34 

r 

STPCH 

1 335 

COMMON /M ATE PL/ vv ALU E*PPF ST *CN T *PCE  ST  S 

STPCH 

1336 

C 

CTR  CH 

1 337 

V*S xoNT*v VALUE * C PPFSTNAF°S1**(F  »OKT-l  • 1 

STPCH 

1 33S 

c 

STPCH 

1 339 

RE  TURN 

CTPCH 

1 340 

ENC 
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APPENDIX  D 


PROGRAM  FOR  THE  ANALYSIS  OF  DEEP  DRAWING  AND 
PUNCH  STRETCHING  WITH  ROUND  DIE  CORNER 

This  program  is  for  the  analysis  of  deep  drawing  with  a hemispherical 
punch  head  and  stretching  with  a hemispherical  punch  head.  In  stretching, 
a round  die  profile  is  considered. 

(I)  Data  preparation  card 

1.  Read  HED  (A  12) 

2.  Read  RVALUE,  T,  ACOEF  (5F  10.0) 

3.  Read  ITER,  NREAD,  ITCONT,  NFORM,  NPUNCH,  NPRINT,  FLIMIT  (615,  F10.0) 

4.  Read  NUMNP,  NDEX  (615) 

NDEX:  2,  if  punch  stretching  is  to  be  analyzed 
3,  if  deep  drawing  is  to  be  analyzed 

5.  Read  PNRAD,  RADIUS,  DIERAD,  RTART  (4F  10.0) 

DIERAD:  Radius  of  the  die  profile 

RTHRT:  Distance  from  the  pole  to  the  die  throat 

6.  Read  FRITNP,  FRITND,  BHFCE  (4F  10.0) 

FRITNP:  Friction  coefficient  between  the  punch  head  and  the  blank 

FRITND:  Friction  coefficient  between  the  die  and  the  blank 

BHFCE:  Blank  holding  force 

Set  0.0  for  punch  stretching  problem 

7.  Read  YVALUE,  PRESTN,  EXPNT,  PRESTS  (4F  10.0) 

8.  Read  TCONTC,  TDIST,  ECONST  (4F  10.0) 

TCONTC:  Criterion  distance  of  the  contact  with  the  die  profile 
To  start  with,  set  this  0.002 

9.  Read  N,  CODE(N),  R(N) , Z(N),  UR(N) , UZ(N),  SLOP(N) , (15,  F5.0,  5F  10.0) 

10.  If  NREAD  = 1,  the  new  input  data  is  to  be  placed  behind  the  nodal 
information  card. 
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1 

? 


PROGRAM  S H*1  FT  ( IN»tlT  .nuTOUT  , TAOf  5*  INPUT  ,T*Pe  f«OUTt>LT  ,PUNf'  ) 
Cr'MMON/GFNC0N/NUMNP.NUMcL  » MED  ( J ? > .OLL.NEO,  NRf?RM,  Y I ELD.  TEST,  ITER. 
1 NO  E AD.NPIINCH.NPP  TnT.CvALUT-.T,  MR  a NT.  , OKIPiO  *P*DIUS*FB!  TNP  , FR I TND  # 
2EC0NST,  FNHE  C • RT  HOT  .OTERAO.TCONTC.TDIST  « *»HFrF 


PROGRAM  BOTH  POP  PUNCH  STRETCHING  WITH  ROUND  PPPC  ILF  AND 
Fpo  OsEO  DRAWING,  PV  j.h.kJm 


COMMON/ M AT F RL /V V AlUF , PR ESTN.E *CNT ,ppp  ctc 
Common/! SOTPv/PVALl 


PROGRAM  SHEET  IS  POP  CONTROLLING  the  DIMENSION  OP  THc  COMPLETE 
OROGRAM.  ITS  OURPCSP  JC  TO  rpry^NT  ASSIGNING  A LACGFR  ^ HAN 
NECESSARY  DIMENSION  POP  ANY  ARRAY  THROUGH  TMF  USE  0*  THE 
POLLOWlNG  ST  ATFMFNT 


COMMON  A ( 500 0 I 


NP  I EL  O IS  THE  D I MS  NS  ION  OF  ARRAY  A.  ITS  VALUE  CAN  PE  DETERMINED 
PRECISELY  PY  PUNNING  THE  crcgcAM  ONCF . 

U*tttM**«*»*******t*:*****t*tt**»****t>**»tK*»*i.>«it»»»»**»*,«it««» 


NF I flO^EOOO 


READ  THE  INPUT  DATA  CONTROL  C *RD« 


TE  S T * 1 . 

RF  ADC  5 , ) 000  ) HED 
READC5.100A)  RV ALUr , T, ACOFF 

RP  ADC  5,  1 0031  ITEP  .NRFaD  , I TC  CNT  , NR  CRM  , NE'JN  CH  , NP  c J n^  , FL  ! M I T 
RFACC5, 10031  NUMNP.wre* 

Rff  A C ( 5 , 100A  ) PNRAC.R  AD  I US. 0 I po AD, R Thc t 
0=40(5.1004)  fc T TNP  ,FR I TND , °HFCC 
RE  AD ( E . 1004)  YVALUE.RRcSTN.E  KPNT ,PRE  STS 
CEADC5.1004)  TCCNTC ,T0 I $T , ECONST 


HEP sOUTCUT  T! TL  = 

RV  ALUE*  VALUE  OF  THF  ANISOTR^P*  DiRAM^T'P 

ACOEFsACCELERAT  ING  DR  CECFlERATInG  CCFFFICI=NT  OF  CONVERGENCE 
NR  F AD  * 0 , IF  TO  PYRASS  THE  dpaDING  STATEMENT  IN  ?U®RCUT*NE  PL  A3T 
I T C ONT  * 0 , IF  CPMOUT  *Tjf  n STARTS  AT  THE  V?=v  «*E  GG  INNING  AND  F [ C ST  / 
SFCCNO  ETFPS  ARE  INCLUDE''  IN  TMp  STE°S  TO  R = COMPUTED 
*1,  OTHER  W T S= 

THIS  INDEX  I c R FL  AT  F D to  T HE  OrTERM  ]*JAT  ION  CF  STE°  SIZE 
NFnRMxNUMRFR  DC  STFPf  ASSIGNED  PE  R RUN 
NPUNCH*|,  IF  CATA  AD?  TP  P = PUNCHED 
=C,  OTHERWISE 

FLIMITsVALUE  OF  (ERROR  NORM » / C S^L UTI ON  NORM)  REQUIRED 
FCR  CCNVERGFNCE 

NOP  I N T a | , IF  NODAL  POINT  DATA  ARE  Tn  FE  PRINTED 
*P,  OTHERWISE 

NlJMNOsNUMF  E°  CF  NODAL  POINTS 
RNRAD*RADIUS  OF  MCW I cowER IC  AL  PUNCH  MEAD 
PAD IUS* RADIUS  OR  THF  FLANK 

Fp  I TNOtF  R I C T I ON  COEFFICIENT  P*TW-EH  BLANK  ANT  PUNCH 
FD ITND=FP ICT TCN  PPTWEFN  “LANK  AND  D T F PROFILE 
01  = CADrRADIUS  CF  DIF  RC0FTLS 
RTMPTrOADIUS  OF  DTF  THF  OAT 

E C P NS  T * S T E P SUE  IN  MAXIMUM  FFRrCTIVf  STRAIN  INCREMENT 
NPE  X*2,  IF  PUNCH  STPe-TCHlNG  NT™  ROUND  ocfc  I L E 
*3,  IF  OF  FP  DRAW  TNG 


YVALUE,  PORSTN,  E K PNT  t PRFSTS  *cr  rr  F XOe  F ? e TH=  wO  RxMAROENTNC* 
CH AO* r TFRl  ST  I CS  OF  THF  "LANK 


SHE  FT 

78 

C 

79 

C 

SHEET 

80 

c 

NFC*M»mbE9  OF  E 0U4T  JON 5 TC-  3E  fOL  V*r 

SHFeT 

81 

c 

NII8EL  *NUfBfC  T*  ELEmEFTS 

S«C?T 

«2 

r 

MO*NO«P*Nr  WIDTH 

SHC«r  T 

e-* 

c 

SHEET 

84 

c 

SHEET 

as 

c 

SH*ft 

8* 

9 VAL1  »0  VALUE 

SHEET 

87 

SHEr  T 

88 

MH AN0«9 

SHEET 

80 

NE08NUMN083 

SHEE  T 

90 

NCNEC 

SHEET 

91 

NEL "NUHEL 

SHp*T 

92 

c 

SME*T 

93 

c 

888* 

SHEET 

94 

c 

CcTERMINE  THE  LOCATION  OF  THE  STARTING  POINTS  OF  niEFgcENT 

SHEET 

9E 

c 

APPAVf  C*N  ARC  * V A 

$mccT 

96 

c 

88*  8 

SHEET 

97 

c 

SHFET 

98 

c 

SHEET 

09 

N 1 ■ 1 

SHEET 

1 00 

K9*M  4-Mi  |MKP 

SHEET 

101 

N3«N2*NUMNP 

SHEET 

1 02 

N4  *N ? ♦NUHNP 

SHPFT 

1 03 

NS»N4*NU®NP 

SHEET 

108 

N6»NS*NU8NO 

EHe’ET 

1 OS 

KT*K6*NU8NO 

SHE^T 

1 06 

N*8N7*NUMEL 

SHEET 

1 07 

N9  *N8  ♦NtIMEL 

SHEET 

1 08 

N 1 0«N9*NUMFL 

SHEET 

1 OO 

N 1 1*N10*NUMEL 

SHEET 

1 1 0 

M2*N1  1 ♦ NU  **  E L 

sheet 

1 1 1 

N 1 3*N12*NUHEL*3 

SHEET 

112 

N1A8N134NUMEL *4 

SHEET 

113 

M 5*N1 4* NU EEL 44 

SHEET 

1 14 

N 1 6 *N 1 SANE  0 

SHEET 

1 IS 

N 17  *N 1 6 *NE  0*8 E AND 

SHEET 

1 16 

N 1 8«N17*NUMEL 

SmE*t 

1 17 

N1 0«N 1 8 ♦NUHNP 

SHEET 

1 18 

N?0*N1 Q8KUVNP 

SHErT 

119 

N2  1 8N204-NUMNP 

<H*fT 

1 20 

K22*N21 ♦NU8NC 

SHEET 

121 

N23*N22*NUHNP 

SHEET 

1 22 

N24.N234.NE0 

SHEET 

1 23 

N25 *N2A ♦NUENP 

SHEET 

124 

N 26  *N  25*NUMN'P 

SH**T 

1 23 

N27«N?6*NUNNP 

SHE*t 

1 26 

N28.N274NUMNO 

SHEET 

12T 

c 

SHEET 

1 28 

CALL  PPELlM(A(Nl)tA(N?l,A(N3),A(NA),A<N5>*A(N6)) 

SHEET 

120 

|F«N?8  .LE.  NEIELO)  60  TO  100 

SHEET 

1 30 

c 

«HECT 

131 

c 

SHEET 

1 3? 

c 

SHEET 

1 33 

•PI TE 16*1 002 1 N28 

SHEET 

1 34 

c 

SHEET 

1 35 

c 

*888 

SHEFT 

1 38 

c 

N28  IS  TOTAL  SPACE  PEOUIPFO.  I®  NEIEL^.  Th*  I N IT  I ALL  V ASSIGNED 

SHEET 

137 

c 

SPACE « IS  NOT  ENO  'GM,  STOP  Tw?  ©PCGPAW. 

SHEET 

1 38 

e 

8888 

Sh*-et 

1 30 

c 

SHEET 

140 

stcc 

SHEET 

1 41 

100 

CONTINUE 

SHEET 

14? 

c 

SHEET 

1 43 

CALL  PL  AST  (A(Kl),AfN?)vA<N3»'ArNA1.AfN5)tA<N*l(A(N7),A(N8),A(N9>« 

SHEET 

144 

1 4 ( N 1 0 ) * A f N 1 1 1 • A ( N 1 2 1 t A f N 1 3 ) « A ( N 1 A ) « A f N1 51 t A ( N1 6 ) * A ( N1 7 )«A(N18)« 

SHEET 

143 

?A(N191vAfN201«AfN?11,A(N?2),A(N?7l,A(N?AlvA(N2S) ,A( N26) #4<N27> 

SHEET 

146 

3,  N0,NEL.FLIM|T.TTC0NT, ACOEP.NrF* 1 

SHEET 

147 

r 

SHEET 

148 

1000 

9 •* AT  M 2 A6  > 

SHEET 

1 40 

1001 

rnOMATf///«  THE  DIMENSION  OF  THE  ABPAV  (41  IS  TOO  ?N  ALL  •/ 

SHEET 

150 

1®  THE  S17E  OF  TH?  AOOAV  (A)  MUST  PE  *,  17) 

SHEET 

1 31 

1 007 

EOP¥AT(//8  THF  NFCCESS apt  S I 75  OF  THE  AOCAV  (A)  ISA.  17) 

SHEET 

152 

1003 

E09MAT( 6T5.F10. 0) 

SHEET 

IS' 

1004 

fpcvat (3F1 0.0  ) 

SHE*T 

154 

1 005 

rnPMATf  4 I S *F l 0. 01 

SHEET 

1 S3 

C 

SHEET 

1 56 

STOP 

SMEeT 

157 

SNH 

SHEpT 

1 SO 

SU9P0UTINE  POEL  T m ( R # , UP • U7 • CCPE • SLOP ) 

SH«ET 

160 

COMMON/GENCON/NUMHP,NUMFL*m?D( 1 21 .OL L t NE 0 • Nr OPM , v IELD , TE ST • I T E© • 

SHEET 

1 61 

1 NO® AP  *N PUNCH  *KPB INT*  PV ALUE,T  f M«  AND, PNPAD, P AOIUS. FP ITN» ,eB I TNO • 

120 


s*-c‘*r'r  16?  ? -CONST  . o^H^n  , OTHC  T f F«AD  , TCP^  ♦r  , rnj  , °MFre 

«MfCT  1 c 

SH*7- ST  1 64  <“*•«**%******•  *M«*(*tt***t«*tti<l  *****  ***********  *****  *******  IN  0 

snrrr  165  r PEAT  AN  O PRINT  PP  CONTROL  I NS  nc  *i  TJ  ON  *Nr  MATfc  P f AL  COCCcorits  «A'N00l3 

she  FT  I 66  C****  ************  ****’*********■**'*******•**  A*  *******  *****•**»**>**•»•  *****1N0014 


SHfPT  1 FT  C 
SHC-^T  1 f-H  C 

CMfST  l*q  OIMSNSION  R C l 1 . Z U ) . CCO*  < 1 I . V>9  < I V ,U7  < 1 > . SLOP  I 1 ) 


SNF^ 

170 

C 

SHcr  T 

1 71 

50  CONTINUE 

?mscT 

1 7? 

WO  I TF  (6,20001  HEO, NUW^ C #NUMFL 

«H?Pf 

1 73 

CALL  HAfiniO .,V IFLC» 

SHCCT 

I7A 

»°1  TE ( 6 ,201 01  VIELT 

SHF^T 

175 

WRITFC6, ?01 1 1 

«H*"P  T 

1 76 

WRITF(6,2012»F4ril»S,PNC4O,P!ECAr.pTHFT#EDyT^D#FpyTNr 

SHE  -T 

177 

w®  1 tc  | 6 , 201  1)  WALL"  ,ocectn,ek  c*  7 , dcejts 

SHS'T 

1 78 

WO  I TF (6, 20 14  1FC0NST 

SMC  e-  r 

t 7<J 

• 9!  TE <n  . 10091  I TF  R 

SHecT 

I 80 

C 

SHSr  T 

1 “1 

• ** 

Smo-t 

1 8? 

C 9E4T  AND  PBIn’  C*  NP^AL  CT  1 NT  pjiTj 

MA  IN0031 

SHF  ST 

I 8.1 

***  wa  T N0032 

SH^fT 

1 «4 

r 

SMeer 

1 95 

C 

SHE  - T 

1 «6 

L*  0 

MA  TN001A 

SMPfT 

1 87 

(P  ( NPPINT.EQ.C1  GO  T*}  60 

Smpft 

1 88 

WC  ITE  (6.1114) 

SHFPT 

1 89 

•PIT'  (6.2004) 

MA I N CC  32 

SHP^ 

1 90 

60  PE*n  (5,1002)  K.COPF(N)  ,P(  N)  ,7(N)  .UMN)  ,U7(N)  .SLCOfN) 

5HPC  T 

1 <51 

c 

SHPP^ 

1 92 

NL  *L*  1 

MA  IN  00  36 

5 HF  CT 

1 9* 

ZK-N-L 

MA IN0037 

«MFF  ■*■ 

1 94 

I F ( L .EG.  0)  CC  TC  to 

1 95 

OPs(R(N  )-P(L)  )/?* 

MA !N003e 

flMTp  T 

1 96 

07* ( 7(N)-?(L ) )/Z* 

MA IN003C 

SHFCT 

197 

DUP*(  UO(N)-UR  (L  ) ) /7Ht 

SHEET 

1 99 

OU7« (UT(N )-U? (L  1 ) /7  * 

SHE  c T 

] 99 

OS*(SLOF(M-SlOF(L>  )/7x 

Shsct 

200 

70  L*L ♦ 1 

MA  TN004 1 

SHE  F T 

201 

C 

<MTeT 

202 

! F ( N— L ) 100.9C.80 

M A 1N00A? 

SHF  FT 

20? 

C 

SHF  S T 

2 04 

80  C0OF( L»  *0.0 

MA  7 NO  04  * 

SHF  FT 

2 OS 

OIL  )*0<L-1>4D» 

M A TN0044 

SHC  FT 

2 06 

SLOCfL) *SL0P(L- 1 )40S 

SMFFT 

20T 

?(L)*7(L-1 )*o? 

SH*FT 

208 

UO(L )*UP (L- 1 > 4PUP 

SHFFT 

20© 

U7<  L > *U7( L-l  ) *PU  7 

SHF  FT 

210 

GO  TP  TO 

M A TN0049 

*mf  f T 

211 

c 

SHPfT 

212 

90  T E f Mjwnc-n ) 1 00,1  1 0 » 60 

SHFFT 

’ll 

100  WOTTE  (6,2009)  N 

MA  I NO  0 5? 

SHF  ST 

214 

C 4L  L EXIT 

MA  T NOOE  2 

SHP-CT 

2 1 S 

110  CONTINUE 

MA IN0054 

<HFFT 

216 

C 

SHFFT 

217 

*°r  TE  (6,2002)  ( K ,rnr>E  (*  ) ,P  (F  ) . Tf  E ) ,uc  ( * ) ,U7(E  ) . SLHP(*c  ) .«  * 1 . 

NUMNP  ) 

Shppt 

218 

c 

*mec  T 

21<5 

r 

SHFFT 

220 

NFO*  34NUMNP 

SHFFT 

221 

we  ITF (6, 1 ! 22  ) NPQ.MFANC 

SHS*T 

222 

C 

SHEET 

223 

• «i 

SME  F t 

224 

100’  FOCMATI 15  .F5.0.5E1 o .0  ) 

SHEET 

2 25 

1001  FORMA’! 1615) 

SHFFT 

226 

1004  EOCMATI  18,211  1, 2F10.6) 

SHFF- 

227 

1005  FOPWATI2T5.4E10.01 

SHFFT 

2’8 

100*-  FOOMATt//  * THF  NOHAt  POINTS  AT  WHICH  EQPCF  CALCULATIONS  ABE 

DESIR 

SHF  E T 

229 

1F0*  //  2015) 

SHPFT 

210 

1007  FOPMA  T1  |M  1 # 1 5 x t JCH  LINEACtr  CIETeiptfTFC  "OUNOABV  STPFSSFS/ 

SHFFT 

211 

1 / 60H  NODE  I NODE  J PPESM  Bf  I PFF5«;U^g  j SHE  AP 

T 

SHEET 

212 

2 14H  ShFAp  j| 

Sheet 

233 

1008  FOPMATI 21 9, 4E 1* .« ) 

SHFFT 

214 

1009  EOCWATI///*  MAXIMUM  NUMPEfc  QF  ITERATIONS  ALLOWED  F0»  EACm  INCREMF* 
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UP! T 1 *UP ( I 14EST  »B 
U7!  I 1*UZ! I I *E  STAR 
CONTINUE 

CONTINUE 
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1 r ! I .LE.  ND  I F 1 rCPFf  11*2.0 
CONTINUE 

COOE!NUMNP  1*3.0 
C 00  E ! 1 1*7.0 

IF ! NOE  * .EO.  21  CTDE!  11*3.0 


NTPMM 1 *K TOUCH- 1 


«Hf  fT 

4!  * 

C 

4 14 

C 

e Hf*’ 

41  5 

C\ 

SH'«T 

41  * 

r 

C0*CVT£  THE  V I FL  r ST  PE  S S AND  THC  40RK  MA  9C  c N * NC.  5 A Tf 

417 

c« 

S He*» 

4|#* 

c 

«H*C  T 

41  9 

DC  220  ►.*!  • MJH6L 

SMf 

420 

CALL  H4P0CTFPSI  4,M  tWlM» 

SH^eT 

4 ? 1 

C4LL  *4QC2  (Tf  PS  14.N1.  vxf  N 1 1 

?MC  " T 

42? 

220 

CPNT1 NUr 

Sm**'* 

4-»T 

c 

«M5«  T 

4 9 4 

r 

SHFffT 

4 2« 

C 

«hfc  t 

426 

r 

SH^*T 

427 

r. 

SM**T 

4?8 

•RITztS.  10071  NSTFP 

«HFrT 

429 

r 

SH«;:T 

430 

COMVt l l 1 l 1 1 • 

S«?fT 

4 3t 

r 

SHFFT 

432 

650 

C3N'T!NUc 

SH?  FT 

4 33 

•PIT*! 6. 1031) 

SH'PT 

4 34 

C 

SHF^T 

4 35 

r t 

4***i 

SHF  FT 

436 

c 

DFT4IL  Oc  THE  PRESENT  C ONS I GUR 4 r T On 

S«FFT 

4 37 

c 

cohIcSTKE  os  ANCLc  Chi 

ShFFt 

4 36 

c 

C°H!« Cn  F I NF  PS  ANGL*  PHJ 

SHF  * T 

4 39 

c 

OL*FLEMFNT  length 

5HCFT 

440 

c> 

ShFFT 

441 

c 

SHF  F T 

442 

OG  6Q0  N«1  , NL'MEl 

SHFFT 

443 

N°  1 *N ♦ 1 

CHCCT 

444 

OPrC(M-CfNPl  > 

SHFFT 

445 

07*71  NO I1-71N1 

SHFFT 

446 

DL  INl.SORTi  DR  *r>F  407*9  7 ) 

SHF  F T 

447 

S°Hl ( Nl *0  R/DL  IN) 

ShFFT 

448 

CPMIIN) «OZ/DLIN) 

«hFF7 

449 

PH! (NlsASlN(SPHlfN) »*180./3.1A1*« 

SHfc  T 

450 

• 9 1 T*(  6.10301K.  PHI  1N|  .THIC*  CM  .OL  «M 

SHF®"T 

451 

690 

cont  inlie 

SHFFT 

45? 

c 

5Hr  FT 

453 

po  1 «o VALUE ♦ 1. 

SHFFT 

454 

RCCNST*3*4PP1/ ( 2.*C  1. 4RVALUE  ♦PVALOF  ) > 

ShFFt 

455 

01  l .1 ) *CO|*RCrNST 

ShF*t 

456 

Oil, 2 | *R  V AL  UE  *R  CON  S T 

SHFF  T 

457 

0(2,  t t *c< 1 ,2 » 

snsr 

458 

Ol 2 * 2 ) *0 ( 1,11 

SH^FT 

459 

CLAHOA*  PPJ*RP1/  f 1 .UPVALUE 4F VALUE  » 

SHFFT 

480 

POI S0N*RVALL£/R9l 

SH*«?T 

461 

c 

«HFCT 

462 

5*0 

S WFCT 

463 

2001 

CONTINUE 

SHEFT 

464 

<«KM 

SHF  * T 

46« 

c 

SHFFT 

466 

c 

SHFFT 

467 

CALL  ST  IPS  (fi,  7. UR.U7.CCCE. SLOP, vv ,r X,S°HI ,CPh| ,OL  *f c s. 

SHFFT 

465 

1 

1 THICK  .ALPHA .GAMMA , ETA .FRNSCF .*F  , A .B.NOl 

SHFFT 

469 

r 

SHE  FT 

470 

Ci 

»***< 

SHFFT 

♦71 

C 

INTRODUCTION  OF  POUNOAPV  CONDITION 

SHFFT 

4 7? 

c« 

SHFFT 

473 

c 

SHF  rT 

4 74 

c 

SHFFT 

475 

CALL  MCPIFY  fCCDE  , A.  ",  AL  PHA,  GAMMA,  FT  A » NUMNP • N’FO.mcano,  ec^FCE  1 

SHFFT 

476 

c 

SHFFT 

477 

c 

°ANCFC  SYMMETRIC  SOLUTION 

SHE*T 

478 

r 

SHFFT 

479 

CALL  TR  IA1NFO.M8ANO.A 1 

SHFFT 

480 

CALL  "ACKSINEC.KPAKr.A.Pl 

SMFE^ 

♦81 

c 

shfft 

♦ 8? 

c 

SHFFT 

♦ 8? 

C' 

4*4*1 

SHF  FT 

484 

c 

3COTl)P«AT  ION  OS  UP  is  COMPUTED  FPfM  f F R T UPB  * T T ON  pr  u 7 *0*1  NODES 

SHFFT 

4 85 

c« 

4***1 

SHEFT 

486 

r 

ewrrr 

487 

c 

SHFFT 

48  8 

C 

SHf  *T 

4 89 

C 

SM«ET 

490 

DC  101  N * 1 , NU*N  C 

SHEFT 

491 

IP  « ”»*n-  2 

Smc*T 

4 92 

T 7 * IR  ♦ 1 

SHFFT 

493 

ls|r0r>F!M  .EC.  4.01  "4  1R)*  — 1 I 7 1 6 AL  FMA 1 M *G*  w A 1 N ) 

SHF»T 

494 

101 

CONT  I Nil? 

CM*  FT 

495 

C' 

4***i 

SHEET 

496 

c 

TO  05  T a | N AN  ESSTCIFKT  CONVERGENCE  A CPeS  IS  CCMPuTSp 

SHFFT 

4 97 

c 

TH«  ACCELERATING  COFSSICIENT  |*  OF  TF  R M | NE  P IN  such  a »Ay 

124 


1 


?H??T 

«Qft 

C 

of  3TIICRATION4L  TFF*  TIMES  4 CCEL  AT  I MG  C^S^lCI'KT  IS  NEVE° 

SHF^T 

494 

c 

Gp c»Tfo  than  the  IM*“IAL  vault  . si  t pe  a *oaCT!Pn  cf  IT. 

CHC  t j 

son 

c 

E.G.  A-»*2  VflAS  HALF 

SHEET 

“■01 

c* 

1*  •*  1 

SHEET 

SO? 

41*10. 

5 N*"  *"  T 

5C3 

5t*reT 

50* 

c.ONrr***r.o 

*Mt  r t 

SOS 

Or  103  1*1 ,KUMKC 

SHfPT 

506 

IP<A6S<UO<IM  ,LT.  .COOOOOI  UjCl M *0 . 

SM-fT 

5 07 

I* C ARS CU7C I M .LT.  .0000001  1U7M  1*0. 

SHCFT 

see 

UU4 ( I \ *U« C 1 > 

SKCT 

509 

UU’( 1 >*u?< 1 » 

SSfFT 

St  0 

I ?*3*T-t 

Shcct 

5 1 1 

10*17—1 

SH'TfT 

Si? 

1P|»<10>  .EO.  0.10°  TC  102 

SHF  FT 

5P 

CTT!«A*SIU»m/A(IP»l 

SMC  »T 

514 

102 

IP|R(r?l  .£0 • 0.1 GO  TC  1 03 

SHce  r 

515 

ccf  ">x  aps  (uzni/prim 

5hF“t 

si  e 

A | * AM  IN  1 ( A 1 , C 0* 1 ,COF?| 

f He  FT 

517 

1 03 

CONTINUE 

cmci-t 

51  E 

c 

SHcr 

* |P 

r 

SM*fT 

5?0 

1 05 

CONTINUE 

SnEEt 

521 

ir(roNcnF  .EO.  |.|A?*A?*5. 

S«FFT 

5 7? 

!c<  M .EO.  1 .0  .AND.  A2  .GT.  10.1C'3P1*5. 

fHcrr 

52? 

I F < COE  1 .EO.  0.1CCFI*?. 

5 MFrT 

5?4 

A l*Cr»F  I / A 2 

SH^CJ 

5?5 

I F < A l .UT  . .OOOOMGC  T.O  ?300 

f HF  z T 

525 

l E < A 1 .GT.  | .Cl  A 1 * 1 .0 

She  ST 

52-» 

ACC***  A 1 

«MF  E T 

*?e 

• «ITEf6#l04 1 a 1 

shee^ 

5 24 

c 

SHEET 

530 

c 

«Heer 

531 

C * 

r**«l 

SHFFT 

532 

c 

OPTAIN  NE»  VALUE 

SHEET 

533 

c* 

SHEET 

*34 

c 

SHEET 

535 

Co  J33  1*1,  K-UMNP 

5HEEr 

*36 

I 2* 3*1- I 

SHEET 

53T 

10*17-1 

SHE  E t 

5 78 

I S*  I ?♦! 

SHEET 

5 34 

UP  I I 1 wUUOf  I 1 

SHE  r ▼ 

540 

U7f  I I-UU7C  t » 

5H*e  T 

541 

UPC  II  *UF<  IMPCIP14AC3EF 

SHE^T 

542 

U?m*UlU  !4e<  \7\4AC<3EE 

SHE  c ▼ 

54  3 

SLOFC I MSLOPC 11  AFC  IS14AC0EF 

SHEET 

*44 

130 

CONTINUE 

5MEeT 

545 

c 

SHEET 

546 

• °1  TE  (6.10161  K 

SHE^T 

547 

• OITEC6,  1006)  K 

5H*  c T 

54  6 

c 

SNC*T 

544 

c 

SHE  ET 

E50 

c« 

i***< 

SHEE  T 

551 

r 

COMCUTC  KCSM  CF  ftOCC  AND  NCO*'  OF  SOLUTION. 

SHEET 

552 

CA 

i»**i 

SHEET 

553 

C 

SHEET 

554 

*NC=M  * 0. 

SHE  *T 

556 

SNOOM  * 0. 

SHEET 

5S* 

DO  134  1*1,  KUMNC 

SHEET 

*57 

T 7 * 3*  T - 1 

c HE  ET 

SSF 

10*17-1 

Smfet 

**9 

I S*I 74 1 

SHEET 

56  0 

ENOPN  * FNORM  4 PCIPIARCIP)  ♦ o<I7»4pC  IZ1  ♦ F<IS»**»<1S> 

sheet 

561 

sncom  * snorv  ♦ ut(  n*unn  ♦ u?cti*u?cti  4 elopc  1 i^sldpc  i i 

SHEET 

56? 

134 

continue 

SHEET 

563 

c 

SHEET 

564 

c 

SHEET 

565 

ENOOH  * SOPTCFNOPMf 

SM*ET 

566 

SNCOM  * SOOTCSNCPVI 

SHEET 

567 

F SN OP*#ENOON/«NnPM 

SHE  E T 

568 

• 0 TTEC 6 • 10 1 5 » SNCPM, ENCOM.fENOCM 

SH^ET 

564 

OO  T?6  1*1.  NUMKE 

SHEET 

570 

I 7*  ?*T-  | 

SHEET 

571 

TP* I 7-1 

SHE  eT 

57? 

I S* I 7 4 1 

SHEET 

573 

»»!T'i6.ioo?i  !.«(  j»i.  ec  1 ?».(■<  is»,o9i« » .u7«  r • • siopi  i i .*<  1 1 ,rt  1 1 

SHEET 

3T4 

776 

C3NTINU* 

SH*ET 

575 

c 

SHEET 

578 

1 31 

FONT INUE 

SHEET 

377 

c 

SHEET 

5T0 

c 

SHE  E T 

374 

c 

SHE'T 

5 50 

c* 

i***i 

SHe  E ▼ 

561 

c 

COMPUTE  STRAIN  FPT.M  Tnr  Nc*  GUESS. 

S^ET 

S82 

c 

rose  1 ,n>*incpenent  cf  *esioian  strain 

125 


$ME?T 

5*3 

C 

F°SC2*N>*  INCREMENT  Or  TANGENTIAL  STC  A I N 

SMFC» 

5-4 

r 

TPSC3 ,N) *INCPEMSNT  CE  THICKNESS  STPAfN 

Sm^ct 

S«5 

C« 

SMte  t 

EES 

c 

S^ET 

EE? 

c 

SHEET 

SE% 

c 

SHEET 

5*0 

no  *00  N«1  • NUEE  L 

SHEET 

500 

NO  |*N4 1 

SHEET 

501 

OLL*nLf N 1 

SHEET 

50? 

cpw*SOHI  < N > 

SHEET 

503 

C»H*CPH1(N 1 

SHEET 

«04 

4U«UO fN »4uaf NCI  1 

SHEET 

*05 

AP»P(N) ♦ O | NO  1 ) 

SHEET 

506 

CC«C(N»-C(NP1 ) 

SHETT 

5 ST 

D»«?(M>1  )-?<N) 

SHEET 

50ft 

DU* UP f N »-UB  f NO  1 > 

•MEET 

500 

Oft*L7  < NCI  1-UZ IK > 

SHEET 

600 

E*1  *1  •♦2.*D0*DU/DLL/r>Ll4  2«*0  7*C'.*/0LL/’0LL4  (mj*DU*0**CK  )-'DLL/CLL 

SHEET 

601 

c 

SHEET 

6 0? 

rPSI 1 ,N)  «ALOG  <F*1  * /2. 

sheet 

603 

EOS  ( 2.N  )=ALOGf  1 . HU/AP  ) 

SHEET 

604 

TPS  (3  ,NI  «-6PS  <1  ,N  )-F°S  I2.N  » 

SHEeT 

6 05 

ftOO 

CONTINUE 

SHEET 

6 06 

TEST  *0  ,0 

SHEET 

6 0T 

c 

SH«-ET 

60ft 

c 

SHE  v T 

6 00 

c 

SHEET 

6 10 

c« 

SHEET 

61  1 

c 

COHOUTE  INCP5*ENT  PE  EEPfCTlVE  STRAIN 

SHEE  T 

61? 

c< 

1***1 

SHEET 

6J3 

c 

SM«r  ET 

614 

«cr TE (6 , 1 026 > NSTEP 

SHEET 

615 

OP  2 22  N«1 , NLMEl 

SHEET 

616 

NO 1 «N ♦ 1 

SHEET 

6 1 T 

TEICOOECNI  .EO,  3.0  .ANC.  COCECNDj)  .EO.  3.01G3  TO  ?22 

S H?CT 

6 1 ft 

ES*tr>S(  l ,N> 

sheet 

610 

ET  * EPS  f 2 , N ) 

SH«ET 

6 ?0 

QHABzCo  1*  (E  S*ES*ET*ET  ) ♦ 2 • *P V A LUF* ES* E T 

SHEET 

6?  1 

c 

SMPCT 

ST? 

EOSf4.Nl *S0CT (?.*PCCNST*CpAC/3. 1 

SHEET 

623 

c 

CM»  T 

624 

c 

SH^ET 

625 

STSf  1 .N1  «Cl_4MD4*f  ES  ♦ Op  1 S^NEE  T 1 * YV  ( N l ✓ fTS  < 4 , N ) 

SHEET 

6 26 

STS! 2 ,N )*CL A*CA4f ET  4 On  I S0N*E  sl*W(M  /E°S  ( 4 ,N  | 

SHEET 

627 

c 

SH*r 

628 

c< 

»***! 

SHEET 

620 

c 

COMPUTE  STOESS  CISTCIPUTlON 

SHEET 

630 

c 

STSf 1 ,N) «MEPirt 4N  STOESS 

Sheet 

631 

c 

5TSI 2.N ) *C IBCUMFECFNT I4L  STPESS 

SHEET 

632 

r 

STS <3  ,N) *EEE*CT  IVF  STRESS 

SHEET 

633 

C< 

»***] 

SHEET 

6T4 

c 

SHEET 

6 35 

lEfESNCC*  .LT.  ELIOTT )TEST*1 .0 

SHEET 

636 

eS*?TSf  1 ,N  » 

SHEE  t 

6 37 

ET*«TSf 2 ,N> 

SHpP-r 

53ft 

?c  STPS*ES*ES*E  T*E  T—  ? , » c C I S0N*B  S*  r T 

SHE  eT 

630 

STSC3.N1*  SOOT ( EBSTPS 1 

SMP'T 

640 

c 

SHEET 

641 

»«MTE<  6,  100  3 )N,  <EOM  T ,N  ) , I * 1 ,4  ) 

SH^ET 

642 

??? 

CONTINUE 

Sheet 

643 

• p I TE ( 6,1027) 

SH*ET 

644 

pn  4 TO  N*  1 , NU  Mt  L 

SHEET 

645 

4 T 0 

•PI  TE  (6  ,2251  IN,  I STSf  ! , N)  , !*1  ,3) 

SHEET 

645 

c 

SHEET 

647 

c* 

r***< 

SHEET 

64ft 

c 

CHECK  WHETHER  4CCEC  IS  TCC  LARGE  t0  CAUSE  A OhvSICALLY 

SHEET 

64  0 

c 

UNACCEOTABLE  SOLUTION.  WWENEVEO  COMCUTEO  H^crriAN  STPFS*  pfco*es 

SMErT 

650 

c 

NCGATIVF  ADJUST  A CTEC  VALUE) 

Sheet 

5ft  1 

c« 

>***! 

ShE  * t 

652 

c 

SHEET 

653 

CTNCPE  *C.O 

SHP*T 

654 

STS  1 *0  .0 

SHEET 

655 

r>c  431  n*1  ,NU*El 

SHEET 

656 

S°  1 *N ♦ l 

SHE  P T 

657 

lEfCOPEfN)  .EO.  ‘*.0  . 4NC  • COOFCNP1)  .EO,  3 • 0 ) S T S f 1 ,N  1 * 0 • 

S^ET 

658 

431 

1 B f STSf  1 ,N ) • L T , STS1  )CCNC0C*1 .0 

SmE  eT 

A«0 

r 

S*E  * x 

660 

c* 

r***i 

SHEET 

661 

c 

CHECK  WHETHEO  fr^rcc  NORM) /( SOL  UT ! ON  NOP*)  1?  LcSS  THAN  ELINIT 

SHPP  T 

6 5? 

c 

1B  YES,  THE  SOLUTION  !«  c INAL 

SMPe* 

66  3 

c< 

SHEET 

664 

c 

SHE  P T 

66« 

c 

sheet 

666 

c 

?H*r  r 

*67 

JEfESNOCM  .LT , ELfMlTICO  TC  4 ’P 

126 


r •“  i 


c.  ►**■*■▼ 

‘fP 

C 

SMe5T 

66C 

fF(«  .Gf  . ITRC  If,"  ▼ r 2307 

cue  e T 

*7C 

r. 

SHFfT 

*7i 

c 

SHF  c t 

F?? 

or  1900  *J«t,N(IMEL 

smt  r ▼ 

673 

T^(*PSI4.M  .fT.  .0OO001IG'*  TC  1 Of*0 

C wc-  er 

* »4 

NO | *N ♦ 1 

Smc  c y 

<T5 

OU*  *UC  1 N ) *U«  c ► O i > 

5h*  Ft 

*T* 

C01?(i° | I * J. 0 

c Mc  r T 

6 77 

COC?fNI*3.0 

*76 

1000 

COSTJ  NUr 

tHfer 

ST<; 

c 

CMC  cy 

FRO 

r,r  tc  TOO! 

Swrcr 

661 

200C 

COM  T I NUF 

CMCCy 

* R? 

2200 

1NUS 

cmtct 

F RT 

C 

SMPCy 

* *»4 

4 IP. 

CONT 1 NUF 

«MC  c y 

6«« 

l*fFS6.rc*  .CT.  PLJHITlC-n  777 

«HS'T 

f“* 

e 

Sm^ct 

R*7 

c 

SH*ccy 

* RP 

»CI  TF  !6 ,2600  1 

SHCCT 

6 RQ 

r 

SMCt'  ▼ 

^on 

c 

E * A*  T M*T  ION  CN  PCUNCACy  ASSl»VC*7  TONS 

CMCCy 

*91 

r •»** 

SMCcy 

* ^2 

r 

CMC  c y 

PCI 

DC  ?50  1 >1  , nChe  f n 

<;hcct 

4A 

NsNTOUCh- 1 

C^  C y 

6 95 

(6  «L  H • 1 ICC  TC  2*0 

$M=cy 

C>9* 

TflUCHf  ! )=(  ;(M  )*U7(  Mppn  = A['-pnh£  r »*(  ? IN  > *u  ? 1 * »*PN;Ar-o\HFf'»4(P(N  » 

<;mcct 

6QT 

1 4UR  IN  1 \ * { R < N 1 4UO  | N 1 » 

SHCC  y 

PS® 

TCUChI  I IsSCfiT  (TCUCH(  1)1  —PNC  AT 

5Hccr 

6 99 

•F I T“ f 6.29001N. TOUCH!  I ) 

CH«T 

TOC 

250 

CONTI NUF 

CMJF  T 

TCI 

WOT  TC  *6 , 2651 

5MFFT 

-o? 

C 

Cmcct 

701 

c 

SMfCT 

T04 

C •*** 

CMC  py 

705 

c 

CP^CUTF  T Pc  CIST  ANC  eS  Pr  THE  FFFr  NT  C E FCQM  TH£  OIF  OPPF  J \_Z 

FiHCCT 

76f 

r **  ** 

ChfCT 

70T 

c 

SMCcy 

7 0® 

cr*IE*CT»-3T*CIECAC 

SHEET 

T 09 

00  2*0  !*1,NCHFCK 

CMC  cy 

71  0 

N * N C 0 NT  C ♦ I 

SHCC-r 

71  1 

CONTACT  I ) *(CI  EC6C-7  ( N 1-U7CN  ) )*  CC.TER40-7CP  1-UM  N 1 !♦(  PDTF-e  IN  1-IIP  IN  1 

CHEST 

Tl? 

1 )♦  ( OOIS-O  (N  1-UC  IN  1 1 

CMCCT 

7|  3 

C0NT4riI>*SCRTICCNTAC(I»)-0ISC4r 

Sweet 

t 14 

760 

ER I TE ( 6.2701N.CDN7AC I I ) 

?MceT 

715 

c 

CMCCT 

71  6 

IF ( A5SI  TOUCH!  1 1 1 ,LT.  .0001  1 TCuCH ( j ) *0 . 

SHE  FT 

717 

IF  ( IRS! CONTAC ( 1 1 1 .IT.  .0001  1CCNTAC 1 1 » *0. 

CMC  c t 

71  fi 

c 

SHCcy 

7J9 

c 6**4’ 

ewrer 

7?  0 

c 

CHECK  *1 N ROUNCARV  OVER  PUNCH  *-fAF> 

SHCcy 

7 ? 1 

r*t*«  i 

SHEET 

7?? 

c 

SHE  c T 

723 

IFITCUCHIll  .CE.  0.1CC  TO  3000 

SHEET 

7?4 

c 

SHEET 

725 

• r |T  = |6 . 3100 1 

CMC  e y 

7?* 

c 

S HE  cy 

7?7 

r ****, 

SHE  c T 

726 

r 

NCOS  AT  NCHW1  is  I NS  I T F PUNCH.  C0V=>UTF  AGAIN 

SHccy 

7 29 

C4*6*> 

SHEET 

**30 

c 

SHEET 

▼ 31 

TCHCOF  *0. 0 

SHEFT 

732 

Tf NSDE*0.0 

SHFET 

773 

GO  TC  2101 

SHCct 

7 34 

c 

SHccy 

735 

c 

SHEET 

736 

3000 

CONTINUE 

SHEET 

717 

TCH'-nti  i . o 

?mcct 

736 

r 

7 30 

r 6*»*i 

eMfCT 

740 

C 

CHECK  ON  ROUNCARV  OfF  OP£)F1Lc 

cmFPT 

74| 

r ****< 

SHFET 

742 

C 

SHE  *"  T 

743 

I F | FONT  AC  ( 1 | .Gf  • 0.1CC  Tn  3500 

SHccy 

744 

•PI  TF  I 6 . 1400 

SMFFT 

745 

r 

SHE  c y 

746 

SMFCT 

747 

C 

NOOF  A T NCONTC*  1 IS  !Nf  T^f  0|e,  CC*»PUTF  AGAIN 

cMcr  y 

74  6 

r***6< 

SHCPT 

740 

r 

CM«CT 

750 

NCONTC  *NCONTC ♦ 1 

cnee  r 

7«| 

NCI *NC0NTC-1 

SHE*"T 

752 

co MFC* I NCONTC  )*FONFCFf  KC , | 

127 


5 HE  5 T 

753 

C 

SHEET 

754 

GO  TO  2101 

SH*-*T 

755 

3500 

CONTINUE 

SHE  FT 

75* 

IEINPEX  ,EQ.  21GC  TC  3501 

SHF  FT 

**57 

IP f CR (NOIE )*UP4 NP1F » » .GT.  OOIF1GC  TC  350! 

SHF  FT 

75  8 

K0IF*NC1E-1 

SHF  FT 

75<J 

*9! TC < 6 ,34501 

SHF^T 

750 

c 

SHF  FT 

761 

ri 

***** 

SH^FT 

Tfc? 

r 

N*>  I F NDDF  H4S  Ppc  H POOUGHT  !MC  CONTACT  *ITh  DIE  PROFILE 

SHEET 

763 

c 

*•**> 

Shf  r t 

76* 

c 

SHEET 

765 

pbnECF  (NCI E >*C.0001 

ShFFT 

766 

GO  TO  2101 

SHEET 

79** 

c 

SHF  c T 

769 

7501 

CONTJ  NUE 

SHF?T 

76Q 

T1NS0C= 1 .0 

S*F?T 

770 

TEST* 1 .0 

SHEET 

771 

r 

SHEET 

7T? 

c 

♦ ***i 

SHEET 

773 

c 

COMMUTATION  CF  ECICTTCK  CCE  Ep  | C ISN'T 

SHEET 

774 

c 

4*4*1 

«hc  FT 

775 

c 

SHEET 

776 

MUDEX*0 

SHEET 

777 

pe»N»-EO*PNHEC 

SHEET 

778 

40  I TE  (6,233  MFPNECE«N»,N*I,NUM1  ) 

SHEET 

779 

401 TE  f 6 ,23  1) 

SHFFT 

780 

c 

SHEET 

781 

DO  ? 30  1*1, NU* 1 

SHEET 

79? 

17*  **  I-  1 

SHEET 

7e3 

IP*  7-1 

SHEET 

784 

IP ( CODE  Cl!  ,E0.  3,01 GC  TO  230 

SHEET 

795 

T F ( I .GE.  NTOUCH1CO  TO  290 

SHEET 

7 E 6 

IP(I  .LE.  NCCKTC  .AND.  I . GT  . NTH  ICO  T 2«  1 

SHEET 

797 

I P ( I 4.T,  NTOLCH  .AND,  I .GT.  NCCNTC1G0  TO  230 

SHEET 

799 

IFCKOFX  .EO.  2 1 GO  TC  230 

SHEET 

780 

Ie  ( I .LE.  NDIF1GO  TC  282 

SMc  ET 

790 

c 

SHFeT 

791 

280 

C0KTINUc 

SHEET 

792 

DU**  1*17(1)  *U7  <11  ♦PNPAr-OPNHFDl/PNF  AO 

SHEET 

793 

OUM2* (R  C I )4UC  ( ! > ) /PNC  AC 

ShEct 

794 

POJ TN rpQITNP 

SHEET 

795 

T p ( PP  IT N .EO.  O.IGC  TO  230 

SHEET 

796 

GO  TO  283 

SHEET 

797 

c 

SHC-ET 

799 

291 

OUMl xDIECAC-Z ( I )— U 7 ( T ) 

SM^ET 

799 

ouM  2* on I E— R ( I >-LO( 1 ) 

SHEET 

900 

pP!TM  = po IT  N 0 

SHF  P T 

901 

lp(PO]TK  .EC.  O.IGT  TC  230 

SHEET 

90? 

GO  TO  2e2 

SHE  c T 

903 

c 

SHEET 

9 04 

28? 

DU* 1*1. 

SMCFT 

905 

0UV?*0. 

SHEET 

906 

PR  I TN*p  C TTND 

5Hp  ET 

907 

Ic ( PR  I TN  .EO.  O.IGG  TC  230 

CHppT 

909 

r 

SMEE^ 

909 

?«3 

CONTINUE 

SHEET 

9 10 

PN*PP ( I 7 )*CUM 1*PP(  IP  >*DUM2 

SHF  c T 

91  1 

OT*pp  C I 7 >*DUM2-FP ( IP  )*0UM1 

5He  ET 

912 

IPfON  , EO . 0 . ) GO  TP  230 

SH^FT 

913 

XMUsCT/ CN 

«HPpT 

0 14 

4®T  TF ( 6 ,232) I • xmu 

SHEET 

915 

XM(J*  XMIJ/ER  I TN 

SHEET 

«1  6 

TP(XMU  .GT,  1.0?  .rc.  x *U  ,LT.  .PP)MUDFX*1 

SHEET 

"17 

PRN  PC  E ( I ) *PPNFCE ( I ) /XML 

<HPpT 

91  9 

JP(I  .EC.  NCONT C )P°M*-CN 

Shppt 

P19 

230 

CONTI NUE 

SHpPT 

920 

234 

CONTINUE 

SHEET 

9 ? 1 

c 

SMcpT 

9 ?2 

c 

SHFeT 

923 

c 

<MepT 

9?4 

r 

4*4*1 

SHEET 

9 2? 

C 

MUOFXxO,  IP  FRICTION  CONDITION  IS  SATISFIED 

SHEET 

926 

c 

*1 , CTFECW ISF 

SHEET 

9?7 

p 

*4**1 

smEbt 

829 

c 

«MFct 

929 

c 

SHEET 

9 ’0 

IPCMtJPEx  .EO.  1)TCHCPP*1.0 

SHr^T 

931 

IPfMiiOEX  .EC.  1 )GP  70  ?001 

prcpt 

93? 

0 

SHEET 

933 

ICHECX  *0 

SHPPT 

0 7 4 

r 

Sh^ET 

9 35 

r 

*4*4  1 

SHCtTT  936  r r,cNF  P *L  I 7FC  NCCA|_  PC^CF  NORMAL  TO  Th?  »UNCh  IS  COWOU^O 

SN^IT  9??  c TO  CH^rK  4HETHEP  TMC  RCUNOAOV  T ? ASSU«FO  TC  *CVE  TOP  PAST. 


128 


sye  ft 

• 1® 

C« 

IM*i 

CMC  er  T 

a-»c 

c 

SHF  FT 

840 

I F I CODE  ( NTOUC  H ) ,EO.  3.CIGO  TP  800 

«Hrc  T 

*4  1 

E 

SHf^T 

84? 

N77i T*N TOUCH- 1 

SHEET 

843 

JP(PF(N7M  .CT.  0.1  GO  TO  500 

SHF  E T 

“44 

lF(*«<l(fF(N77M  .LT.  .000001  .OP.  TTNSTF  ,F0.  O.IGO  Tp  500 

«HcFt 

*45 

TO  I ST  * TOUCH?*. 8 

SHEET 

846 

ESITF  f 6,5 10  ITC  1ST 

FHF'S'T 

84  7 

NTnCCH*NTC(jCH*1 

shf  rr 

HAP 

TCHCOF* 1 .0 

Shpft 

F4Q 

I CHE C*  * 1 

SHF  FT 

850 

500 

CONTINUE 

«H**T 

SSI 

c 

SHEET 

?E2 

c 

S hp  FT 

8«=  3 

c 

SHEET 

054 

c< 

It*  A*l 

S HF  FT 

ESS 

r 

GENERALIZED  nopal  PPOCF  NOOM4L  TO  DIE  15  COMPUTED. 

SHE*  T 

856 

C« 

► #**! 

SHEET 

«57 

c 

S Hc  eT 

858 

!F C CODE (NCONTC 1 .FO.  3.C1G0  Tr  EEC 

SH- E T 

“59 

I r ( NC ONT  C .EO.  NCIF1CC  TO  E50 

SHEET 

8*0 

I*(oon  .GT.  — • OOOCO 1 1 GC  TP  EEC 

SHEET 

861 

c 

SMC  FT 

862 

TCOKTC  tTrONT?* . e 

She  ET 

863 

w°lT:f C,*4C ITCONTC 

SHEET 

8 64 

NCCKTC*NC0NTC-1 

SHC  - T 

865 

ICH^CK  * 1 

S HE  e"T' 

8*6 

*50 

CONTINUE 

SHF  c T 

8*7 

c 

SHcfT 

868 

k*  **i 

SHE  FT 

866 

c 

ICHBCK*1.  IP  F0UNC4RV  ASSIIHOTION  KEEPS  TP  nr  mODIPIcD. 

SHF  E T 

870 

c< 

******  ******** ******  ****** *************** 8 ******** ***************** 

SHEET 

871 

r 

SHEE  T 

8T2 

IPilCHECE  .EO.  1 ICO  TC  446 

SHEET 

873 

NNCPNC *NC ON  T C 

SHC  E T 

6 74 

C 

SHF^T 

675 

C< 

»***! 

SHEET 

8?6 

c 

M4<c  BOUNC4PV  ASSUMPTION  *0°  Kck  STcc  PASE^  UPON 

SHE  E T 

877 

c 

THE  DISTANCE  A A A V FBOM  PUNCH  f.  F DIE 

SMFET 

878 

C' 

SHEeT 

870 

r 

SHEET 

cop 

DO  560  ! * 1 ,NCN  EC* 

SHE  eT 

88  1 

I * ! CON  T A <*  ! I 1-TCONTC  1 ef  1 ,56?,  *6? 

SHE  c T 

8*2 

66-  1 

KCCNT  CsEKCPNC* I 

5Hc-7 

8P3 

7CONT  2* CON  T AC  f I ) 

SHEET 

884 

560 

FONPCP ( NCONTC  lePONF Cc < NNCONC  » 

SHE  E T 

005 

*6? 

CONTINUE 

SHE  E" 

8 86 

C 

SHEE  T 

8*7 

r 

SHE  CT 

oop 

C 

SHF  E T 

ao  g 

r 

«MCCT 

“CO 

C 

SHFE’ 

«ci 

NNTChsN TOUCH 

SHC  C J 

80? 

D"  240  I*  1 , NC ► E CX 

5 Hee  ▼ 

803 

I P ! TCUCH<  ! I-TDI  ST  » ?4r  ,246,?4f 

SHEET 

oca 

?4  6 

NT  OUC  H * NN  T C H—  I 

SHE  c T 

OC8 

TOUCH?* TOUCH! ! » 

SHEET 

8 06 

P c NFC  c ! N TOUCH  1*PPNP  CE  (NNTCHl 

SHEET 

807 

540 

CPNT  I NU  e 

SHEET 

■ eg 

246 

TCHC08  *1.0 

Ehcct 

600 

c 

4MEE  T 

000 

»*  **< 

SHEET 

901 

c 

COMPUTE  TOTAL  STRAIN 

SHEET 

002 

r* 

k***i 

S HE  c T 

003 

r. 

sneer 

004 

WP I TE  ! f , 1 04  3 1 

«HC  F T 

005 

00  444  N *1  ,NUMEL 

SH-  E* 

006 

OO  447  1*1*  4 

SHEeT 

007 

443 

T C°S  C T , N » * T £P5  f I,MAFPSCI,N)*TEST 

chcct 

<50* 

I p ( E SNO  C m , u T . PL  IMTT ) THl CX CN 1 =THIC* (N ) *E*P( EPS!  3,N  )*TrsT 1 

SHC  eT 

0 co 

*9 I TE ! 6 , 10071N,  ! TFP  5 ! I ,N>  , I *1  ,41 

SHC  e r 

0 1 0 

4 A 4 

CONTINUE 

5 Mc  C T 

o 1 1 

r 

SHC?T 

Ol? 

c 

fWC  c r 

SIT 

c 

SHprT 

014 

c 

SHC  c t 

015 

z* 

«***i 

SHC  cr 

016 

r 

COMOUTC  OPNHF  P POC  6CKT  STEP  AhICm  AfULO  GTVF  ^ hF  JNC»EM FNT 

SHCfT 

«17 

c 

np  maximum  EPPFCTIVP  fTPAIN  ap  PC  0 x I m A T EL  v EOUAL  PRESET  VALUE. 

SHF  E T 

<51  8 

c* 

SMfCT 

019 

c 

SMC  e T 

O?0 

CMAXsO .0 

SmcpT 

9?1 

00  ? 7 5 N* 1 , NU mc  l 

SHEET 

o?? 

77S 

IP!-P<!4,N1  .GT  . F M A x 1 FMA XaEPe f 4 ,N 1 

129 


1 


EHCC-r 

SmE^t 

cviCfT 

SMeET 

SWEET 

SHEET 

SHF?T 

sheet 
SHEET 
S uc  eT 
SHEET 
SHEET 
SHEET 
Sh*E- 
ShEct 
SHCCT 
SHccT 
SHEET 
SHCET 
SHEET 
SH”r  T 
SHEET 
cwE^T 
SHEET 
SHEET 
SHff  ET 

ch<"Et 

SHEET 

«h?ct 

SHEET 

CHCFT 

SH?eT 


97<* 
q?? 
97  R 
Q-»Q 

qeo 

qsj 
OR? 
qr3 
q °4 
qps 
q®*> 
q °7 
qpo 
qpq 
qqo 


cw«cc  » 

9-73 

E^  4CT4  s: CONST /r  vt> 

SHEET 

q?4 

JcJK'STFr  • LE  • ? • 4 NP  . IT'T*n  . SO.  OlGC 

S«ECT 

q?5 

Of'Khrn*?,/tFifT  A/l/7  4 N«  l**NP  )-  1 ./'»  4CT  /rrN<  ’ r 

SMCc-r 

q?6 

0°NH“">*  | . /RO^KCp 

CM*-  =T 

q ?7 

CCT  |Os  r)ONM  = c /L7  4 NL«No  » 

CMC  c T 

s 

0°N  e rz  *L  7 ( KU«M  1 

SHF«=T 

q?q 

T7«» 

CON  T T ML” 

C Hc  c T 

qio 

cc  *CTsEF#CT 4 

C wr  C ▼ 

q ii 

** 

C MFCT 

032 

C 

chcc-  t 

97? 

c 

SH*cct 

q -*4 

77-* 

C*^N  7 I NUC 

SHE  ST 

o ■»* 

IF(F5NoCm  .C-t  . PL  |H  TT  ) TC»-COc*  1 . 0 

ShFct 

976 

r 

CHr'T 

037 

Ci 

• •**■ 

SMErT 

9 7P 

r 

K=  k CPHC  IGUE4*  iru. 

cwFET 

<3  7<3 

C> 

**»*' 

c M=  c T 

04  0 

c 

CHCC  T 

04  1 

Op  4 30  1 *1  , MMKC 

5^=T 

Q42 

I 7*  ?*  I - 1 

CH-  E T 

C43 

1 Or  T 7-1 

SHE  = T 

0*4 

04  I |iO(  * »*ua< I 1 * TE  S3 

cmF  — T 

045 

?(!»=?(  1 l*U?(  1 1 *T  S S T 

SHEET 

045 

A ?C 

COK TI NUr 

SHEET 

047 

c 

5«ceT 

04« 

c 

S Hc  ST 

040 

r 

SWFCT 

050 

r 1 

***»! 

SHE E T 

051 

C 

°UMfH  thc  SClUTICK 

FMFCT 

0 = 2 

C1 

**  J 

SHCET 

053 

c 

S hFct 

0S4 

l=4NOUNCH  .CO.  C)  GC  7C  310 

SHEET 

055 

? 307 

C3NT INUE 

SHE0  T 

055 

OUNCH  10l7#(UO4T).i»7CIl,  CC^EC  I ) , ! * 1 ,NU«NP  ) 

e hcct 

057 

PUNCH  1017,  ( 04 I) , T 4 J ) , ! * ] , NL^KP) 

fhFCt 

055 

PUNCH  1017,  4 !TECSn.N»,l*1.41,  »:i),  NUV*L ' 

S Hr  CT 

030 

PUN  <“m  1 0 1 T , ( TM  I E K f * 1 ,MjK:U 

SHEET 

950 

PUNCH  2 223.0NHE3,  NTPMCH,  TChcop  f ccr  4C-  .NCr^’TC 

cmEE  t 

9M 

«\>NCH  235#  IEOHS  C“  t K > • h*l  , NU»M  > 

ShEpt 

0*2 

PUNCm  2"*3fcej^r  , r»c>-  ETK  , OPNmcP  ,Tr«1 

shE~  T 

OS? 

3 1 C 

CO  N’T  1 NUE 

Sheet 

054 

r 

ShE  ET 

OC5 

!P(ES  MO=m  ,GT,  CL  I h IT  ) GC  T“*  2300 

cmc  :r  t 

055 

C 

SHEET 

0*7 

r 

SHEET 

950 

KC’tc  (5 , 1040  1 

cmC  crT 

0*0 

Pn  P40  T*l,NUHN,o 

SHEET 

OT0 

10  = 3*  I-? 

cmcc  t 

971 

I 7*  TOM 

SHCCT 

972 

tL=I7*l 

SHEET 

073 

WO*TECS,  104  11  < I ,CFf  IC),FF(I7),ff,  II  > > 

s HE  e t 

OT4 

840 

CONTI NUE 

SM-Cct 

075 

r ■ 

c * *r-  :***►***** 


£lJV«r  sO  • 

P«r  T=?fNUWNC 
I 7***3- 1 
10=17-1 

8 50  SUMCsSUMF+FP ( 17  )*u 7 ( I HPP ( I 0 ) *U® ( I ) 
«UMPsSUMP/U7(  HV.IHKC  4/TEST*? .*7 . 1 4 

K°  ITE (5  « 10*?)SUMF 
TEfNSTPC  .EG.  11  GC  Tn  2300 
W®ITC(6,  1025  1PPNHEP *NN7CH 


001 

*STTc(«,fio2  9 ) EH  AX  ,UT<  NUMNP  1 , CE  ACT 

002 

jCfNCTCD  ,LT.  Ncgq*)  GC  T C ?1CC 

003 

2700 

CPN t INUE 

c 04 

2 3 0 1 

CONTI NUE 

095 

C 

905 

C 

997 

100? 

ppdmaT(I«,3C12.7,?X,3E13.7,5X,3C13.7» 

095 

1 007 

FTC VAT (IT,  1 IF  1 1 .6  1 

0 99 

l 004 

f 0G  v A t ( 16x51 

1000 

1005 

PPCMATC  in  1,4  ST94  I^.-ST0C  SS  «PVA  7!  CN  AT 

STE» 

N\l*Rcc 

S*  , 14// 

1 001 

l * eL • NC, • .R-STC A IN . . , 7-STOA !►  . .TH-STO AIN . . ,c 

7-STPA  IN 

• • ,EP—  STc  A I N 

1 002 

?...o-STCcS...7_cTce«... Th-stcp?. ,.=?-ST 

S?5..rF 

-STCrSt. 

A vr-STRES  • , » 

1 007 

? 1 

1 004 

1 006 

FpPVAT ( ///  30X , A OISCLACEHCNT  STUlTION  AT  IT  £D  A7  I CM 

NUMOfo  r 4 , I 4 

1005 

1///  20  X , * OLPTUOnef)*,  25X,  * TOTAL*, 

?0X  , * 

0EC  C°HpO 

CCCRP*/ 

1 006 

?/  * NP  CU  OW 

o°eta 

U 

1 00T 

■*  «r  PETA  9 

7*  » 

130 


SHC-J-T  loop  1007  f .1**  *A  T ! |h  ] ,70*  ,*  |TrF*T  JfK  cijPCPef  FCt  ?Tf  P»,  I*  I 

100Q  100*»  rn  6M|T(  to*.  * TOTAl  c-LOAC  *»,  F \ 2 • 7 

SMcrx  10tO  1 / <*.0*.  • TPTAl  7 — L~  A r‘  = . f 12."’ 

inn  ? / *CX.  » T r T A L •'-ITA''  *r  f|2#71 

SMCCT  | Ol  ? 1010  FPOM»T|  //  • K'pni  L Pf  FCC.C*  STEP  -•  , 14// 

•H^ft  1013  1* N.R ..t-FfSCr. 7-*C=  Cc  . . • • . .7-STBFS5  ON  DIE  5U? 

NHFFT  1014  ,c4Cc...*) 

ch^'t  i oi  ^ toil 

loi/  loi?  foaM»rnc,^n,5) 

SMt -▼  | 0 1 t 10!f  F0DM1T  l^px,  • VtLPf  !T»  CONVF  ■>r.c\Cr*  t/ 

SMr"T  101°  1 NO*  , * KTC4  SOUl’  ! CK  VfCrOD  **.  F11.fi 

Shfc-t  101Q  1 / 6TX,  * N^RV  Oc  FCBOf  VECTOR  **  , Fj3.fi 


SHE  E T 

l o?o 

2 / e>0X,  * FRACTIONAL 

NORM 

* • , F 1 3.fi| 

Shfft 

102  1 

1016 

f^omaT!  * f T SPLACF  MpNT 

SOLUTION  AT 

1-^1  PAT  ION  NUMBfF  * * , ]4  ) 

SHE  ET 

1 02? 

10)7 

format (flFJO.71 

SHF  p r 

1 023 

101  B 

FORMAT!////*  op cc  kf*  CCNVcCGf*// 

SHFFJ 

i o •>* 

1*  TRY  AGAIN  WITH  DECELL^RA Tl ON 

CFf  FF  I C IE NT 

*ACCE5*  LESS  THAN*. 

«H”F  J 

1 025 

?FA . ’) 

SHFF- 

10*£ 

1020 

FORMA  7 ! per  4.  1 I 

SHFFT 

1 0 27 

1C  25 

format ( Ax,  15. 3 x ,F  1 2 .6  , 1 c X , ! 5,  3 

x ,»=  12.  t • 1 OX  • 3 

fs.ax.Ei?.#) 

Shfc-t  1 0?4  ??  5 1 FpqK*T(  n,4F20.7l 

tourer  lO^O  1 0?  5 eO  B M A T ! / / / * I N CC  E *F  N T A L 5 TP  A l N — TFT  AL  STRAIN  »'r  STpP  Mlwf'C**,  M// 

chcft  1 030  1 *CL  no.  . . .S-ST  O A IK’ • . . . . ,tw--c;tc  A in  . . • . . . .-M!- STPA  I N .E^-STCAN 

CHcrsrr  1031  1*1 

5HJFT  1C-*?  1 027  FO°vU(///*EL.  fiF  . . . . S-  ST  « ESS  • . . .T  HE- STC  E*  S . . . . Ec-  $TOc  < S • • „ • * ) 

Smf-t  1*33  1 0 A ? FnoMA^i*  PUNCH  FPOCE=*,F15.7| 

cmc-c-t  lOH  1 04  ■»  fnc  v at  < /// *EL  NO  . . . .S-^T0  A IN  . . . . . .Twr  . 5jc  A T N th  I - $T  C A I N . . . . 

SHF-T  10'‘*  lE’F-fTfAIN...*) 

«H«rcT  10^*  1041  FQOV  AT  ! 5 * , l 1 C , E * . n ?0  .7  | 

twf'T  1 037  1040  rnoMAT!//*  NF.  -F  NODE  POPCE*-* 

Swr«rT  1030  510  T<//»  KTOIJfH  IS  E HP  CEO  TO  TOUCH.  fOMOLITf  AGAIN 

s r*re  T I C 5 I*,/*  TTIST*4  , F 1 C .7  1 

SM=cr  1 A A 0 1030  foovaT(/I7,3F1C.S» 

«’"c£T  1041  10^1  c07maT(*  GEOMETRY  r>F  PPOcILe*// 

cmcb-t  104?  1 oc-  l NNO ANCLE...  ..THICKNESS *) 

S^P?T  104  * 104  FPSMATf  / /*  ACPEP  CALCULATE''*  »F  J 0.7  ) 

SmcET  1 044  265  ^CF^ATl//*  CHFC*  INC-  ClSTANCr  A * A Y Fonv  DIE*) 

cnee--  1 04*  270  fq;wa  T U ? ,F2C  . 7 1 

S*-FCT  I r*4/  ?Q 0 0 Fncu»T(/|  |0,f?C.7) 

e^'-e-r  1 047  ->«0C  F^SMA’I//*  CFFCKfNr  DISTANCE  A A AY  FO^m  PUNCH*./ 

SMPcr  1040  1*  ELFM  NO.  TOUCH*/) 

fKFeT  1 040  '*100  *HCt/*T(*  NOTE  AT  NCNMl  IS  TN*IOE  FUNrH,  CO«3  AGAIN*) 

SHEET  1050  34  0 C format!  /*  NODE  AT  NCCNT C ♦ 1 IS  I NS  I OF  ^IF,  COMP  A'AJN*) 

SHEET  1051  3450  FORMAT!/*  NDIE  NDPP  HAS  BEEN  BROUGH T I NTC  CONTACT  *!Th  CORNER*/) 

SHEET  1052  231  Fflc**T(//4  NCCAl  OCINT  COEFFICIENT*,/! 

SMf  ft  | If  3 540  PP3MUI  /*  NCCNTC  IS  ECRCcD  T*~  TOUCH,  T CC  NT  C *•  . * 1 O ,-  > 

Shf  e T 1054  ?•*?  FC^?*  AT  < I \0  ,F  i 0 .*  ) 

ewee-r  1055  P33  FDPMATC  *F  15. 7) 

SH* ET  10E5  ? 3 * FOCMATC 4F  15  .?  ) 

SHE e T 1CE7  ???■»  FCBVATfFis.r,  15,2^15.7. 15) 

SHEET  10**  102P  FOP  M A Tf  * PUNCH  HEAT  0 I SOL  A CE ME NT* , f j o • * /*  NTCUCn*  *,I*> 

* HE  F T 1 050  1020  FOOWATC/-V*  max  EFFECTIVE  STRAIN  I NC RE mf N t * * ,F  j o . 7 /*  Punch 

Shf  er  t 1 060  1 Hf*n  T NC  PE  M£  NT  * * . f i o . 7 . /*  PUNCH  HEAD  ADJUSTING  H^OtcN'l  ' . "* 

SHEET  J061  ?> 

SHE  ft  106?  C 

SHEET  1063  RftlRN 

SHE*tt  i 064  EN p 


SHEET  1066  SUBROUTINE  ST  I ff  ( c , ? , lR , U ? , CODE . EL 3P , y V , Y x , S«m l , C Rh I , Hi  , f p 3 , 

«HE*T  1057  1 TH  I CK  , ALBHA , G AMM A,  ET A, FDNFCE .FF  , A ,B ,N0 1 

SHEET  106B  C 

Shffj  10  6B  r*M**3***************4*********** ************  *******  ******  ***•••*•****• 

SHEET  1070  C CALCULATION  OE  STIFFNESS  MATE!*  E TO  FNTIOF  SvSTFM 

Shf  ET  1 071  c*** ***************** *********** *********r* *******•»*••**••*• *•**•••**• 

SHE f t 1 07?  C 

SHEET  107  3 rOMMQN/r-FNCON/NUMNP.NUMEL  . H^(l?)  ,PLL  • NFC,  NFQPM,  VJELO.TEST,  I T Ec  . 

SHEET  10  74  1 NRF  An,  NOUNCF.Nro  fNT.R  VALUE  » T,  *B  AND  , P NF  AC, RADIUS. FPI  TigO  ,FO  J TNC  . 

Shf e r 1 075  2ECONS-, CNHEC, PT HOT , r T EC AD.TCTNTC. TOIST, PHFCE 

SHEET  1076  MMfnSION  Q ( 1 ) , 7 ( ) ) .CCDE  ( 1 » . UC  f 1 > VU  7 < 1 > . SLPO  f 1 ) • **  ( 1 ) , 4 « NO  , I ) . 

SMSpT  1077  1 FCS(4, 1 >,PR(2),Z7(?), UU(6J,YY<  1 >, THICK!  1) .PLl  1)  ,fPH! ( I)  ,CPHJ  ( 1 1 , 

SHEET  107«  ? YX ( 1 > . ALPHA!  1 1 .GAMVA ( 1 ) ,E- A ( 1 1 .FCNFCe ( 1 > . EFf  1 J 

SHEET  1070  CPMMON/STFM AT/H! € ) ,D! € , £ » ,TE X , TE V ,TE? . TM» L 

SHEET  10P0  COMMPN/CPNOUAC/SS ! A ) ,WT ! A ) , 0! 2, 2 1 .SOFT  1 

S HF  FT  1 OB  1 COMMHN/ATOLCH/NTOLCH.NCIF , NCCNTC 

«HEET  1 OB?  C 

SHF  FT  10B"»  C 

SHFCT  10  34  r*4****A*44*********«4*3«  ******** ** ***************** A****************** 

SHF  e t 10*5  C INITIALISE  a ANT  P AATCIK  fqc  FOUATIPN  AXtO 

SPFc-r  10B6  C BE  c A U eP"  BANDET  SVMMFTPIC  ocCOfr^V  OF  THE  STIFFNESS  MATPI*  A, 

SHE e t 1 0B7  C THE  STORAGE  OE  A T«  IN  A SOUARF  MATRIX 

«Mp -T  1 0°B  f ********  **********  *********  ******  **************  *****4**  * A************** 


131 


chF  FT 

1 0*9 

r 

SM£C  T 

1 090 

or  «o  n*i  , kf  c 

SHE  FT 

1 091 

* C • 

SHCCT 

1 09? 

pn  Msi.mranO 

SMFCT 

109”* 

FO 

AC  K ,M) »0# 

SHF  = T 

1 

WT ( 1 > = 0. 347Ef 4P4E 1 

SHe  F "* 

l nos 

• T l?l*0, 652 1 AS  1 549 

SHeeT 

1 0C6 

• T ( 3) s»T<  11 

SMFFT 

1 09T 

«T<4»=*T(2  » 

CH'CT 

J OSS 

r 

SH?FT 

1 099 

c 

cmF  c T 

1 1 00 

c 

**«*< 

SM^F  T 

1101 

c 

CO*.'  STeLC  T P ANT  H A t FLEVF^T  L E VF  L 

Smc-ct 

1 1 02 

c 

*44* 

cmFET 

1 1 03 

c 

SHcFT 

1 1 04 

OO  100P  Ns  1 , NUMFL 

SHF  FT 

1 1 OF 

NCl =N+ 1 

cHEfT 

1 1 OP 

IF(FOOF(N)  .£0.  3.0  .ANO.  rOPECN^ll  .FG.  3.01  GC  TC  1000 

SHF  FT 

l 1 07 

0LL*OL (N  ) 

5HF- T 

1 1 OP 

COM  *COh 1 < N) 

SHFFT 

1 109 

COmsTOHI  IN) 

c w«F  T 

1110 

QO( 1 ) >Q (N  ) 

SHF  - T 

1111 

72  f 1) =?(N) 

SHF  FT 

111? 

PR ( 2 1 = 0 CNP1  ) 

CMC  try 

111? 

LUl  1 » *UR  IN) 

chEct 

1114 

UU< 2) *U7» N> 

SHF  FT 

lllF 

UUI  »)*SLCP<M 

Smc  FT 

1116 

UUl  4)  *U9  f NP1 1 

Shcc^ 

1 1 17 

UU(*I=U7(NP I ) 

chcpt 

11  1 * 

UUI6)  *SLCFC NF1  ) 

SHFFT 

1119 

THKL*TMICK<  N) 40 ll 

SHF  F T 

1 l ?o 

ZZ( 21*7<NP1 1 

SHFFT 

1 1 ?! 

YG*VX<N> 

sheft 

1 1 22 

VHs VY (N  ) 

SHFC  T 

1 1 23 

c 

shf^t 

1 I ?4 

CALL  OUAOI  PP  . ??  ,UL  «CLL  « .C  PH  • Vf,  , VH  ) 

SHFFT 

1 1 ?S 

r 

5ht;^ 

1 126 

C 

SHErT 

1 1 ?T 

c 

SMCCT 

1 1 ?P 

f *******************  **************  *ft  44  ******-****-*•  ******  ****Sr»r4**T* 

S^FFT 

1 1 29 

c 

PCCC09M  the  ASSEmoL  V OPE  ° A T I ON • BECAUSE  MATRIX  A IS  SYMM^t^IC 

CMC  C T 

1 1 30 

c 

CNLV  UPPER  HALF  rc  THE  FATCjv  IS  CREATEO.  ANO  THE  STORAGE  pps 

S HC*T 

i I 31 

r 

matrix  a is  a souaff  appav  fecalsp  OP  ranoft  symmetric  opcpectv 

SHFFT 

1 l 3? 

C 4*  ******  4 *44*  ********  *******4  » *****  4 444  4 4 4** 4 **** **4 **4 **************** * 

« MC  C Y 

1 1 33 

C 

SHFFT 

1 1 3* 

00  ?O0  I*If  6 

CMC  p Y 

1 1 

II*N*3  - 3 4 I 

SHF^T 

1 1 -*e 

«M  I ! I**M  T I I4M  ! ) 

CMC  r y 

1 1 37 

or  2oo  j*if  6 

‘HEFT 

1 1 ■»« 

r 

SH'Ff 

1 1 39 

J JcN* 3— 3 4 J— I T 41 

SHFFT 

1140 

TP  < J J .LT . 11  GC  TC  200 

SHEffT 

114  1 

* ( 1 1 *jj >*a< ! i ,jj>4Pf y .j) 

SHFFT 

114? 

?00 

CONTINUE 

1 14? 

1000 

CONTI NUF 

SHF  rT 

1 144 

C 

SHC  £ y 

1 1 4F 

C 

4*4*4 

SMfrT 

t 1 4P 

c 

TO  HXMDLS  MIX50  OrilNOAOY  CONOITION,  f=OLLO»’!N0  HTJIC'S  1»E  EVAlUiT' 

C M«  C f 

1 1 47 

r 

***** 

SHFFT 

1 14? 

r 

SHF  ^T 

1 1 49 

C*ono ao-pnhEC 

CMF  c t 

1 1 FO 

NMFl rNUMNP-1 

s^r 

1 1 FI 

PONRA0  *PNRAO 

Shc  pt 

1 1 5? 

R 0 I e * RT  HPT  40  I EB  AO 

T 

I 1*3 

c 

SH?  *T 

1 1 S4 

00  1200  N*  1 ,NVP  1 

SMFFT 

1 1 55 

I* ( CCOE ( N » .NF.  6.0IC0  TO  1200 

SMcey 

l l«P 

OUM  ?*POtE-P(  N)-»UP<  N) 

SMfFT 

1 1F7 

OUM 1*0! F°AO“ 7 f N )-U7(Nl 

cmc  c r 

1 1 F* 

DNP  AD *0? f RAC 

SMFFT 

1 1 59 

?P(N  .LF.  NCONTClGr  TO  1100 

cmc  r t 

1 1 60 

CNR  AOcCCNR AC 

1 IPt 

0 U*4  ? * P f N ) ♦ UP  f N) 

SHC  r t 

1 16? 

RUM 1*C4?CN> 4U7 (N1 

«MC  f? 

1 1 *3 

r 

SHFFT 

1 166 

1 100 

CONTINUE 

SMFFT 

1 IPS 

ALPHA (N1 *-CUM2/CUMl 

SMCCT 

I 1 P 6 

GAMMA  IN)  *<ONR  Ar*PNPAO»OUM?4-OUM?-OUMl  *PU4*1  )/2./0 UM2 

SMcrj 

1 167 

FTA(N)cCRNFCECN) 4PNRA0/OUM? 

Smcct 

1 1 P« 

I * f N #GE.  NTOUCN )GC  TC  1200 

SMFFT 

1 1P9 

GAMMA |N ) *-G AMMA 1 N) 

Smfet 

l 1 70 

FTA(M)i-fTAlN) 

SMFT’T 

! 1 7 | 

1200 

CONTI NUE 

Smcct 

1 17? 

C 

SHC®  T 

1 1 73 

c 

***** 

132 


«Tr*o<?  G^nFRALIPEO  NCOAL  EQPCE 


m 


SM-FT  * 174 
SHEET  1175 
SMf®-T  I17f 
St-EST  1 1 77 
«H«fT  1 1 7P 
SMPCT  1J7V 
SMFFT  1150 
SHF^T  11«1 
SW^FT  115? 
SHE  FT  l 1 " 5 
SWEET  1154 
SM*«TT  1155 

SHpfT  1 1 P6 
SHF?T  1157 
S«F FT  1155 
ShF'T  115Q 


MU*'  'iNllWND*  * 

r>C  1500  1 *1  (NU»? 

1 ■■00  55  III  —51  1 » 

C 

C 

f 

c 

1001  •»r3M*t(///,»  T WE  ''lAGPKAL  V^rrec  c E MATejK  r»*  STI^MFfS*/) 
100?  POOVtT ( x ?ei 1 .31 

1005  eOU^K//  ?VH  FL5VFA-  ft  ITW  NEGATIVE  AE  FA  ».  7£  » 

C 

0*TU0N 

FNO 


SHCfT 

1 I PI 

SU°0 OUT  7NE  CUAPIPR,  77*  UU.  DLL  .S  Rw,  C«>w,  S2.  S 1 1 

SW5CT 

1 1 O? 

COMMON /STRM AT /w (61 ,0(6.6) ,TF*.lFV.Tr7,TM*L 

«Hf?T 

1 I 03 

COW5PN/C0NCU AC/SS (4),ftT(A),0(?,7» , SCF ▼ ! 

SwE  E T 

1 1 04 

PI  MENS!  CN  oc  ( ?>  ,E7(  ? » f CU46  » .5(  2 .6  1 ,2*  (6  .6  1 , «‘7f  Of'(  6 1 .0*  ( 2, 6 1 

SWEET 

I I®* 

common  / I S'lTPv/P  VAL  uE 

SHEF  T 

1 l Sfc 

C 

SWE^T 

1 107 

9f«fB9| 1 •♦©»» 2) )✓?• 

cmFET 

1 1 OA 

c 

SwcpT 

1 1 CO 

or  2 1*1,6 

S*-ceT 

1 200 

W(  I 1*0  • 

SHCC  r 

I 20 1 

OC  2 J*1  .6 

SHEET 

1202 

2 Of  I , J 1 *C. 

SHcfT 

1 203 

TE*  *0  • 

SHEET 

1 20a 

TE  V * 0 • 

SWEET 

1 205 

TF  7*0, 

CHF  c t 

1206 

c 

1 207 

c 

SHC^t 

l 205 

07*77(2  I-77C  1 1 

SHEET 

1200 

ric*cp  ( l i-Po  ( 2 I 

SHEET 

1210 

DU*UU( 1 l-UUi A ) 

cwFF  T 

121  1 

P**UU(S l-UU (2  I 

SHE  Ft 

1212 

AU*tU(  1 ) ♦ W U ( 4 \ 

SHF  FT 

1217 

AC  sCO  ( J J *QR (21 

cmF  E T 

1 21  A 

r 

SWEET 

1215 

Cl*  2. *D» /DLL /DLL 

SHEET 

1216 

C2*?.*0U/CLL/CLL 

SWE*"T 

1 217 

C 3* 2. *07 /DLL/OLL 

SHE  ET 

1215 

C4*  P.AOft/OLL/CLL 

SHE  E T 

l 21  o 

C«  =AU/AO/2 • 

SHEET 

l 2?0 

C 6*  1.40O4C  ?A07*C4A(nu*OUAOft*Oft  1 /OLL/P-LL 

«HEFT 

l 221 

C7*2./DLL/CLL 

SHFCT 

122? 

C «*  2. /AS/AR 

SHFFT 

1 223 

CO* 1 ./SOFT{ C6 1/2. 

SMFr 

1 2 2 A 

Cl 0*CO/C6 

SWEET 

1 225 

Cl  l «C  14C2 

«M«ET 

1 226 

Cl ?*C7*C4 

ewcc  r 

1277 

r 

Swcer 

1 ??« 

C 

cnee  T 

1 2 ? C 

c 

SMfPT 

1230 

OE  S 1*S0»7( CC) 

SHEET 

4?31 

OFT  1 *2.*C«M  • 

chERT 

1 27? 

o 

SWEET 

1 2.17 

Fl*COftr i l/OESl 

SHEET 

1 2.14 

c2=-cOftCl 2/0ES1 

Swe  FT 

1 275 

E3*-F 1 

SHf«T 

1 ? IF 

SA  *— E? 

SWE«T 

12  77 

FK  * 1 • /AO/DET 1 

CH**T 

1 »15 

ES=(-C10*C11*CI  1/2.4C0AC7)/0FS1-R!*r 1 

SHEET 

1230 

r?s-e 5 

eweer 

1 740 

E5*  C 1 0*C 1 1*CI 2 /?• /PE Sl-E t *E  2 

cmc  c t 

1 241 

FQ*-F5 

* we-er 

1 24? 

El  0*Fa 

SHCC-r 

1 741 

C1 1 *(-C104C124C12/?.4CSAC7>/CESl-E2*E? 

CMCCT 

l 244 

F I 2 *-E E*EE 

e>.epT 

1 PAS 

c 

OE?*MpOTriAN  STRAIN  JNCPF^FNT 

0*T*C  fFCt'MFEPfNT  I AL  STRAIN  IN‘CRFMFNT 

rnMctjTATjON  or  erFRC*iv£  s^oatn  INCRF*Sn* 

FleCFPfVATIVc  PR  MF&7DIAN  «TOAJN  INCREMENT  ft|*N  RFfPECT  TC  UU(I  » 

*c  c cf  s i /c  cm  nil 

“2*01  O- S >/D(  UL ( ?)  ) 

F 7*  0 ( r*FS  I/O  f UU  ( A 1 1 

c4*onpp»/c  ruirr>> » 


?4|S 

24  T 
245 
2AO 
2FO 
2*1 
7*2 
?*1 
2?  A 


13.-5 


<HFC-’ 

«w|crT 

SHTCT 

S HP  FT 
SMFFT 
SMFF’ 
SHFeT 

CMfP  ” 
SHCCT 

SMEC^ 
SHFCT 
SHC -▼ 

sheet 

SHTFT 

cMcpr 

enecT 

SHEET 
SHF  FT 
SHF  FT 
«H«eT 

ch'CT 
SHEET 
«HEEt 
SH^ET 
?mEFT 
ChF  e T 
SHEET 
$He  F T 
SHECT 
SHEET 

£H«'F*r 

SH*?FT 
SHF  FT 
SHEET 
SHF  FT 
SHF  c T 
SHFFT 
«HPpT 

SHFe"r 

CM-  «T 
SHFFT 
SHFFT 
SHFFT 
SHerT 
SHFFT 
SHFFT 
S HFeT 
SHFeT 

5MCTT 

SHFFT 

SHFFT 

SHFFT 

SMPCT 

SHCCT 

SHE  c T 
SH^F'r 
SHFFT 

SHFFT 
SHFCT 
SHFFT 
SHFFT 
SHFFT 
SHFFT 
SHFFT 
SHFC  T 

SHF  FT 
SH^FT 
SHFFT 
SM^FT 
SHFFT 
S HFe‘T 
SMeet 

FMf  f t 

SHffCT 
SHFFT 
SHF^t 
S«F  c T 
SM**1  T 

SHFF’*’ 
«HFCT 
?H=e’ 
Sneer 
SHC  F T 


?5» 
2 5* 
?«•* 
?5* 
?«c 
?60 
2*  I 
2*2 
263 
2*  4 
?6F 
266 
[ 267 
26* 
269 

?70 
2T1 
27? 
273 
2 74 
27b 
?7fr 
?77 
27p 
?79 

2 SO 
2*1 
7*2 
2*3 
2*4 
2*F 
2*6 
2*7 
2** 
?*<J 
?90 
7 9 1 
?Q? 
293 
204 
29  * 

296 

297 
99* 
?QO 
30C 
301 
30? 
303 

3 04 
30F 
305 
307 
3 OB 
309 
MO 
31 1 
’1? 
31  3 
7 1 A 
3 1 F 
3 1 fc 
717 
31  « 

319 

320 

321 
32? 

323 

324 
3 2* 

325 
327 
32* 
3?  c 

330 

331 
7?? 
33’ 
3 34 
375 
33* 
’37 

3« 

33C 


F F5*OI*>F^l/D<ULf  1 II 

C F6  * C ( El  1/CCUUCl  I 1 

C E7*OC  F 1 » /OC ULCA  1 ) 

r F*«CIF I l/DCUUC?  » | 

C C9*D(E3 » /0< UUC? 1 1 

C E 10-DCF41/DCUUC ■» 1 

C 51  1 *C  IF5  l/DU'Uf  ? 1 I 


c 

c 

OFS»4LOGI Ot  S 1 > 

OFT  * ALOC I TFT  1 » 

P V° 1 *® VALUF* 1 • 

PVP?«  SOPT( ? • *P  VALHF ♦ 1 « > 
RVC3«RVP| /F  VP? 

PVDA«9.  ACVALLiE/P  VP  1 


Er  * ECT 1 VF  STPAIN 

pi  *ceb 1 vat rv*  r*  FPF*fTtvr  stcatn  ik cf ehfnt  with  cfs=*fct  to  uuc  i » 

F ?*  A I Th  BFSPFCT  TC  UUC  21 
F 3 * * I T H BESPECT  TC  UU  <41 
p 4 * A I Th  PFSPECT  TC  UU  (5 1 
F l l*r>(P  1 l/OCUCC  1 » 1 
PI  ?*DfPl  l/CCUUC 2 > » 
p 1 ? *0  C p 1 ) /D  < L L C A I 1 
f 1 4 *0 ( p 1 1/CCULM5) I 
p?2*0(p?» /DCUUC 21  I 
P?7*0rc2 l/CCUUC  A I I 
P?A*n(P?)/C(l'UC5l  ) 

/MUUC  41  1 
P34«c* P31/0(UUC5I > 
ca4*D(P4)/P(UL(«) > 


FP  S sOFS*r>E  S40ETA0FTAP  VPAAOE  S*TFT 
PFS 1 *PVP7*S0FT  (EPS  I 
FP  S2rCVF3/SCPT<FPSl  /2 . 
pp  S ?*  — q vp  3/£P  £ /SOP  7 f FP  S I /4 • 


01  * r ? .*OES  VO A* PET  )*P 14  < ?.4DF  T4P VO A^OESl *FF 
02 *( 2.*0FS4PVPA*DCT > *F2 

0 3*  ( ?.*PES*9VQ4*p>FTI*E?*(  2**Dp  T*p  VP4*0FS1*FF 
OA* (2 ,*rcS*PVCA*r!’ >*E4 
*PI * reS?*0 | 

F2*  FFS2*C2 
F 7 *EF  S2  *03 
fa*fps?*oa 

f 1 1 *FPS7*o  1*01  4EFS2*(  C 2 .*OES4CVPA*OFT  >*E54f  2 ••'PET  ♦BVPA*OES  > *E1? 
1 4(9. *F1 4CVPA***  1 *F 1 ♦(?.*EE*°vd«<Ei  »*F5I 
p 1 2*FPS3*P 1*0?4EP  ??*C  C 2 .*0FS4  9VP4 aOET  )*£*+(?••  cl ♦ p VP A * EF 1 *E2 I 
FI  3* EPS  3*01 *0?4EPS?M ( 2 .*0ES4BVP4»0FT»*P7  4f  2.*fcE  3*P  VPA*E  5 1 *£  l 
1 4( ?»*ES4BVP4*F3 )*E*4( ? . » 6 F T ♦ o v p A * C F S 1*F1? ) 

F 1 A =c PS 3*0  1*D44PPS?*( ( 2 .*0" S4P vO4»0FT) *F°+( 2.*F  1 49VPA*F* I *EA 1 
C ( 2 .♦CFS4BVP4A0FT 1 *Z 1 1 4 ? « *E 2*E  ?» 

?,*OFS4SVP4*OETl*F34(  ?.*E34BV04AFF  1*E?  » 

F?4*PPF 3*02*044FcS?*(-( ?,*0£S4P VPAAPFT) *e 1 !♦?•*= A*F2I 
P3  3*EPS3*03*0  34FPS?*(  C 2 .*DES4CVP4  *DFT  I *E6  4 ( ?.♦£ 3*PVO< *E  E 1 *E  34  ( 

1 2. *FF4PVP4*E 3 )*EF4f ?. *DE T45 voa*OE ?) *E 1 2 1 
P7 A = E *S 3 *07  *04  *EpS?*f  ( 2 .*PF S 4PV»4 »OET I *E  1 04f  ?.*E  ?4S VP 4*E * ) *£  4 1 
P44  *FP  S3*04*04  4FP  F?* 1 < 2 .*0= S4C VP4 *0ET » AF 1 142.*E4*E4> 


°r  1 ♦ 1 )*(  I Sl4S?*EF*n*P  1 1 ♦ S?*F  1*P  1 1*PC*THKL 
°C 1 ,2  1* ( «S1 ♦S2*FFS1  >*P124«2*P?4P1  »4BC*Thkl 
P(  1 .41  *<  < S1*S?*FFS1 1*F 1 ?♦ S?*F 1 *pt»*cc*-hkl 
P(  1 .*  >*  f ( S 1 4S ?*FPS 1 >*F  J4  4S?*p 1 *p 4 )*CC*T hk L 
P(  2 .?  1 *(  ( SI  4S2*  *PS  1 1*F??4S?*p?*F?  lACCA’-HKL 
P(  2.4)  *H  S1*S?*PPS1  >*F  2 34  32*P  ?*P3  )*PC*ThkL 
P(2*5  »*  < fS  I *S2*FPS  | l*F244S?*P9*P4  »*PC*TMttL 
of  4,4 »*( C SI  4 S?*  Fr  5 I >*P37*?2«p**P3)*RC*Thkl 
0(4,51*1 C S 1 4S2»ppc 1 )»r?4*s?*P*«r/ ) *C  C*rHKL 
Pf 5.5)  *( C s 1 ♦S?*fPSI  » *P4A 4S ?AC4 4*4  |*rc»vHKl 

p( ?• 1 1 *P< 1,21 
p (A  ,1  )*Pfl  .A) 

P(  4 . 2) *P( 2 .4) 

OC3.1  )*PC  1 . £ 1 
PC  5 %? 1 *PC? .5) 

0(5,4) *°( 4 ,£) 


HC  1 )*— ( S 1 4S?*FF  ? 1 )•*  1*BC*TWSL 
M(  ? ) «-(  SI  4S?**P«  1 )*P?4Cr*TH«fL 


134 


SHEET 

1 ’4C 

HI  4 l*-l  5 1 1 )**  •■•tC^'THKL 

c T 

1 7AJ 

1 i*r4*cr*TH«(L 

1 1*? 

C 

« *.b  cr 

1141 

C 

«MF£  T 

1 ’44 

71  r^KTIMUF 

SHF*’ 

1 '*45 

oc  Tl_°N 

CMSfT 

1 34* 

CNf* 

fhFCT 

1 34« 

SUR’OU'"  INc  rr  N^ENf  A ,n  ,KFO  ,H«4jkKO  ,N  #U) 

SWfPT 

1 3 A <3 

r 

SHf'T 

1 750 

Ci 

k***i 

k»* 

SHF  FT 

1 351 

r 

OBcpnovi  THE  MIT3t»  C^NOFNSA  t1  rw  NHEN  THE  VALU?  OF  A CO*o*"»KFNT 

SHF  - T 

1 *52 

c 

X IK  A**P  IS  «PFCt*leC  FOUAL  Tf  ▼ffcC 

5MCF' 

1 353 

Cl 

k«**i 

M« 

SH£  ft 

1 ’64 

c 

fHfC’ 

1 35S 

0 ? HCN  c I ON  R(NEO) .41 K’C.l 1 

SH^f-r 

1 356 

c 

SHF  ~ T 

1 3*7 

no  230  Ms?.*6ANC 

Mrrv 

2 

F_t_F_Fr 

1 3*« 

x 1 

vrr 

4 

SHF  F t 

1 359 

jciifi  2 3* • 2 35  » ? 3 0 

VO'- 

5 

ShE  F^ 

1 360 

230 

RIK 1 *«(  * >-4<  * .m  >*L 

«:r 

c 

SHF  FT 

1 361 

A I*  . M 1 »0  .0 

Mrr 

7 

SHF  ” T 

1 36? 

236 

K 1N4V.I 

Mnr 

P 

S Hp  CT 

1 363 

fC(veo-K)  25C.240.24r 

Mrr 

9 

SHFF  T 

1 364 

24  0 

RIF  »r»(K  )—  A ( K , M ) e'J 

nor 

1 0 

SHF  FT 

1 36? 

A IN. Ml  *C.O 

Mpr 

1 1 

SHE  FT 

1 366 

250 

CPNT  INI1E 

Mcr 

1 2 

SHFe”r 

1 35’ 

AIN  , 1 »*  1.0 

MOP 

1 2 

SHFFT 

I ’66 

°t  V >*u 

MC^ 

14 

Sheet 

1 ’69 

QF TLRN 

MOP 

15 

ShFF^ 

1 370 

C 

m«p 

1 t 

S^F-  T 

1 371 

rNO 

wnp 

I 7 

SHEET 

l 37  3 

SU6B3UT  INE  MOCIEVI  CHOE  . A , P . ALPHA  .C-AHMA  ,E’A  • AUMNC  . K * G • MR  ANC  • F BNP Cr  1 

SHEET 

1 374 

0 I MENS  I CN  CODE ( 1 1 • A (NEC. 1 1. « 1 1 1 . AL  ®MA I 1 »• OAMma 1 11, E’  A ( 1 1 .FSNFCf  C 1 1 

SHfe-T 

1 775 

on  1PJ  1*1,  nlmno 

cmcct 

1 376 

IL*’*  T 

SH«?et 

1 777 

I 7ML-1 

sheet 

1 37S 

!»*  !7-  1 

Shcct 

1 ’7Q 

C* C COr ( I 1 

SHEET 

1 3*0 

rr  ic.eo.  i.i  oc  th  ioi 

SHEET 

1 361 

IP  fC.”0.  ?.l  GC  TO  102 

JMC-p  T 

1 ’»? 

IP  (C.EC.  3.)  GC  TC  103 

SHEET 

1 3*7 

r 

SHEET 

1 3 64 

CALL  CHNOENC  A.P.NEf).  HP  AND,  IL.O.  1 

SWCffT 

1 365 

PON STA  *ALPHAI  I » 

SHEET 

1 766 

CONST  r * eT  A ( I 1 

SHE**T 

1 3*7 

consti;«gammac  i » 

sheet 

1 ’6« 

IP(C  .EO.  41GG  T3  1C4 

SMB  £ T 

1 360 

GC  TC  121 

SHEET 

1 3 QC 

C 

SHEFT 

1 3 91 

101 

CPNT  INije 

sheet 

l 392 

CALL  CONOEM  A .P.NFC,  M6A6D,  19.0, 1 

SHPCT 

1 303 

CALL  CONOENI A ,B .NFO.mf ANO , tL .C. 1 

SHEET 

1 ’94 

GC  tc  l?l 

SMEE’ 

1 395 

C 

SMB  ej 

1 306 

1 02 

CONTINUE 

SHEET 

1 397 

CALL  COAOcNI A ,P,Nf C.MPANC.  17,0.  > 

Sheet 

1 396 

CALL  CONOENI  A.Ft.NPO.MBANP,  IL.C.  1 

sheet 

1 3 99 

p|  yc)*o(  IRmfcnfcEC  I 1 

sh*et 

1400 

GO  TC  121 

Smb  bt 

1401 

r 

SHfCT 

140? 

1 03 

CCNTINtjr 

SHEE*r 

140’ 

CALL  CONHENI  A .6  ,NF0.  *4>  AKO.  IB,  f . 1 

«MC?T 

14  04 

CALL  CONCeNf A« °. NFC, HP  AND,  I?,0.1 

SMEF'P 

1405 

CALL  CCNOEMA ,p ,Nc0, MBANO,  ILfO.  1 

SHEET 

1406 

GC  TO  121 

SHfE  T 

1407 

z 

Sheet 

1 406 

104 

CONTfNu* 

SHEET 

1400 

CALL  PC*!* ( A, F .NEC, *»PANC, CONST A . CPN*T£ ,CPMSTG, I» .171 

<Hcer 

1410 

r 

SHEET 

141  1 

121 

C3NTTNU* 

CMCC  T 

14|? 

c 

SHPET 

1417 

OP TL9N 

Sneer 

1414 

ENC 

135 


mmrnrnm 


SH^ROIIT’NC  Tri4|MN,VM,l) 
0IMeNS  If!K  A I N K , 1 > 


T°  I *NGUL  I 7 At  I CK  OF  GAUSSIAN  CL  I **  ! H A*  I CG  F~»  T*-?  SOLUTION 
of  PANOEn  CVMMFTOI r VAIS J » 


Nj  0 

100  N *K  * 1 

|F(N.:O.MN|PeTtO\ 

JPIAJN,  ] I.NE.n.)  CT  TO  150 
GO  TO  ICO 

c 

1*0  T*N 

MB  =M  TN0(  MM  ,NN-N*  \ ) 

r 

oo  ?60  L*?.Mp 
I » I ♦ 1 

C*AfN.L 1 / A ( K » 1 1 

I^< C.FO.0.0IG0  TO  ?00 

J * C 

r 

DC  2*50  KsL.MP 

J*  J 41 

?5  0 Af  I . J ) « A < I , J »-C* A (K ,K  ) 

AIM  ,L  1*0 
26  0 CONTINUE 


SU^COMJVTN^  RACKEIK'N.MV.A.oi 


SwC  eT 

1 Af3 

c 

°*0'  SlIP  ST  I TUT  I ON  rpo  SOLUTION  ?f  BANDEO  5W»>c  TCJC  matc]* 

«H«T 

IAEA 

rr+  ** 

SMCct 

1 A s* 

C 

e mc  * r 

) AS** 

C 

SHEET 

1 A57 

01 MENE1 CK  All), Fill 

Si-efT 

1 AS" 

c 

eH=e- 

1 AES 

mmm.mm-j 

SHEeT 

1 A00 

M = 0 

SHEeT 

1 A*>  1 

?T  0 

K*K  41 

SHEET 

MS? 

C *° INI 

St-sfT 

1 A6T 

I F ( A ( M I .ME  .O.C  1 P 1 N 1 *n (M 1/AIM) 

<Mr  - T 

1 A *-4 

IF  | K.EC.KKIGC  -A00 

5«e'T 

1 AS? 

I L *M 4 1 

SHEeT 

1466 

I H * M I MO  1 NN  * N 4 M M M 1 

£Me  - T 

1 A *•  7 

VrM 

IAS  « 

;»c  T * IL  • I ► 

Smc^t 

1 A OO 

MsM4KK 

Shcct 

1 A70 

?S5 

PI  I >=*l  I >-A|M)*c 

CMFPT 

1 A71 

GO  TO  ?7  0 

cue  c j 

1 A 7? 

c 

SH*fT 

1 A^3 

r 

SHEET 

1 A7A 

TOO 

IL*N 

S MF  c T 

1 A-»5 

M«M-  1 

SHEET 

1 A**6 

IFIM.EO.O)  PETHPM 

«M«.T 

1 4*7 

IHsMTMOINNfMAA-MMl 

SHEBT 

1 A->* 

M«M 

euce  r 

1479 

c 

SM*-*-T 

1A«0 

00  A00  I*IL,!H 

* wcST 

IASI 

M*M  4MM 

ct^E  E * 

|A«? 

A00 

«|  M*n  |K  ).A  I P )*«  ( T * 

SHEET 

1 4*3 

GO  TC  ?00 

chcct 

1 ASA 

0 

CM*CT 

1 A*S 

r 

SHEET 

1 ASft 

EM  r 

SH»E*r 

1 A®S 

StP^OcTTME  HAPO(eoC'V) 

CHCCT 

1 ABQ 

Cnwn*4/M»TCRt  / v v v AL  . OPE  STM,  FE  wrMT  ,ocr  fT  « 

5HerT 

1 A Q 1 

C*A*A 

1 4-71 

e**« 

SH'PT 

1494 

c 

SHCjT 

1 4<5* 

VV4LUF  iWVAL 

SHE*T 

1 40* 

5 XDNT *cc VPN  T 

5M?CT 

1497 

VsVVALUF*(PRFF7N4prsi»4'K3»J’*rcc^*5 

•sM^T 

iasp 

RP TtPN 

SH?*T 

1400 

FNT 

1 501 

SUOPfU-TINE  H4Rfi?(5C«,Yl 

«M?eT 

1 50? 

CO**OK/MATcGL/VYV AL.PPESTN, EE  XCNT.DFf STS 

«Hee* 

1 50-» 

C 

* 

sheet 

1 504 

C««i 

SHfPT 

1505 

c 

CC**UT*  *CF*  ►icrEMKG  FAT® 

SME?T 

1 5 Of 

<:♦*« 

SHCCT 

1 *07 

c 

SHEA" 

ISO® 

EXPM*>p  *PM 

S*-?CT 

150S 

WALU5*  WVAU 

SHrrT 

1 51  0 

yaEXPKT*TVALU?* t OFFS^NApPS  I* A (P»PNT. \ . » 

S HEAT 

1*11 

RET1.3N 

SHEET 

1*12 

cNC 

S ME  *t 

1514 

S«jF*»rUTINE  «CM  IV!A,5,NEO,H0*5:C.ErN«TA.CrKSTA.CCNSTr,rA.F7l 

cmcpt 

151* 

OIMENSTCN  A (Nfo, ! >• PC  1 1 

SMtr“T 

1515 

C 

<«F*T 

1517 

c* 

***** 

SHE  c T 

1 31  P 

C 

THIS  SlV*»OCUTINF  HAKXLFS  TH:  Hjvfn  PCUNCAOV  COKOITICN 

«;MrpT 

1510 

r i 

1570 

r 

SHEET 

1 S 2 1 

S^ET 

1 5?? 

A(  ! ?,  11  *A(  I 2 #1  >42.5A(  ie  , im  /CC5S7-A4AMC  , 1 > ✓ C^NS  T a ZCPN«  T t 

SHCPT 

1 *23 

OC  I7I*«C  171  4F  f rc  1 /C^NATA-C.7VSTO*IAf  I C , J ) /C^N  *T  A 4 A C I C . J ir  » » 

cMfE  7 

1*’4 

PIT  71  *P(  I 7 MCCNST* 

SHE  2* 

!«?« 

4< IF, 11*1.0 

<H«?T 

1 326 

A ( r C.  r<  1*0.0 

€M«* 

1527 

SC  I Cl *0. 0 

S*-e5T 

15  25 

C 

SHPCT 

1 

c 

SHPPT 

IE  30 

or  300  M*2 ,np ANC 

SHrc-r 

1511 

K*  I 7-M  ♦ 1 

cmcc  r 

1 E’2 

JEf  * .LT.  1 1 GC  TC  300 

SHE  p7 

1533 

!*«*  .50.  TRlC-r  TO  ?0C 

SHEET 

1334 

JFCK  .GT.  IF  ICC  TO  350 

SHEET 

1E3- 

MM  s {P-K4  1 

SHE  at 

l 535 

A'K  ,M|«A(K,Ml4A(K. MM  1 /CONST  A 

sheet 

1*37 

°(«  I *“<*  l-A(  v,w»*cCNSTC 

SHAFT 

15  3** 

ACS  .MM  1*0. 

CMC  E ▼ 

15  30 

GC  TO  300 

c me  cT 

1540 

350 

CONTINUE 

Shact 

154  1 

r 

CMC  A T 

1 542 

MM  SK-!  F4 1 

ShF'T 

1 54-* 

A(K,M)*MX  .MUMIC.YM)  /C  ON  a T A 

CMC  p Y 

1544 

«<*C|*ef<!-A<!e,  MM  I * CON  ST  G 

SMTPT 

1 54* 

AC  I F .MM  | *0. 

SHEET 

1545 

300 

CONT INUE 

SHEET 

1547 

C 

Shcct 

1545 

mb  I XMOANO-  1 

SHEET 

1 540 

OC  200  M*?.MF1 

SHEET 

1550 

MMxM 4 ( 

SHEET 

155! 

Af  IT. Ml* A f I7.M144C  IO.MMI/C0NST* 

SHAFT 

15*2 

I 77  *IF4M 

Shaet 

1553 

“f  !77»*PC  I 771-4  (1°, -MM)  *CC\’cTG 

SHEET 

1 554 

*00 

AC  ic, MM  1*0.0 

SHAFT 

15*5 

r 

sheet 

1556 

C 

SHEET 

15*7 

PPTIPK 

She  at 

1**5 

FNO 

137 


REFERENCES 


1.  Lee,  C.  H.,  and  Kobayashi,  S.,  "New  Solutions  to  Rigid-plastic  Deformation 
Problems  Using  a Matrix  Method,"  Trans.  ASME,  J.  of  Engrg.  for  Ind., 

Vol . 95,  p.  865,  1973. 

2.  Lee,  C.  H.,  and  Kobayashi,  S.,  "Deformation  Mechanics  and  Workability  in 
Upsetting  Solid  Circular  Cylinders,"  Proc.  North  Amer.  Metalworking 
Res.  Conf.,  Hamilton,  Canada,  Vol.  1,  p.  185,  May  1973. 

3.  Shah,  S.  N.,  Lee,  C.  H.,  and  Kobayashi,  S.,  "Compression  of  Tall,  Circular, 
Solid  Cylinders  Between  Parallel  Flat  Dies,"  Proc.  Int.  Conf.  Prod.  Engr. , 
Tokyo,  p.  295,  1974. 

4.  Shah,  S.  N.,  and  Kobayashi,  S.,  "Rigid-plastic  Analysis  of  Cold  Heating  by 
the  Matrix  Method,"  Proc.  15th  Int.  MTDR  Conf.,  p.  603,  1974. 

5.  Lee,  S.  H.,  and  Kobayashi,  S.,  "Rigid-plastic  Analysis  of  Bore  Expanding 
and  Flange  Drawing  with  Anisotropic  Sheet  Metals  by  Matrix  Method," 

Proc.  15th  Int.  MTDR  Conf.,  p.  561,  1974. 

6.  Shah,  S.  N.,  and  Kobayashi,  S.,  "A  Theory  on  Metal  Flow  in  Axisymmetric 
Piercing  and  Extrusion,"  J.  Prod.  Engrg.,  Vol.  1,  p.  73,  1977. 

7.  Chen,  C.  C.,  Oh,  S.  I.,  and  Kobayashi,  S.,  "Ductile  Fracture  in  Axi- 
symmetric Extrusion  and  Drawing,"  USAF  Technical  Report  AFML-TR-77-96, 

June  1977. 

8.  Hill,  R.,  "A  Variational  Principle  of  Maximum  Plastic  Work  in  Classical 
Plasticity,"  Quart.  J.  Mech.  Appl.  Math.  1,  p.  18,  1948. 

9.  Prager,  W. , and  Hodge,  P.  G.  Jr.,  Theory  of  Perfectly  Plastic  Solids, 

Dover  Publications,  1951. 

10.  Hill,  R. , "On  the  Problem  of  Uniqueness  in  the  Theory  of  a Rigid-plastic 
Solid,"  Part  II,  J.  Mech.  Phys . Solids,  Vol.  5,  p.  1,  1956. 

11.  Miles,  P.,  "Bifurcation  in  Rigid-plastic  Materials  Under  Spherically 
Symmetric  Loading  Conditions,"  J.  Mech.  Phys.  Solids,  Vol.  17,  p.  303, 

1969. 

12.  Chakrabarty,  J.,  "On  Uniqueness  and  Stability  in  Rigid-plastic  Solids," 

Int.  J.  Mech.  Sci. , Vol.  11,  p.  723,  1969. 

13.  Hill,  R.,  Mathematical  Theory  of  Plasticity,  Oxford  Press,  1950. 

Hill,  R. , "Eigenmodel  Deformations  in  Elastic/Plastic  Continua," 

J.  Mech.  Phys.  Solids,  Vol.  15,  p.  371,  1967. 


138 


14. 


15.  Malvern,  L.t  Introduction  to  the  Mecnanics  of  a Continuous  Medium, 
Prentice-Hall,  1969. 


16.  Hill,  R. , "On  the  Problem  of  Uniqueness  in  the  Theory  of  a Rigid- 

plastic  Solid,"  Part  III,  J.  Mech.  Phys ■ Solids,  Vol.  5,  p.  133,  1957. 

17.  Hill,  R.,  "On  the  Problem  of  Uniqueness  in  the  Theory  of  a Rigid- 

plastic  Solid,"  Part  IV,  J.  Mech.  Phys.  Solids,  Vol.  5,  p.  302,  1957. 

18.  Hill,  R.,  "Stability  of  Rigid-plastic  Solids,"  J.  Mech.  Phys.  Solids, 

Vol.  6,  p.  1,  1957. 

19.  Strang,  G.,  and  Fox,  G.  J.,  An  Analysis  of  the  Finite-element  Method, 
Prentice-Hall,  1973. 

20.  Zienkiewicz,  0.  C.,  The  Finite-element  Method,  McGraw-Hill,  1971. 

21.  Odell,  E.  I.,  and  Clausen,  W.  E.,  "Numerical  Solution  of  a Deep  Drawing 
Problem,"  ASME  Paper  76-WA/prod-3,  1977. 

22.  Dahlquist,  G.,  Numerical  Methods,  Prentice-Hall,  1974. 

23.  Hill,  R.,  "A  Theory  of  the  Plastic  Bulging  of  a Metal  Diaphragm  by 
Lateral  Pressure,"  Phil.  Mag.,  Vol.  41,  p.  1135,  1950. 

24.  Woo,  D.  M. , "The  Analysis  of  Axisymmetric  Forming  of  Sheet  Metal  and 

the  Hydrostatic  Bulging  Process,"  Int.  J.  Mech.  Sci.,  Vol.  6,  p.  303,  1964. 

25.  Yamada,  Y.,  and  Yokouchi,  Y.,  "Elastic-plastic  Analysis  of  the  Hydraulic 
Bulge  Test  by  the  Membrane  Theory,"  Manuf.  Res.,  Vol.  21,  p.  26,  1969 
(in  Japanese) . 

26.  Wang,  N.  M.,  and  Shammamy,  M.  R.,  "On  the  Plastic  Bulging  of  a Circular 
Diaphragm  by  Hydrostatic  Pressure,"  J.  Mech.  Phys.  Solids,  Vol.  17, 

p.  43,  1969. 

27.  Shammamy,  M.  M. , and  Wang,  N.  M.,  "Comparison  of  Experimental  and  Theoreti- 
cal Results  for  the  Hydrostatic  Bulging  of  Circular  Sheets,"  SESA  Fall 
Meeting,  1970. 

28.  Iseki,  H.,  Jimma,  T.,  and  Murota,  T.,  "Finite-element  Method  of  Analysis 
of  the  Hydrostatic  Bulging  of  a Sheet  Metal,"  Bull.  JSME,  Vol.  17,  1974. 

29.  Wang,  N.  M. , "A  Variational  Method  for  Problems  of  Large  Plastic  Deforma- 
tion of  Metal  Sheets,"  General  Motors  Report,  1970. 

30.  Budiansky,  A.  B.,  "Nonlinear  Shell  Theory,"  J . of  Appl . Mech . , Vol.  35, 
p.  393,  1968. 

31.  Mel  lor,  P.  B.,  "Stretch  Forming  Under  Fluid  Pressure,"  J . Mech . Phys . 
Solids,  Vol.  5,  p.  41,  1956. 


139 


32.  Loxley,  E.  M. , and  Freeman,  P.,  "Some  Lubrication  Effects  in  Deep  Drawing 
Operations,"  J . of  Inst . Petr. , Vol . 40,  p.  299,  1954. 

33.  Keeler,  S.  P.,  and  Backofen,  W.  A.,  "Plastic  Instability  and  Fracture 

in  Sheets  Stretched  over  Rigid  Punches,"  Trans.  ASME,  Vol.  56,  p.  25,  1963. 

34.  Swift,  H.  W.,  "Plastic  Instability  Under  Plane  Stress,"  J.  Mech.  Phys . 
Solids,  Vol.  1,  p.  1,  1952. 

35.  Hill,  R. , "On  Discontinuous  Plastic  States  with  Special  Reference  to 
Localized  Necking  in  Thin  Sheets,"  J.  Mech.  Phys.  Solids,  Vol.  1,  p.  19, 
1952. 

36.  Marciniak,  Z.,  and  Kuczyr.ski,  K. , "Limit  Strains  in  the  Processes  of 
Stretch  Forming  Sheet  Metal,"  Int . J . Mech . Sci . , Vol.  9,  p.  609,  1967. 

37.  Marciniak,  Z.,  Kuczynski,  K.,  and  Pokara,  T.  , "Influence  of  the  Plastic- 
Properties  of  a Material  on  the  Forming  Limit  Diagram  for  Sheet  Metal 
in  Tension,"  Int.  J.  Mech.  Sci.,  Vol.  15,  p.  789,  1973. 

58.  Gosh,  A.  K.,  and  Hecker,  S.  S.,  "Failure  in  Thin  Sheets  Stretched  over 
Rigid  Punches,"  General  Motors  Report,  1974. 

39.  Kaftanoglu,  B.,  and  Alexander,  J.  M. , "On  Quasistatic  Axi symmetrical 
Stretch  Forming,"  Int.  J.  Mech.  Sci.,  Vol.  12,  p.  1065,  1970. 

40.  Chakrabarty,  J.,  "A  Theory  of  Stretch  Forming  over  Hemispherical  Punch 
Heads,”  Int,  J.  Mech.  Sci.,  Vol.  12,  p.  315,  1970. 

41.  Woo,  D.  M. , "The  Stretch-forming  Test,"  The  Engineer,  Vol.  220,  p.  876, 
1965. 

42.  Wang,  N.  M. , and  Gordon,  W.  J.,  "On  the  Stretching  of  a Circular  Sheet 
by  a Hemispherical  Punch,”  General  Motors  Report,  1968. 

43.  Wang,  N.  M.,  "Large  Plastic  Deformation  of  a Circular  Sheet  Caused  by 
Punch  Stretching,"  General  Motors  Report,  1969,  and  J,  Appl . Mech., 
p.  431,  1970. 

44.  Wifi,  A.  S.,  "An  Incremental  Complete  Solution  of  the  Stretch-forming 
and  Deep  Drawing  of  a Circular  Blank  Using  a Hemispherical  Punch," 

Int.  J.  Mech.  Sci.,  Vol.  18,  p.  23,  1976. 

45.  Lee,  C.  H.,  Masaki,  S.,  and  Kobayashi,  S.,  "Analysis  of  Ball  Indentation," 
Int.  J.  Mech.  Sci.,  Vol.  14,  p.  417,  1972. 

46.  Jarvinen,  P.  A.,  "Representation  of  High  Temperature  Plastic  Behavior 
of  Austenitic  and  Ferritic  Stainless  Steel  by  Empirical  Equations," 

Scand.  J.  Metallurgy,  Vol.  6,  1977. 

47.  Voce,  E.,  J.  of  Inst,  of  Metals,  Vol.  74,  pp.  537-562,  1948. 

48.  Gegel,  H.  L.,  Private  communication,  January  1978. 


140 


F 


49.  Chakrabarty  ' <d  Mellor,  P.  B.,  "A  New  Approach  for  Predicting  the 

Limiting  Lh  i<atio,"  IDDRG  5th  Biennial  Congress,  September  30, 

1968,  la  Metaxiurgia  Italiana,  p.  791,  1968. 

50.  El-Sabaie,  M.  G. , and  Mellor,  P.  B.,  "Plastic  Instability  Conditions 
in  the  Deep  Drawing  of  a Circular  Blank  of  Sheet  Metal,"  Int.  J,  Mech. 

Sci ■ , Vol . 4,  p.  535,  1972. 

51.  Budiansky,  B.,  and  Wang,  N.  M. , "On  the  Swift  Cup  Test,"  J . Mech . Phys . 
Solids,  Vol.  14,  p.  357,  1966. 

52.  Chung,  S.  Y.,  and  Swift,  H.  W.,  "Cup  Drawing  from  a Flat  Blank,"  Proc . 
Instr.  Mech.  Engrs.,  Vol.  165,  p.  199,  1951. 

53.  Woo,  D.  M. , "Analysis  of  the  Cup  Drawing  Processes,"  J.  Mech.  Engrs., 

Vol.  6,  p.  116,  1964. 

54.  Lee,  P.  K.,  Choi,  C.  Y.,  and  Hsu,  T.  C.,  "Effect  of  Drawing  on  Formabil- 
ity  in  Axi symmetrical  Sheet  Metal  Forming,"  J . of  Eng.  Ind . , p.  925,  1973. 

55.  Hsu,  T.  C.,  Dowle,  W.  R.,  Choi,  C.  Y.,  and  Lee,  P.  K.,  "Strain  Histories 
and  Strain  Distributions  in  a Cup  Drawing  Operation,"  J.  of  Eng.  Ind, , 
p.  461,  1971. 

56.  El-Sebaie,  M.  G.,  and  Mellor,  P.  B.,  "Double  Operation  Deep  Drawing," 

Int.  J.  Mech.  Sci.,  Vol.  15,  p.  945,  1973. 

57.  Chiang,  D.  C.,  and  Kobayashi,  Shiro,  "The  Effect  of  Anisotropy  and  Work- 
hardening Characteristics  on  the  Stress  and  Strain  Distribution  in 

Deep  Drawing,"  ASME  paper  66-prod- 3,  1966. 

58.  Woo,  D.  M.,  "On  the  Complete  Solution  of  the  Deep  Drawing  Problem," 

Int.  J.  Mech.  Sci.,  Vol.  10,  p.  83,  1968. 

59.  Levy,  S.,  Shih,  C.  F.,  Wilkinson,  J.  P.  D.,  Stine,  P.,  and  McWilson, 

R.  C.,  "Analysis  of  Sheet  Metal  Forming  to  Axisymmetric  Shapes," 

General  Electric  Report  No.  77CRD257,  December  1977,  Schenectady, 

New  York. 


141 


<*U.S.Qov«rnm«nt  Printing  Offlc#:  1979  — 657*002/490 


